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LXVII. The Production of Hyperons and Heavy Mesons 


By G. D. Jamzs and R. A. SatmERon* 
The Physical Laboratories, The University, Manchestert 


[Received February 21, 1955] 


ABSTRACT 
A full dynamical analysis of the reaction 


is given and it is shown that some of the results apply to the production 
of hyperons in any collision with a free nucleon at rest. In particular, 
difficulty is found in accounting for the low-energy hyperons that are 
frequently observed in both cloud chambers and photographic emulsion. 
Their existence necessitates a mechanism by which hyperons can lose 
most of their energy. It is suggested that this process is scattering by 
nucleons. 

19 cloud chamber photographs showing two V-events associated with 
the same nuclear interaction are discussed. In no case do the V-events 
represent the decays of the two secondaries of reaction (1). 

The observed frequency of associated V-events provides no evidence 
for or against the hypothesis that heavy unstable particles are always 
produced in reaction (1) or in similar ‘double production ’ reactions. 
Because the secondary particles from these reactions are scattered in 
heavy nuclei it is concluded that the reactions can best be studied in 
light materials. 


$1. INTRODUCTION 


In this paper we discuss the production of hyperons and heavy mesons 
considering, in particular, the evidence available from a cosmic ray cloud 
chamber experiment. 

Fowler et al. (1953, 1954) first observed the double production of 
hyperons and heavy mesons in collisions between energetic negative 
z-mesons and free protons. Several events observed in the gas of a 
hydrogen-filled diffusion cloud chamber were interpreted as examples of 


the reaction oe PY EK 


where Y stands for a hyperon and K for a heavy meson.{ Since their 
discovery of this type of reaction a number of examples of the production 
TEESE 

* On leave from the Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro. 


+ Communicated by G. D. Rochester. | 
+ Throughout this paper we use the nomenclature proposed by Amaldi et al. 
fo) 


(1954). 
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of Y- and K-particles in a single nuclear interaction have been reported. 
These are considered briefly in § 2. 

There is a good theoretical reason for believing that Y- and K-particles 
should be more readily produced in pairs than singly (Nambu ef al. 
1951, Pais 1952). Their long lifetimes cannot be reconciled with their 
copious production if they are produced singly. In fact, it may well 
be that every Y- or K-particle observed experimentally is produced 
in association with another such particle which generally escapes detection. 
The experimental evidence for and against this hypothesis is discussed 
in the later sections of the paper. 

When it is realized that some thousands of unstable heavy particles 
have now been observed in cosmic ray cloud chamber experiments, 
it appears surprising that so little can be said with confidence about 
the mechanism of their production. The reasons for this situation are 
considered in § 5. 


§2. PUBLISHED EXAMPLES OF ASSOCIATED PRODUCTION 


The most important evidence for double production of hyperons and 
heavy mesons comes from the work with the Brookhaven accelerator 
summarized by Fowler et al. (1955). In interactions between negative 
m-mesons and protons in hydrogen nine events have been observed that 
are interpreted as examples of the reaction 


a-+P+Y-+K. 


Only in a few cases could complete dynamical measurements be made 
and there is some uncertainty in the masses that must be assigned to the 
secondary particles. In particular, Fowler et al. show that in some 
events the neutral hyperon produced is a A°-particle while in others 
the measurements are better fitted by assuming a larger mass for the 
Y°-particle. It is suggested that a A°-particle could subsequently arise 
from the rapid y-decay of the heavier hyperon. 

Details have also been published of a further four events that are 
dynamically compatible with the general reaction 


a+ N->Y-+K PNET Evie att 


where N represents a nucleon. In each of the four cases the reaction 
occurs in a complex nucleus. One of the events (Fowler et al. 1954) 
showed the decays of a A°- and a 6°-particle coming from a nuclear 
interaction in the wall of a cloud chamber, A similar event was reported 
by Thompson et al. (1954) but in this case the nuclear interaction was 
caused by a cosmic ray particle. The remaining two events were found 
by Dahanayake et al. (1954) in nuclear emulsion. These events are 
stars from which a charged hyperon and heavy meson emerge; measure- 
ments on both events show that they are compatible with the dynamics 
of reaction (1). 
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In addition to the above well-measured examples of double production 
in reaction (1) there are many examples of associated* Y- and K-particles 
coming from a single nuclear interaction. In nuclear emulsion experi- 
ments such association particles have been described by Lal et al. (1952) 
and by Debenedetti et al. (1954). In cloud chamber experiments associated 
pairs of V-events are quite frequently observed but the only detailed 
analysis of such pairs has been made by Ballam ef al. (1955). In this 
paper we report measurements on 19 pairs of associated V-events. 

In cosmic ray experiments on heavy unstable particles production 
occurs in complex nuclei where scattering might be expected to mask the 
dynamical relations of the production process. The four examples of 
reaction (1) that have been observed in heavy nuclei show that this is 
not always so and, for this reason, we give in §3 a full analysis of the 
dynamics of reaction (1). In later sections we compare our experimental 
observations with the results of this analysis. 


§3. Dynamics oF DouBLE PropucTION PROCESSES 
(a) The Interaction 7-+-N+Y+K 
Some features of the dynamics of reaction (1) have been described 
by Reynolds and Treiman (1954) and by Dahanayake et al. (1954). We 
wish to have a complete description of the reaction in both the laboratory 
and centre of mass systems of coordinates. 


Fig. 1 


(a) Cb) 


Diagram explaining the nomenclature used in the dynamical analysis of the 
~~ reaction n+N— Y+K. (a) represents the collision in the laboratory 
system where the nucleon is at rest and ()) shows the collision in the 


centre of mass system. 


istingui i nd ‘ associated ’ production 
* We distinguish between double production and * associated P sapietis 
restricting the former to describe production in a specific reaction such as (1). 
By associated particles we mean only particles which appear to come from a 
single nuclear interaction. 


Z2R2 
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Figures 1 (a,b) show the quantities used to describe the reaction 
in the two frames of reference. In the laboratory system (in which the 
nucleon is assumed to be at rest) there are three momenta, p,, py. and 
p,, and two angles, dy and ¢x. These quantities are, of course, related 
by the vector equation py+px=p,. In the centre of mass system the 
momentum (p*) of either secondary particle is determined uniquely by p, 
and, for a given p,, the only parameter needed to describe the reaction 
is @*, the angle between the trajectories of the incident 7-meson and 
the emergent Y-particle. 

Figure 2 (a) shows the variation of dy with 1/py for different values 
of p, and 6* while figure 2 (b) shows the corresponding variation of dx 
with 1/p,x. In calculating these values the masses of the hyperon and 
heavy meson have been taken to be those of the neutral A°- and 6°-particles 
(2181 m, and 965 m,). 

Figure 2 shows some features which demand particular attention. 
Most striking are the limits within which py is confined. Under no 
circumstances can the Y-particle have a momentum of less than 195 Mev/c 
in the laboratory system nor can its trajectory ever make an angle 
greater than 57° with that of the incident 7-meson. For low values of 
p, the limits imposed on py and ¢y are even more restrictive. 

A second feature of the dynamics which is apparent from fig. 2 is that 
the Y- and K-particles produced in reaction (1) cannot have a small 
angle between their trajectories unless one of their momenta is small. 
This is better shown in the type of diagram used by Dahanayake et al. 
where the values of py and px are plotted for a given value of the total 
angle ¢(=¢y+¢x). This diagram is drawn in fig. 3 for the interesting 
cases where ¢ is small. It shows, for example, that if ¢ is less than 30° 
either py or px must be less than 6 x 108 ev/c. 


(b) More Complex Production Processes 


The limitation of py shown in fig. 2 (a) is due to the difference in mass 
between the target nucleon and the secondary hyperon. It arises because 
increasing the value of p, never makes the velocity of the hyperon in 
the centre of mass system equal the velocity of the centre of mass system 
in the laboratory. Letting yy and y~y correspond to these two velocities 
we find that 

Yy/Yom—>™My/My aS p,—->00 
where my and my are the nucleon and hyperon masses. To find the 
minimum value of py we take the difference between these two velocities 
which leads to the result 


Py > kmy(my/my—My/My). 
This quantity is equal to 195 Mev/c if my is taken to be 2181 m,. 
In this extreme relativistic approximation the masses of the z- and 


K-mesons do not enter the calculation of the hyperon momentum and 
the limiting curve in fig. 2 (a) holds for all reactions of the type 


A+N>Y-+B Fe ee at ee fey 
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(b) S 
Momenta and angles of emission of the secondaries from the reaction 
at+tN>+Y+K. The full lines are curves of constant Pas the broken 
lines are curves of constant 6* (see fig. 1). The value of p, is given in 
Gev/c against each solid line and the values of 6* are given in degrees. 
(a) shows the variation of dy with I/py and (b) the corresponding 
values of ¢x and 1/px. The masses of the Y- and K-particles have been 
taken to be 2181 m, and 965 m,. py 
The solid lines marked E.R. show the limiting angles and momenta 
as px ~ ©. These lines represent the limits for any reaction of the type 


PEN BO. oe. (see text). 
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whatever the masses of the particles A and B. A full discussion of 
reaction (2) for all values of p, is complicated but it is not difficult to 
show that the minimum value of py occurs when p, has an extreme value 
—either the threshold value or infinity. The minimum value of py is 
then either its value at the threshold of the reaction or 195 Mev/c. 


Fig. 3 Fig. 4 


° 2 4 6 8 10 12 | 
P [108 ¢) : p, [1% ewe | 
Fig. 3 
The variation of px with py for constant ¢(=¢x+¢y) in the reaction 
a+N-—-Y-+K. For values of py > 1-15 Gev/c the angle is never less 
than 30°. 
Fig. 4 
Differential momentum spectrum of A°-particles at production. The circles 
and full line are obtained from our data ; the triangles and interrupted 
line from the measurements of Gayther and Butler (1955). The two 
curves are made to coincide at a momentum of 0-7 Gev/c. The total 
number of events used to obtain this spectrum is only 47. The full line 
corresponds to a power-law spectrum with an exponent of —2-0, 


A further generalization is possible to reactions of the type 
A+-N-Y + Bs-C52 D295) sete enn y 
where A, B, C, D, are arbitrary particles. For a given momentum of 


particle A the maximum value of p* for the Y-particle occurs when the 
other secondaries (B, C, D,...) have the same velocity ; that is, they 
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behave as a single particle of mass m,-+-mp+mo+.... Since we are 
interested in the maximum value of p* when finding the limiting values 
of dy and py, reaction (3) is, for this purpose, equivalent to reaction (2).7 
This extension of the argument is of great practical importance since 
_an example has recently been reported of hyperon production in a collision 
between two nucleons (Block ef al. 1954). The reaction is interpreted as 


Nee NS Vee K oie OE ay 


In this reaction the momentum of the hyperon is restricted to the 
region inside the curve (K.R.) shown in fig. 2 (a). 

The production process discussed by Fowler et al. (1955) in which 
A°-particles result from the rapid y-decay of heavier hyperons (2) is 
also included as a reaction of type (3) since the reactions 


m+N>2°+0 and 3°+/1°+y 
are together equivalent dynamically to 
a7+N>A°+(+y 
provided that the 2°-particle decays before moving a detectable distance. 
Thus very slow hyperons cannot be produced even in this type of * cascade ’ 


process. 
(c) Effect of the Fermi Energy 


The results shown. in fig. 2 apply, strictly, in the nucleon frame of 
reference. If account is taken of the Fermi energy of a bound nucleon 
a further Lorentz transformation is required to find the angles and 
momenta in the laboratory system. This has an appreciable effect 
on the minimum value of py. If a nucleon is assumed to have a kinetic 
energy of 10 Mev, it is possible to produce a Y-particle whose momentum 
in the laboratory’ system is only 30 Mev/c (corresponding to a kinetic 
energy of 0-5 Mev). However, the chance of producing such a low-energy 
Y-particle is extremely small since it demands a head-on collision between 
a very energetic 7-meson and a nucleon of 10 Mev energy followed by 
the emission of the Y-particle at an angle 6* greater than 177°. 


§4. SumMaRy OF EXPERIMENTAL Data 


The observations that are reported here were all made with a large 
magnet cloud chamber situated at the Jungfraujoch (3580 m). The 
apparatus has been described by Newth (1954) and a typical system of 
counter control is shown in the same paper. During the past four years 
some 70 000 photographs have been taken and about 40° of these show 


+ We are extremely grateful to Mr. B. Powell for making this point clear 
to us. The argument can be outlined as follows. Let the secondary particles 
B, G, D, ete. (Y excluded) have a total energy U* and total pigeon p* in 
the centre of mass system. The quantity U ae ie is equal to U®, where U is 
the total energy of the particles B, C, D, etc. in their centre of mass system. 
Thus, if U* is fixed, the maximum value of p* is obtained when U is a minimum, 
This leads directly to the required result. 
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evidence of nuclear interactions in the lead transition layer above the 
cloud chamber. Most of the photographs were taken with the chamber 
free from solid material but on occasions a 3cm lead plate (or, more 
recently, a 1-2cm copper plate) was placed across the centre of the 
chamber. 

In table 1 a summary is given of the number of V-events observed. 
The V°-events include both A°-decays and 6°-decays and the V+-events 
include, presumably, decays of charged hyperons and heavy mesons. 
Table 1 gives also the number of photographs which show two associated 
V-events. 

Table 1. Yield of Single and Associated V-Events in 
Jungfraujoch Cloud Chamber | 


No. of No. of 
associated | associated 


yo V= | 


Site of nuclear | No. of single No. of single 
interaction | neutral V-events | charged V-events yo. po 


Transition 
layer 106 13 


Inside the 
chamber 


Table showing the yield of V-events in the Jungfraujoch cloud chamber. The 
events originating in the 25 cm lead transition layer above the cloud 
chamber are distinguished from those originating inside the cloud 
chamber. The latter group includes three neutral V-events coming 
from stars in the chamber gas. The event marked with an asterisk 
is not discussed further since one of the V-events appears to be neither 
a A°. nor a 6°-decay. 

A distinction is drawn between nuclear interactions occurring in 
the transition layer and inside the chamber, since the probability of 
observing the decay of an unstable particle depends strongly upon the 
distance of its origin from the illuminated volume of the chamber. Among 
the particles produced inside the cloud chamber there are three V°-particles 
produced in stars in the chamber gas: the remaining particles were 
produced in the lead or copper plates. 

Most of the V-events observed in the Jungfraujoch experiment cannot 
be classified as decays of a particular type. This is because the momenta 
of many of the particles involved in the decays are high and the maximum 
detectable momentum in the chamber is only 4 Gev/c. Newth and 
James (1953) made a statistical survey of the high-energy neutral 
V-particles observed and showed that approximately equal numbers of 
A°- and 6°-particles decayed in the chamber. From this they showed 
that comparable numbers of the two types of particle were produced in 
nuclear interactions above the chamber. In itself, this result is suggestive 
of some common production process—in fact, it was the starting point 
of the present work. 
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Concerning the charged V-decays, it is only possible to say that they 
include no events that cannot be interpreted as K-meson decays. The 
majority of the events have been discussed by Buchanan et al. (1954) 
who found no evidence for hyperon decays. This does not necessarily 
mean that charged hyperon decays are rare; it may well be that few, 


if any, satisfy the rather stringent criteria used for classifying charged 
V-events. 
§5. ANALYSIS OF RESULTS 


(a) The Nature of the Associated V-Events 


The 19 associated pairs of V-events (table 1) merit particular study 
and the details of these events are set out in table 2. Where possible 
the neutral V-events have been classified from angle and momentum 


Table 2. Measurements of 19 Associated V-Events 


Vax Vo 
ree a d Pab Ba by 
erent nae Momentum Nature Momentum |(deg.) Bel deg.) | (deg.) | (deg.) 
pare (aev/c) ae (Gev/c) 
MJ29 | Not/A® >0-°8 ¢ Sal 15 9 |}46+5 |82+5 |538+5 
NJ42 Q >4 2 ~1-7 3 2 |86+5 |20+5 | 804 
NUI1 2 =} Q Sail 10 | 11 | 442 |11+45 |10+5 
OC31 2 q Not/A° ~3 35 8 |10+5 | 04+5 | 104 
OC56 2 ~3 2 ~3 10 | 12 |85+5 |80+10| 10-4 
PH161 | NotA? >4 Q >4 10 | 20 |85+2 |85+5 | 104 
QL759 2 ~5 @ >4 10 | 10 | 25+15|79+5 | 63-4 
RJ451 2 ~2 @ >2 10 6 | 57410] 53+10| 73+10 
RT851 Not/A°® Sal D Q 2 8 a 43 ae 2, 30 ele 5 13 = 
SF1841 @ >2 2 >2 Oo it — — — 
SG75 Vip. 0-7 Not/A® 2°2 21 0 |88+2 |45+10]48+ 
SH232 2 ~2 NotA° ~3 8 | 10 |89+5 |14+5 | 8942 
SL325 2 Sl Not/A° ~2 8 2 | 31410] 55+10| 45+ 
RC941 | NotA® ~2 Me ~1:3 15 | 10 | 33+5 SY +5 | 294 
RK829 Q ~0:8 7+ ae e 2145 |27+5 |45+5 
RQ31 | Not/® >2 Ve = 2 = == = 
: 2 ~) Vax >0:3 16-5| 3 |72+10/70+410| 
SJ388 2 
SL101 ¢ >2 V+ ~0-5 31 2 53 +10 37 +10} 86+ 
SN2069 Y 0-7 Vr >0:5 20 8 |61+5 |88+5 | 204 


The table gives details of the 19 pairs of associated V-events. The first 13 
are associated neutral V-events ; the remaining 6 are associated neutral 
and charged V-events. Columns 2 to 5 give information about the 
individual V-events. ¢ is the angle between the lines of flight of the 
two unstable particles. 2, is the number of associated fast charged 
particles seen in the cloud chamber. yap is the angle between the 
planes of decay of the two particles. %, and y are the angles that each 
plane of decay makes with the plane containing the lines of flight of 

- icles. 
Pl aoa of the V-particles have been obtained by the method 


described by Buchanan et al. (1954). 
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measurements as compatible or incompatible with A®- or 6°-decays. 
In every case a lower limit to the momentum of the unstable particle 
could be set. This information is given in columns 2 to 5 of table 2. 
The sixth column gives the angle (¢) between the trajectories of the 
unstable particles. 

From the values of ¢ and the momenta of the unstable particles, it 
is evident that not one of the 19 pairs is a simple example of double 
production in reaction (1). Figure 3 shows that for the small values 
of ¢ observed, one or other of the secondary particles should have a 
low momentum ; all the particles listed in table 2 have high momenta. 

Two other processes could account for these associated events. They 
could result from ‘ plural’ production, each being produced in a separate 
reaction inside the same nucleus. Little could then be said about the 
nature of the individual reactions. Alternatively, they could be examples 
of double production in reactions more complex than (1). In favour of 
the first interpretation is the high multiplicity of the nuclear interactions 
with which the events are associated. The number of associated shower 
particles seen in the cloud chamber is given for each case in column 7 
of table 2. These figures show clearly that our events are very different 
from the example of double production observed by Thompson ef al. 
(1954), where no additional products of the nuclear interaction were 
seen. Two of the events are reproduced in Plates 12 and 13. 

Nothing can be said with certainty about the second hypothesis 
that the events are examples of more complex reactions than (1). 
However, Ballam et al. (1955) have found an apparent correlation between 
the planes of decay of associated V-particles. The existence of such 
a correlation for our events would be a strong argument for their being 
examples of double production processes. 

In fact, we find no correlation of the decay planes. The last two 
columns in table 2 give the angles between the plane of decay of either 
V-particle and the plane containing the paths of the two V-particles 
{the angle less than 90° is quoted). Ballam ef al. found that, for the 
11 pairs of events they considered, these angles all satisfied the condition 
that one was more than twice as great as the other. Among our 17 
measurable pairs 9 satisfy this condition ; this is just what we should 
expect if the planes were randomly orientated. 

We conclude therefore that none of our 19 pairs of V-events shows 
the decays of the two products of reaction (1). Nor is there any evidence 
for the associated events resulting from a complex * double production ’ 
process. The most plausible interpretation is that each V-particle 
was produced in a separate reaction within a nuclear interaction of high 
multiplicity. Whether or not the individual reactions involved in this 
plural production are of type (1) we discuss below. 


(b) The Frequency of Double V-Events 


The significance of the fact that we do not observe any simple example 
of reaction (1) can only be assessed if we calculate how many examples 
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we might expect to observe if all Y- and K-particles were produced 
in this type of reaction. An exact calculation of the experimental bias 
involved is very difficult, but we give in this section the results of some 
approximate calculations that are sufficient to indicate the orders of 
magnitude involved. 

Neutral V-decays are much more efficiently observed than are the 
decays of charged K-mesons. This is because the mean lifetimes of the 
AA°- and 6°-particles are such that a large fraction of the particles decay 
within a distance of 10 to 100 cm from their origins (assuming momenta 
-of a few Gev/c). The charged K-mesons have much longer lifetimes and 
relatively few decay in even the largest cloud chambers. For the charged 
hyperons the situation is still obscure since the number of identified 
‘decays in cloud chambers is very small indeed. 

The relatively high probability of observing neutral V-decays suggests 
the following approach to the calculation of the expected number of 
-double events. We consider only the reactions 


Sa N10 0 ae ee ee 17,3 
and arte N llores ee eee eee ur (Lb) 


For reaction (la) we find the probability that, having observed one 
V°-decay, the other should be seen. For (1b) we find the probability 
‘that, having seen the decay of the K-meson, the /°-decay is also seen. 
These two probabilities both refer to chance of seeing a neutral V-event 
-and can be calculated from the dynamics of reaction (1) using the known 
properties of the A°- and @°-particles. The first probability we call 
N(2V°)/N(V°), and the second N(V°-+K)/N(V°). 

Values of these two probabilities have been calculated for various 
parameters of reaction (1) and the results are shown in table 3. Figures 
_are given both for production in the transition material above the chamber 
and for production in a plate inside the chamber. The advantage to 
be gained from studying nuclear interactions inside the chamber is 
obvious. 

From table 3 the values of N(2V°)/N(V°) for production above and 
inside the chamber are of the order of 2°, and 20°, respectively. This 
suggests that 2°% of the observed V°-events from reaction (1a) occurring 
in the transition layer should be accompanied by a second visible V°-event. 
Of the 750 single V°-events observed we might expect one quarter to be 
‘produced in reaction (la), the remainder being produced in (16) or in the 
reactions which result in a charged hyperon. The absence of any pair 
-of V°-events compatible dynamically with (1a) thus suggests either that 
the reaction is not as important as we have assumed, or that the secondary 
particles produced suffer scattering sufficient to destroy the dynamical 
relations of (1). A similar conclusion is reached from considering the 
-other probabilities given in table 3 although it should be emphasized 
that the statistical weight of all the observations is poor and the possibility 
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of A°- and 6°-particles decaying into neutral secondaries has not been 
considered. 


Table 3. Expected Frequency of Associated V-Events from 
a+N-Y+K 


Interaction Site of interaction 


parameters Transition layer Inside chamber 


N(2V/N(V% | N(V°+K)/N(K) | N(2V9/N(V2) |N(V°+K)/N(K) | 
(Gev/c) (%) (%) (%) (%) 


1 0-2 40 40 
vo 


1 to 3 5 16 
(1 0-2 
0°2 1tol-5 0 
1:5to3 0 


‘al 
180° 
iis to3 


The calculated probabilities of seeing double events produced in the interaction 
a+N—+Y-+K for various parameters of the reaction. The table 
gives a probability of seeing the second neutral V-event when one 
neutral V-event has been seen and also the probability of seeing a 
V°-decay when a K+-decay has been seen. The probabilities are given 
as percentages. 


(c) Slow Neutral Hyperons 

The importance of observations of slow hyperons has been pointed 
out by Reynolds and Treiman (1954). From the discussion in § 3, it is 
clear that slow hyperons can only arise from reactions (1) or (2) if the 
hyperon is emitted backwards in the centre of mass system. 

Among our events there are three decays of remarkably low-energy 
A°-particles. One of these has previously been described by Millar and 
Page (1953). The A°-particle had a momentum of (190+15) Mev/c. 
The other two events are decays of A°-particles whose momenta are 
(230+40) and (150+20) Mev/c. Reference to figure 2 (a) shows that 
the conditions necessary to produce such slow A°-particles in reaction (1) 
are highly improbable. Moreover, the particle of momentum 
(230+40) Mev/c comes from an interaction in the copper plate across the 
middle of the chamber and is emitted at an angle of 90° to the path of the 
primary particle causing the interaction. This makes interpretations 
of the event as a simple example of reaction (1) almost impossible. 

Any or all of these events can, of course, be accounted for by assuming 
production according to reaction (1) with subsequent loss of energy by- 
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scattering of the A°-particle. To obtain more information on this point 
we have estimated the momentum spectrum of the A°-particles produced 
in our experiment. Confining attention to particles with near-vertical 
trajectories, the momenta of a random sample of A°-particles were obtained 
and, knowing the mass and mean lifetimes of the particles, it was possible 
to calculate the efficiency of detecting particles of any momentum and 
hence to find the differential momentum spectrum at production. This 
spectrum is shown in fig. 4; to our results have been added the data of 
Gayther and Butler (1955) referring to particles of very low momentum. 

Figure 4 shows that the momentum spectrum of the A°-particles is 
similar to that of the charged shower particles in penetrating showers 
(Barker and Butler 1951). While the observations contain a large 
statistical uncertainty, there is good evidence that the spectrum rises 
with decreasing momentum down to the lowest momenta accessible to 
observation. This is strong confirmation of the results of Ballam et al. 
(1953) and is in complete disagreement with what would be expected from 
reaction (1) with isotropic emission of the secondary particles in the 
centre of mass system. The significance of this result is considered in 
§6 where the observations of slow A°-particles in photographic emulsion 
are discussed. 


§ 6. Discussion 
(a) Hyperon Production 


The most interesting and important point that arises from the foregoing 
sections is that it is extremely difficult to produce low-energy hyperons 
in any type of collision with a nucleon at rest. There is, however, strong 
evidence for the production in nuclear reactions of a large number of 
slow A°-particles. In addition to our evidence and that of Gayther and 
Butler which has already been quoted, the results of the Princeton group 
(Ballam et al. 1953) also show a large fraction of low-energy A°-particles. 

Photographic emulsion experiments have produced the best evidence 
for the existence of very slow A°-particles. Danysz and Pniewski 
(1953 a, 1953) first observed a nuclear fragment that disintegrated 
after coming to rest in the emulsion. They interpreted this as a fragment 
in which a bound A°-particle took the place of a neutron. Since their 
discovery many similar events have been observed (Padua Conference 
Report 1955). The binding energy of the A°-particles in these fragments 
cannot be accurately found but is certainly not more than a few Mey. 
These bound particles are therefore good evidence for the formation of 
very slow hyperons in nuclear interactions. 

If every hyperon is produced in a collision with a single nucleon at rest 
the observations of slow hyperons necessitate some mechanism of energy 
loss. The simplest mechanism is scattering of the hyperons by the 
nucleons in the parent nucleus. (In a single elastic scatter nearly all 
the hyperon energy can be lost. The maximum angle of scatter is, 
from §3, 57°.) By assuming that such scattering usually occurs we can 
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explain all the experimental observations of slow hyperons, the existence 
of bound A°-particles, and the results of Gayther and Butler (1955) who 
find that there is a considerable fraction of backward-moving slow 
A°-particles from nuclear interactions in lead. 

The four examples of reaction (1) that have been observed in heavy 
nuclei (see § 2) are clearly cases where neither the Y- nor the K-particle 
has been appreciably scattered after production. If we assume that 
both Y- and K-particles have a mean path 2 free from scattering in 
nuclear matter, the probability of both emerging from a nucleus without 
scattering is exp —(ly+/,x)/A where /y and /, are their paths in the parent 
nucleus. 

Lock and Yekutieli (1952) have discussed the analogous problem of 
the scattering of recoil protons in nuclear interactions. They found 
that their experimental observations were fitted by //A~2. If the heavy 
unstable particles are scattered by the same process as protons we should 
expect //\~2 for them also (neglecting differences in nuclear radii). 
Thus the probability of both particles emerging without scattering from 
a heavy nucleus would be about exp (—4) or ~2°%. In this way we 
can account for the few recorded examples of reaction (1) occurring in 
heavy nuclei. 


(6) The Frequency of Double V-Events 


It was shown in §5(b) that in our experiment several examples of 
double V-events compatible with (1) should have been observed if all 
Y- and K-particles were produced in this type of reaction. The calculation 
ignored the possibility of the secondary particles being scattered. If 
we assume that only in about 2°, of examples of (1) do the secondaries 
emerge without scattering our failure to observe such an example is not 
surprising. 

In general, the hypothesis that scattering nearly always destroys 
the dynamical relations of the production processes makes the study of 
these processes in heavy nuclei almost impossible. Only a study of 
nuclear interactions in light materials could give a high proportion of 
valuable dynamical information. 


(c) Anisotropy in the Reaction 7+N>Y+K 

To account for the large number of slow A°-particles observed in the 
Princeton experiments Reynolds and Treiman (1954) suggested that 
there must be a strong anisotropy in the emission of hyperons in reaction 
(1). If the hypothesis of hyperon scattering outlined above is correct 
it is unnecessary to assume anisotropy to account for the emission of slow 
hyperons. 

The direct evidence available from completely measured examples of 
reaction (1) suggests very strongly that anisotropy favouring values of 
é* near to 180° does exist. Table 3 shows that there may be some 
experimental bias against observing examples of reaction (1) with small 
values of 6*. It seems unwise, at present, to combine clear examples 


Production of Hyperons and H eavy Mesons 585 


of reaction (1) in which scattering is known to be absent with incomplete: 
observations where only py is known. 


(d) Polarization Effects in Associated V-Events 


Since our associated V-events are interpreted as the products of 
separate reactions the absence of any angular correlation between their: 
planes is not surprising. What is difficult to explain is the difference 
between our results and those of Ballam et al. (1954). Among the 11 
events that these workers considered there were 7 pairs of V°-events 
obtained in a cloud chamber experiment very similar to our own. Ballam 
et al. do not publish sufficient information to show whether or not their 
events could be examples of reaction (1) but, from our experience, it 
is unlikely that they are. If their 7 examples are added to our 12 mea- 
surable pairs of V°-events we find that 13 pairs satisfy the condition 
$1 > 2, while 6 do not. If the planes of decay were randomly orientated 
we should expect half the pairs to satisfy this condition. 

It seems again preferable to discriminate between double events 
that are dynamically compatible with reaction (1) and associated events 
that are not. At present, the evidence for anisotropic distribution 
of the decay planes is meagre for both categories. 


§7. CONCLUSIONS 

We can summarize the conclusions drawn from this study in the 
following way. 

(i) To account for the copious production of slow hyperons in nuclear 
interactions some mechanism of energy loss is essential. Scattering 
of hyperons by nucleons is the simplest mechanism. 

(ii) If the hyperon—nucleon scattering cross section is comparable with 
the nucleon-nucleon scattering cross section, the observed features of 
hyperon production can all be accounted for without difficulty. 

(iii) The study of production reactions such as 7+-N->Y-+-K in heavy 
nuclei is extremely inefficient since scattering of the secondary particles 
will usually destroy the dynamical relations of the reaction. The 
anisotropy of emission in such reactions can only be studied in completely 
measured examples where scattering is known to be unimportant. 

(iv) The statistical analysis of the frequency of associated V-events 
in our cloud chamber does not provide evidence for or against the 
hypothesis that Y- and K-particles are always produced together. We 
interpret those associated V-events that we observe as examples of 
‘ plural’ production in separate reactions inside the same nucleus. We 
find no evidence for correlation between the two planes of decay. 
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Event PH 161 showing two neutral V-decays among a shower of twenty fast 
charged particles. The event is described in table 2. 
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Event RC 941 showing associated neutral and charged \-particle decays. 
The event is described in table 2. 
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LXVIII. Measurement of Ionization in Nuclear Emulsions 
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SUMMARY 


The gap-length distribution along the tracks of ionizing particles in 
nuclear emulsions has been investigated, and shown to be exponential for 
all values of specific ionization. The coefficient of the exponential, g, 
is suggested as the most useful measure of the ionization. In the first 
part of the paper, practical visual methods of measuring g are discussed. 
Accuracies not far inferior to the maximum possible are obtained by 
counting the small proportion of gaps of lengths exceeding two or three 
times the average value. The ratio of the coefficient g in a given track, 
to that in the track of a relativistic electron, is shown to be independent 
of degree of development, for particles below six times plateau ionization. 

In the second part of the paper, a detailed model is proposed, which is 
capable of explaining the dependence of g on the velocity and charge of 
the ionizing particle, and on the development and sensitivity of the 
emulsion. 

List of Symbols Employed 


Coefficient of the exponent of the gap-length distribution. 

Coefficient of the exponent for track at plateau ionization. 

9/9o- 

Blob density (i.e., total number of resolvable gaps per unit length). 

Plateau value of B. 

Parameter appearing in the empirical relation between 6 and g, 
i.e., B=g exp (—ga). « is nearly equal to uw! (see below). 

Minimum chosen gap length, as measured between inside edges of 
developed grains. 

Density of such gaps per unit length. 

Projected distance between centres of successive crystals along the 
path of the particle. 

Path length of ionizing particle through a crystal. 

Path length of particle through gelatin between adjacent crystals. 

k, 5, & Mean values of k, s, and x. 

I/k, Coefficient of the assumed exponential distribution in s. 

ui Mean distance between centres of two crystals, which after develop- 

ment can just be resolved in the microscope. 
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u(n) Distance between centres of two developed crystals containing n 
undeveloped crystals between them. 

p Probability of development of a particular crystal traversed by the 
ionizing particle. 

p Mean probability of development of all crystals traversed. 

A Number of acts per unit length in silver bromide. An ‘act’ is 
defined in § 13. 

g Probability that an act leads to development of the crystal; it 
depends only upon degree of development. 


§ 1. IyTRODUCTION 


Tue distribution in length of the gaps occurring in tracks of ionizing 
particles in nuclear emulsions, has formed the subject of a number of 
recent papers. This distribution may be specified by means of some 
suitable parameter, e.g. the mean gap length, the blob density or the grain 
density. It gives a measure of the velocity and charge of the ionizing 
particle producing the track. It was shown by O’Ceallaigh (1954) that, 
at least over a limited range of the ionization, the distribution in observed 
gap length was exponential and that this fact leads to a great simplification 
in the analysis of the experimental data. 

The present paper is divided into two parts. In the first part, observa- 
tions on the gap-length distribution are described and it is shown that 
this distribution is exponential over the entire range of ionization which we 
could measure (between once and several hundred times the minimum 
value). The coefficient, g, of this exponential appears to be the most 
useful measure of the ionization in the track. g has been determined as a 
function of residual range of protons and 7-mesons, and as a function of 
charge for relativistic particles. For values of ionization below about 
ten times minimum, it is shown that g is nearly proportional to the rate of 
energy loss of the ionizing particle in the silver bromide crystals which it 
traverses. Further, over this interval of ionization the value of g/g, 
where g, corresponds to minimum ionization, is practically independent of 
the degree of development of the emulsion. It is shown that g can be 
determined most rapidly from the ratio of the total number of observable 
gaps (i.e. the blob density), to the number of gaps of length greater than 
some optimum value. The value of g so obtained is independent of the 
mean developed grain size. The statistical errors associated with the 
measurement of g are discussed and the method is illustrated by a deter- 
mination of the mass of K-particles. In particular, the determination of 
masses by this method is shown to be at least as accurate as, and much 
more rapid than, the photoelectric method over the entire range of 
ionization considered. 

In the second part of the paper, a detailed model is proposed of the 
formation of developed crystals along a track. It is shown that this model 
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is able to account for the dependence of g upon both the specific ionization 
and the degree of development of the emulsion or the type of emulsion 
used. 


PAn De 
§ 2. EXPERIMENTAL ReEsutts on Gap Luneru DistRIBuTION 


_ Figure 1 shows the integral gap length distribution observed on the 
tracks of relativistic cosmic ray particles of different charges. Within the 
limits of experimental error, these distributions are exponential, and thus 


Fig. 1 


GAP LENGTH DISTRIBUTIONS on TRACKS 
of INDIVIDUAL RELATIVISTIC PARTICLES 


i 
Sees 270 


i 


1 
At gap |> ; 
tangth 
1S units 


fe) 1 2 3 4 5 6 y) 8 9 10 
. Gap length in scale divisions (div=-834) —> 
Gap length distribution in tracks of individual relativistic particles in emulsions 
of Stack 1. As described in the text, the length of each gap is measured 
between inside edges of the bounding grains. The number per unit 
length of all visible gaps corresponds to the blob density. Within 
experimental error, all the points lie on exponential curves. 


confirm previous measurements carried out by O’Ceallaigh on lightly 
ionizing particles. It follows that the ratio of the numbers H, and H, of 
fo} 
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gaps of length exceeding J, or /, per unit length of track, yields a measure 
of the coefficient of the exponent g : 


1 BS 
=< 2" _- Jog? ——= | ee ce 
a (le 10) e. (77 | 


The length of gap lis the distance between the inside edges of the developed 
grains bordering the gap. It should be noted here that we have supposed 
the gap length to be measured along the track. If however, the track 
dips at an angle @ to the plane of the emulsion, the projected gap length, J, 
as observed in the eyepiece scale of the microscope, is related to the true 
length J,, by 1,=l sec 6. Thus, the corrected value of g can be found from 
eqn. (la) simply by multiplying by a factor cos @ (O’Ceallaigh 1954). 

In general, we have chosen /, =0, corresponding to the smallest resolvable 
gap, so that H,=B, the blob density as customarily ee the density 
H of gaps exceeding length / is given by 


H=B exp (—gl). ee eo ote ee os LY 


Let g, be the value of g for the track of a high energy electron (corresponding 
to the plateau value of the ionization) in identical conditions of develop- 
ment ; and let g/g,=g*. Figure 2 shows the observed relation between the 
quantity g* and the residual range for protons in G5 emulsion. The 
measurements were carried out on protons and 7-mesons, the equivalent 
range of a proton being calculated assuming a mass for the 7-meson of 
273 m,. The observations were made in two separate stacks of stripped 
emulsion, each of approximately homogeneous development. In stack 1, 
the average value of gy=270/mm; in stack 2, only 180/mm. It will 
be seen that the values of g* obtained from the two stacks are practically 
indistinguishable down to proton ranges of the order of 1 em, but differ 
considerably at lower velocities. 

It is worth noting here that if instead of measuring the normalized 
coefficient, g*, as defined above, one chose to measure the probability of 
development, normalized to minimum, as described by other workers, the 
discrepancies between the two stacks would be much greater (~10% at 
1 cm range). 

The value of the exponent of the gap length distribution for relativistic 
heavy primaries is shown in fig. 3, plotted against Z?. It is clear that g 
tends to a limit as Z?— and hence the specific ionization of the particle— 
becomes very large. This limiting value of g, ~5000/mm, is the same for 
both stacks despite the large difference in development. In a detailed 
model to be described in Part II, we make the assumption that this value 
corresponds to the exponential distribution of gelatin gaps between 
the silver bromide crystals, when all these crystals have been developed. 

In fig. 4 are shown the experimental values of blob density B, plotted 
against the density H of gaps longer than 4:17 1. Each point corresponds 
to an observation on a single track in a particular emulsion of the stack. 
The number of blobs counted for each point was ~400, the number of 
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gaps ~60. The range of g covered by these measurements was from 
270—1200/mm. 


§ 3. ANALYSIS OF THE EXPERIMENTAL RESULTS 


With the help of the results of fig. 4, we find an empirical relation 
between B and g: 
B=g exp.(—ga) = Vo ee eee 
% is a parameter which is determined largely by the developed grain size, 
and to a lesser extent by the optical resolution of the microscope and the 
convention employed by a particular observer. 


Fig. 3 


COEFFICIENT of EXPONENT for 
RELATIVISTIC PRIMARY PARTICLES, CHARGE Ze 


' 2 3 5 10 2 3 5 100 2 3 s 


The measured values of the coefficient g, observed on individual relativistic 
particles, plotted against the square of the nuclear charge Z, as estimated 
from 6-ray counts. 


In fig. 4, curves, 1, 2, and 3 show the relation between B and H calculated 
from eqns. (1b) and (2) for values of “20-60, 0-65 and 0-70 « respectively. 
It will be seen that the shape of the calculated curves is in good agreement 
with the observations and the best fit is obtained for a—0-65 yu. The 
spread of the individual experimental points in this diagram is consistent 
with the statistical errors arising from the number of blobs and gaps 
counted (§ 5). 
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In Part I, with the help of detailed models of the distribution of crystals 
along the tracks of ionizing particles, we show that « can be identified with 
the mean separation of the centres of two crystals, which upon develop- 
ment are just resolved. Further, we show that eqn. (2) will be accurate 
to about 1% for all values of g<1200/mm. 


Fig. 4 
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: : h 
i lotted against the density of gaps exceeding 4:17 p. EKac 
aA Berle ee ponds a an individual track in a single emulsion, the mean 
number of blobs per point being ~400, of gaps ~60. Curves (1), (2) 
and (3) represent the empirical relation between B, g and H (eqn. 2) 


for «=0-60, 0-65, and 0-70 p respectively. 
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Thus in stack 1, from which the above measurements were taken, « is 
independent of both the coefficient g and of casual variations in develop- 
ment from one emulsion to another. «, however, can depend on the 
conditions of development, for a similar analysis carried out in stack 2 
gives a value of 0-74 for this quantity. 


§ 4, MEASUREMENT OF Low VALUES OF g AND NORMALIZATION 


In order to determine g when g is less than ~550/mm it is sufficient to 
measure the blob density B alone. The reasons for this are as follows. 
Blob counting is both reliable and rapid for lightly ionized tracks, since 
the proportion of developed grains which are close to each other is small. 
Further, any small variations in the developed grain size and/or optical 
conditions have, for the same reason, a trivial effect (§ 5(c)). The 
appropriate value of g has been calculated from B using eqn. (2) and the 
best value of «0-65 » given above. 

In order to determine g, throughout the stack, blob counts were made on 
relativistic electrons in cascade showers, in each emulsion, as a function of 
depth. The value of B, in a given emulsion varied by ~7°%, between top 
and bottom. Approximately 4000 blobs were counted in each emulsion of 
the stack. g, was calculated from B, and so determined to a statistical 
accuracy of less than 2%. No significant differences were found in the 
variation of B, with depth in the various emulsions. Furthermore, 
variations of B, over the surface area (150 cm?) of a given emulsion were 
found to be less than the statistical error of the blob counts (3°). 


§ 5. STATISTICAL ERRORS ASSOCIATED WITH DETERMINATION OF g 


In order to calculate the accuracy expected in individual determinations 
of the coefficient g, we quote here an approximate form of the result 
derived in the Appendix. Let o be the standard deviation associated 
with the number of gaps counted, V,,. Then 


f? l 


where (= 1/H —1/g—(l+-«), and is the mean distance between the centres 
of those crystals which upon development form the end of one gap and 
the beginning of the next. 

lis the length of the gap as previously defined, and « is the ratio of the 
number of gaps exceeding the minimum length / to the total number of 
developed grains in the track. The statistical error in a determination of 
blob density B can be obtained by substituting /=0 and B=H in (3). 

Expression (3) is sufficiently accurate in practice for all values of 
g and /, but the complete expression is given in the Appendix. 

As can be seen from (3), when the mean separation of the gaps f becomes 
very large compared with 1/g, then o\/(N)—1; the gaps are distributed 
at random. However, for the values of B and H commonly measured, 
o/(Njq) and oy/(N x) are significantly less than unity. o/(N,q) in fact 
passes through a minimum of ~ 0-52 when (/--a)g=1-0. 
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(a) Determination of g from B and H 


In order to determine the optimum value of the minimum gap length, 
a detailed calculation must be made. This must take into account not 
only the statistical fluctuations of B and H given by eqn. (3), but also 
the correlation between the fluctuations of H and those of B when taken 
over the same section of track. 

As a result, we find that the theoretical optimum value of J occurs 
when gl ~ 2-0, for all values of g. The accuracy does not vary appre- 
ciably in the interval 1-5<gl<2-5. The actual expression for the 
statistical error in g is somewhat complex, but is represented approxi- 
mately by l 


This expression is in fact a slight overestimate of the error, but. is 
sufficiently accurate in practice provided B/H>4, and N,>4N,,. It will 
be noted that the expression (4) is that which would have been obtained 
by assuming that o/(N)=1 always, and not as given by (3); that 
the correlations between B and H could be ignored; and that the 
statistical errors in B are negligible. 

Tn considering the optimum value of /g, we must also take into account 
the following practical considerations. In scrutinizing gaps of length 
close to the limiting value /, there will be a number of decisions to make. 
One cannot therefore exclude the possibility of a small temporary bias 
which might vary from track to track, and possibly be different for blobs 
and gaps. The proportion of such decisions is of course independent of 
the gap length chosen. The effect which the actual assignment of 
borderline gaps produces on the value of g obtained is clearly inversely 
proportional to/. Furthermore, the larger the value of / chosen the fewer 
the gaps that have to be counted and the more rapidly the measurements 
are completed. 

Thus we feel that the best value of / to use in practice is that which 
gives Ig ~ 2-5. 1 is of course chosen to be an integral number of eyepiece 


scale divisions. 


(b) Determination of g from B alone 


When the determination of g depends on a blob count alone, (§ 4), and 
if Nz is the number of blobs counted, then for low values of g, 


ov/(Nz) > BY 
Cee SPU) eee eee 
g V(Nzg) (l—ge) 
For such low values of g, as can be seen from (5), the fractional error in 
g is close to 1/\/(N,). For ga>0-35, the error in g rapidly becomes 


larger, and as mentioned previously, much more sensitive to real 
fluctuations in « and more dependent on decisions regarding borderline 


and 
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gaps. Hence one should not rely on blob counts alone for the deter- 
mination of values of gx greater than 0:35. We recommend that any 
tracks for which the dip correction exceeds the standard deviation in g, 
should be measured using B and H and not B alone, for all values of g. 


(c) Comparison of the Present Method with a Mean Gap-Length 
Determination 


In discussing the accuracy attainable by the present method, we must 
compare it with the maximum precision which it is possible to attain. 
The greatest precision would presumably be achieved by measuring the 
length of each individual gap along the track. Since the gap length 
distribution is exponential, the precision with which the coefficient g 
could be determined (e.g. by finding the mean gap length) would then be 
dg|g=1// (Np). 

On the other hand, the precision in g obtained by the method described 
in this paper is given with sufficient accuracy by eqn. (4). With 
log, B/H=2-5, we thus obtain as much as 55%, of the total information 
available in the track by counting (as opposed to measuring) the long 
gaps, amounting to only 8%, of the total, and by counting a sample of 
the blobs (such that the number of blobs counted is large compared with 
the number of long gaps). 


§ 6. SUMMARY OF PROCEDURE FOR DETERMINING g* 


It is now perhaps convenient to summarize the procedure that we 
have adopted for determinations of g*. 

(1) The plateau blob density B, must be measured in each emulsion, 
by counts on relativistic electrons. The dependence of By, on depth 
must be investigated and this is more conveniently achieved by gap 
counts on tracks of relativistic «-particles rather than on those of electrons, 
since more information is available per unit length of track. We have 
found that the depth dependence in a given batch of emulsions processed 
in a similar manner is quite uniform, but this point must be checked. 

(2) Measurement of the parameter « (or better, w!, as defined in 
Part II, eqn. (14)). « or w! can be most readily found by measurement 
of blobs and gaps on protons of 1 em range, where B is near its maximum 
value for average conditions of development. We find that « is nearly 
constant in a given batch of emulsions. It may differ slightly however 
from one observer to another. These two points must also be investigated. 

(3) From the values obtained for « (or wt) and By, gy is obtained 
from eqn. (2) at a given depth in a particular emulsion. 

(4) Measurement of g : 

(a) For lightly ionizing particles, when gx<0-35, g and hence g* is 
most conveniently determined from blob counts alone. 

(6) For values of gx>0-35 we have determined g from counts on blobs 
and gaps. The gap-length is chosen such that lg ~ 2-5 and the number 
of blobs counted is about four times the number of gaps. 
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(5) For blob counts alone, the fractional standard deviation in g is 


=1/V/(Nz). For blob and gap counts combined, the corresponding 
error is I 


— V(Nq) . log, B/H’ 


provided that the number of blobs counted is more than about four times 
the number of large gaps. In addition, normalization errors must be 
taken into account. As a check in estimating the accuracy in mass 
determinations on long tracks by this method, the final error in mass 
should also be obtained from the internal fluctuations of the individual 
measurements about the mean value. 


§ 7. DETERMINATION OF THE Mass oF INDIVIDUAL PARTICLES BY 
MEASUREMENT OF g AND REestpuAL RANGE 


Figure 2 shows that for values of the residual range R of a proton 
in the interval 1-40 cm, g*a R-°. Thus, neglecting errors in deter- 
mination of the range of a particle and straggling, the statistical error in 
mass M resulting from a determination of g* is 

dM Tage 
0g 

In order to illustrate the accuracy of the technique, we describe 
measurements on three heavy mesons observed in stack 1. Their 
ranges at emission were 9-1, 7-2 and 7-1cem. Some fifteen determinations 
of g* were made on each track, and the mean value M-°* was computed 
from the various mass values, suitably weighted according to the number 
of blobs and gaps counted and including statistical errors in normalization 
of individual emulsions. We obtained : 


M,=942+24m,; My=1022+30m,; Mz=1020+4 24 m,,. 


On coming to rest, particle X decayed into two 7+ and one 7 ~-meson, and 
wasa7zt. Particles Y and Z were K-mesons each giving rise to a singly 
charged relativistic particle, too steep for identification. 

The time required for the determination of the mass of each particle 
was less than two days for a single observer. This does not take into 
account time required for calibration of the emulsions. 

Methods so far described for the visual measurement of the gap 
distribution have effectively involved the estimation of the lengths of 
all individual gaps, and have proved very time-consuming. In the 
present method, the proportion of gaps which must be scrutinized care- 
fully is small, and for a given statistical accuracy the actual number of 
gaps counted is reduced by a factor of more than five. The increased 
speed of measurement that results appears to be one of the chief advantages 
of the present method. 

The photometric method gives a measure of the density of a track, 
which is substantially free from subjective errors. It has been chiefly 
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advocated on this account. The accuracy claimed for mass deter- 
minations on long tracks by the photoelectric method is no greater 
than that obtained by the technique we have described ; but the time 
required for these measurements is at least ten times as long. 


Parr 
§ 8. Basic AssuMPTIONS FOR MODELS OF GRAIN PRODUCTION 


The primary factor which determines the characteristics of a track is 
the number of silver bromide crystals developed. In this paper, we 
shall consider only those grains which lie along the path of the ionizing 
particle, hence no account will be taken of grains produced by 6-rays. 

Let p be the probability of development of any individual crystal, and 
let p be the average probability of development taken over all crystals 
traversed by the particle. If the crystal is not traversed by the particle 
we assume p=0. We assume further that for a given value of the average 
probability p, the probability of development p of any one crystal is 
independent of whether or not its immediate neighbours are developed. 
We must emphasize that for the moment we have made no assumptions 
as to how p might vary from one crystal to another, with degree of 
development or with the specific ionization of the particle producing the 
track. Thus, the probability per unit length of obtaining » undeveloped 
crystals bounded by two developed ones is 

m(n)=pe— 2" 


(6) 


where 5 is the mean spacing between crystals traversed by the particle. 

The second important factor determining the track characteristic is 
the distribution of the crystals along the path of the particle. The 
simplest model which can be advanced is that in which the distance 
between the centres of consecutive crystals is assumed to be a constant, §. 
It will be shown that this model is inconsistent with the experimental 
evidence. 

In a model which we advance later (§ 10) account is taken of a distri- 
bution in spacing s of neighbouring crystals traversed by the particle. 


§ 9. Constant Spactne MopeEn 

Provided that the centre of a developed grain coincides with the 
centre of the crystal from which it was formed, the distance between 
the centre of two developed grains, containing » undeveloped crystals 
between them, is (n+1)5=L, let us say. 

Thus, the differential distribution of gaps per unit length is a discon- 
tinuous function of the form 

p(1—p) 28 
§ 

which exists only for integral values of L/S. This distribution has an 
exponential ‘ envelope ’ of coefficient log, (l1—p)/é. 
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The gap length which has actually been measured is between the edges 
of developed grains. Fluctuations in developed grain size (and any 
displacement of the centre of a grain from the centre of the crystal from 
which it was formed), result in a smearing out of the above distribution, 
to a considerable extent. Such an effect is easily capable of accounting 
for the fact that the observed gap distribution is not noticeably irregular. 
Thus for the larger gap lengths, the integral distribution changes from a 
step-function to a function with an exponential envelope with the same 
coefficient as above, with only a slight ripple. For the values of gap 
length displayed in fig. 1, we can for all practical purposes identify the 
experimental slope with the coefficient of the exponent of the envelope ; 
thus we have 

pee ie p) (7) 

We have seen, from fig. 3 that g appears to tend to a limiting value, 
5000/mm, at high ionization. Taking 5=0-40u we have, according to 
- the above formula, p tending to 0-86. This value of § is taken to agree 
with the mean diameter of silver bromide crystals measured with the 
electron microscope. 

It is of course possible to employ a model in which 14% of the crystals 
are completely incapable of development. Calculations which we have 
made suggest that such a model would be capable of explaining all our 
experimental results. However, we prefer to explain the saturation of g 
entirely in terms of fluctuations in the individual values of the spacing, s. 


$10. VARIABLE Spactnc MopEL 


O’Ceallaigh (1954) was the first to make allowance for the fluctuations 
in the spacing between individual crystals. He concluded, under reason- 
able assumptions, that such fluctuations would not appreciably affect 
the exponential nature of the gap length distribution. He gave no 
-details of any effects which these fluctuations might have on the actual 
value of the coefficient, g. As we shall show, these are appreciable. 

As indicated above, we assume that in both the stacks employed, 
practically all crystals traversed by particles of Z>12 will have been 
developed. It follows therefore that the gaps observed in these tracks 
must be those between consecutive crystals. The distribution of such 
gaps is observed to be exponential, with a coefficient which we shall 
call 1/kj=5000/mm. We therefore assume that the distribution of 
length s between centres of successive crystals is of the form 


(9) ds=Uhy exp| — Jas ate. (3) 


for all values of s>Z, the mean path of the ionizing particle through a 
crystal. Thus the mean path length between centres of successive 


erystals will be S=k)+ 2. 
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It is known that at normal humidities the emulsion consists of nearly 
equal volumes of silver halide and gelatin. Thus the mean path of the 
particle through a crystal is very nearly equal to the mean value of the 
path & in the gelatin between two crystals, i.e. <= k=k j=0-20 ». The 
value of % so found is consistent with that expected for crystals of diameter 
0-30 4, as observed in electron microscope measurements. These values 
of ky and Z are to be compared with those of 0-1 » and 0-3 p respectively, 
taken by O’Ceallaigh. O’Ceallaigh’s model was based on the simplifying 
assumption that all traversals of crystals were along a diameter. 

For large values of s, formula (8) will result in an exponential distri- 
bution of the gelatin gaps between the developed crystals, as is observed. 
The true distribution however will not cut off sharply at s=#, since 
individual values of both x, the path of the particle through an individual 
crystal, and k, the path of the particle in the gelatin between adjacent 
crystals, will fluctuate. From consideration of the geometry involved— 
assuming no aggregation of crystals—the fraction of crystals along the 
path of a particle for which s<% would be unlikely to exceed 5%. Thus 
for the purpose of the detailed analysis which follows, it will be assumed 
that the distribution in s is rigorously given by eqn. (8). 

The distribution F',(w) du of the distance uw between the centres of 
two crystals which after development bound a gap containing n 
undeveloped crystals, is then 


F(u) du=(w"/n !) exp (—w) du for w>0 
==0 for w<0 Pee ek ED. 


where w= ulky— (n-+-1), and is the path length of the particle through 
the gelatin in the gap in units of ky. The probability per unit length 
of obtaining n such undeveloped crystals bounded by two developed ones 
is 7(n)=p?(1—p)"/s. Thus the contribution to the total gap distribution 
for each value of n is (n)F',(u) du. The form of this expression, for 
different values of n and for p=1—1/e, is shown in fig. 5 (curves 1, 2, 3 
and 4), 

The total differential gap length distribution is obtained by summing 
over all values of n 


su, pdu={ZECP” pwhau, eee Hy 


The integral distribution 4(u, p) of gaps exceeding wu can be obtained 
in a similar manner, though the expressions are more cumbersome. 

Both ¢(u) and ¢’(w) may be shown to become exponentials as w+ co. 
Thus the integral distribution becomes 


—, p* exp [—g(u—hy)] 
Duco (ts p)= rah (Oa See for rN ic (11) 


where the coefficient g is given by 


(l—gky) exp (—gky)=1—p. . . 2 2 (12) 
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ie generally, if we longer assume #=k,, we have instead of (11) and 
1 

—,__ p* exp [—9(u—2)] 

u->colU, Sas = = Be te eee }) 
eae as gk( 1+ &/ky—ga) Ce 
and (Gio exP ge) l== pes tee s- (120) 
Expressions (12) or (12a) are to be compared with (7) where the spacing 
is assumed to be constant. 


xe) 
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Theoretical gap-length distributions 


Probability of development = 1-4 


oY 


Probability —> 


fe) 
n 


-O3 


O02 
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Distance u between centres of developed crystals 
in units of Ko (=0-2p) 


Theoretical gap-length distribution (see § 10). Curves 1, 2,3 and 4, show the 
differential gap-length distribution corresponding to gaps containing 
0, 1, 2 and 3 undeveloped crystals. Curve (a) denotes the total differ- 
ential distribution, (b) the integral distribution. For u>4k,, both are 


sensibly exponential, 
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The full-line curves (a) and (b) in fig. 5 represent the asymptotic 
expressions for 4’(w, p) and ¢(u, p) normalized to unity at w=ky—eqn. (11). 
For values of u<4ky, the points represent the result of the numerical 
evaluation of eqn. (10). For w>4k,, corresponding to nearly all visible 
gaps, the analytic expression (11) and the results of the numerical calcula- 
tion are practically indistinguishable. They depart significantly for 
smaller values of uw, however. 

Let us now neglect all fluctuations of the size of the developed grains, 
and of the displacement of their centres from those of the original crystals, 
as in §8. Hence, all separations between developed crystals which 
exceed a certain value w1, will give rise to a visible gap. Then the blob 
density will be ¢(uw!, p) equals B,, say. Using the asymptotic expression 
(11), for this quantity, we obtain 

B= 2) | pe Oe Petes 13 

=gexp (gu) [ 1+ AEs Toes... (08) 
The series in brackets rapidly converges, and reaches a maximum value 
of e/2 when gk)=1, as can be readily seen by substitution in (11) and (12). 
For small values of g it is clear that the expression (13) tends to the 
empirical expression (2) for the blob density, and furthermore that « is 
to be identified with w1, which must be very nearly equal to the diameter 
of the developed grains as measured under normal optical conditions of 
the microscope. 

A further effect which should now be taken into account arises from 
fluctuations in developed grain size. Such variations are clearly visible 
in any track. The distribution in the observed diameters of single grains 
appears to be approximately Gaussian. The standard deviation p in 
stack 1 was found to be close to 0-14 u. The technical details and analysis 
of such measurements have been described by Della Corte, Ramat and 
Ronchi (1953 a). 

It can be shown that the effect of the distribution in developed grain 
size is to increase the predicted values of both B and H by a factor of 
approximately (1-++g?p?/4). We have no means of measuring the 
component along the track of displacements of the centre of the developed 
grain relative to the centre of the crystal from which it was formed ; 
but clearly, such displacements will result in an effective increase in the 
value of p to be used. The new expression for the predicted blob density, 
taking into account such fluctuations, is 


By=Bi(ict gp 4) gee eo eee Tas 


B, being given by (13). 


§ 11. COMPARISON OF THE VARIABLE SPACING MODEL WITH EXPERIMENT 


Figure 6 shows the experimental values of blob density B, plotted 
against the coefficient g of the exponent of the gap-length distribution. 
The observations cover singly-charged particles of velocity 8 between 1-0 
and 0:18, relativistic and slow «-particles, and a selection of heavier 
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relativistic cosmic ray nuclei. The points labelled ‘ fast protons’ are 
mean values taken from fig. 4. The points for slower protons and that 
for slow «-particles are also averages over a number of tracks. The 
points corresponding to multiply charged particles refer to individual 
tracks averaged over several emulsions. 


RELATION BETWEEN BLOB DENSITY AND 
COEFFICIENT OF GAP DISTRIBUTION 
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Relation between blob density B and the coefficient g, for the two stacks 
investigated. Curves (1), (2), and (3) show the predicted variation of 
these quantities, according to the variable spacing model (§ 10). In 
Stack 1, the points for heavy primary particles are taken from the 
results of fig. 1. 


The curves (1) and (2) in fig. 6 show the predicted variations of B, 
and B, with g according to (13) and (14) respectively, for a value of 
ui=0-67 ». Both curves are in good agreement with the experimental 
data from stack 1. ‘The best value of wu is thus slightly different from 
the value of 0-65 used for « in the empirical formula (2), Part I. 


SER. 7, VOL, 46, NO. 377.—JUNE 1955 20 
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In the same diagram we show some representative points taken from 
tracks in stack 2. Curve 3 shows the variation of B, with g, for a value 
of wi=0-76 w. 

The average values found for v1 in each stack are consistent with the 
mean values of the diameter of individual developed grains, as directly 
estimated. Such an estimate is, in fact, quite difficult to make and it 
appears that the most satisfactory way of determining the effective 
grain diameter, u1, is to measure the values of B and g for protons of 
~ 1 cmrange, where B is close to its maximum value for normal conditions 
of development. 

It should be pointed out here that the values of B and H observed along 
the tracks of particles of Z>8 are considerably lower than those predicted 
(see fig. 1). This we attribute to obscuration of parts of the track by 
the presence of 5-rays, which becomes serious for high values of Z. 


§ 12. INFLUENCE OF FLUCTUATIONS FOR SMALL VALUES OF g 


Taking into account the fluctuations described above, it can be shown 
that for low values of p the coefficient g can be approximated by the 


expression 
x y . Kee 


where y=log, (1—p)/8, and is of course the value which g would have had 
in the absence of fluctuations. This expression has a more general 
validity in that it is true whatever the form of the distribution of path 
length of the ionizing particle through either the silver bromide crystals 
or the intervening gelatin. /,)? should now be considered as the mean 
square deviation in individual values of the crystal spacing about the 
mean value 8. 


§ 13. VARIATION OF THE PROBABILITY OF DEVELOPMENT p WITH 
IONIZATION 


So far we have not considered how p might vary with ionization or 
degree of development, except in so far as we have assumed that it tends 
to unity as the specific ionization becomes very large (see § 10). 

In considering the variation of gy with the velocity and charge of the 
ionizing particle, it will have been noticed (i) that the value of g* for a 
relativistic «-particle is very nearly equal to four in both stacks investi- 
gated (fig. 3); (ii) that for proton ranges exceeding 1 cm the values of 
g* for the two stacks are almost indistinguishable, despite the large 
difference in development ; (iii) that for singly-charged particles g appears 
to vary approximately as dH/dR, the total rate of energy loss of the 
ionizing particle. 

It is natural to attempt to explain the variations of g in terms of the 
energy lost by the ionizing particle in individual silver bromide crystals, 
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i.e. the local ionization as distinct from the total rate of energy loss of 
the particle. 

On traversing the crystal, the charged particle produces a certain 
number of knock-on electrons, and excitations. If such knock-on electrons 
are below a certain critical energy, we shall call them ‘acts’. These 
acts will be distributed very nearly at random, since the number of 
atoms traversed per crystal is large. From any one act, one obtains a 
distribution in the number of free electrons liberated in the same crystal. 
d-rays of more than a certain critical energy will contribute very little 
to the ionization in the crystal from which they originate. We shall 
therefore neglect their contribution. This critical energy for 5-rays will 
be about 2 kev for emulsions with crystal diameters of about 0:3 4. This 
is much less than the maximum 5$-ray energy for the slowest particles 
investigated (G=0-18). Hence we feel it is reasonably close to the truth 
to assume that the distribution of the number of electrons resulting from an 
act, 1s independent of the charge and velocity of the primary particle for all 
tracks considered, and hence to assume that the number of acts per unit length 
is proportional to the local ionization produced. 

Let the average probability of development of the crystal as the result 
of a single act be g. q depends upon the degree of development and 
emulsion sensitivity. We assume however that it is independent of the 
size of the crystal and the position of the act init. This seems reasonable 
in view of the fact that the spatial distribution of the secondary ionization 
in the crystal about the path of the particle will be sharp. In any case, 
variations in q from this cause or from variations in intrinsic sensitivity 
of the crystal can be allowed for in terms of an extra variation of the 
path length x through the crystal. 

Let A be the average number of acts per unit length in the silver 
bromide. For a given value of x the distribution in the number of acts 
per crystal will be very nearly Poissonian, and therefore the fraction of 
those crystals with path length x which are developed is 


p=[1— exp (—Azq)]. 


We can obtain the average probability of development p for all crystals 
traversed by integrating this quantity over all values of x—or more 
accurately all values of the product xq. . 

We shall now calculate p for three possible distributions in z, as follows : 


(a) x= and is constant for all crystals traversed. 


(b) The crystals are spheres of equal diameter 32/2, distributed at 
random about the path of the particle. 


(c) The emulsion consists of spherical crystals of different diameter, 
all values of the diameter up to 2@ being equally probable (per unit 
volume of emulsion). 


2T2 
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All three models have the same value, #, for the average path length 
of the particle through a crystal. Let us write A®q=y and let f(x) dx 
be the distribution of the possible values of x. Then : 


(a) 1—p=exp (~Y), 
Sa eee 3 
(6) f(z) de= sad and 1—p= 5| 1 G y+ 1) exp (—3y/2) | 


ys 
= Cxp | —y—v8/16+ Evi tec | ; 


Ss 1 1 
and Dimes Ms - rey exp (—2y) | 
2 
—=exp | - CG a --terms in y* and higher) | : 


Model (a) has been proposed by a number of authors, but would appear 
to be an over-simplification. Model (b) first proposed by Demers (1947) 
appears at first sight to be very plausible ; for it is known that the 
silver bromide crystals are indeed nearly spherical, and that the dispersion 
of their diameters is not large. However this model must certainly 
represent an underestimate of the fluctuations involved which are those 
of the product ag. We therefore feel that model (c) is preferable, bearing 
in mind that there may be a distribution in the intrinsic sensitivity q 
of the crystals. 

It is now possible, at least in principle, to decide between the above 
three models by comparing their predictions of how g should vary with 
the development factor q, for a given value of ionization and therefore i. 

In fig. 2, the full line (1) is a smooth curve drawn through the experi- 
mental points of stack 2, for R<1lcm. From any given value of gs, 
the coefficient g in stack 2, we can calculate the corresponding quantity 
g, for stack 1, for each of the three models. The results are shown in 
the dashed curves (a), (b), and (c) of fig. 2. It is apparent that model (c) 
predicts differences between the values of g* in the two stacks in good 
agreement with experiment, though neither model (a) nor (b) could be 
excluded with confidence. The reason for this is of course that in the 
region of high ionization, where the three models diverge considerably 
from one another, the coefficient g becomes insensitive to variations in 
log, (1—p), as a result of the saturation due to the gelatin gap distribution. 
If however, one knew with confidence the variation of A with velocity 
one could decide among the three models quite easily. 


§ 14, VARIATIONS OF A witH RANGE oF PRoTons 
Using model (c), we have determined y=Aq# as a function of range R 
from the experimental points for stack 1, fig. 2. Making use of an extra- 
polated range-energy relation which takes into account the long-range 
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point of Heinz (1954), we have found that y and therefore A, the number 
of acts per unit length in silver bromide, varies very closely as 1/f?. 
The curve (2) in fig. 2 shows the relation to be expected between g* and 
_ & assuming A to be proportional to 1/8? in the range interval 1-60 cm. 
Models (a) and (b) give essentially the same result. The curve has been 
normalized at the energy used by Heinz (R=24-2 cm). It is interesting 
to note that the local ionization in the silver bromide crystals can be 
computed under reasonable assumptions, and also turns out to be very 
closely proportional to 1/6? in this range. 

From this experimental fact that A is proportional to the local ionization 
in the crystal, it appears that the processes of recombination of ions 
resulting from separate acts have little effect. 

For values of B>0-7, \ is certainly no longer proportional to 1/62. 
For if it were valid up to B=1, it would predict a value of g*=2-12 for 
R= 24-2 cm, as compared with the observed value of 1-59; and further, 
there would be no ‘ trough’ at B ~ 0-95. 

To investigate the variation of \ at low velocities, measurements must 
be made in less sensitive enulsions, such as C2. Observations that have 
been published (Peters 1952, Dainton, Fowler and Kent 1952) suggest 
that indeed still varies as 1/6”. 


§ 15. VARIATION OF A WITH CHARGE Z 


From the foregoing discussion, one might expect that the quantity A, 
the number of acts per unit length in silver bromide, should be propor- 
tional to Z?.. This is indeed approximately true. Careful measurements 
on the tracks of two extremely energetic «-particles gave a value of 
g*=4:10+0-05. Hence the value of A must lie between 4-3 (model (a)) 
and 4-55 (model (c)) times the value for a relativistic electron. This 
appears significantly greater than 4-0. Similarly, the value of A for an 
extremely energetic lithium nucleus was ~ 10 times the minimum value. 
Furthermore, in numerous measurements on «-particles of whatever 
energy, no case has been found in which g*<3-9. be 

From these provisional results, it would appear that there is little or 
no ‘trough’ in the curve of ionization versus energy for a-particles and 
that A/Ay is significantly greater than 4-0. Further observations are being 
made to investigate this point. 


§ 16. COMPARISON OF THE PRESENT RESULTS WITH THOSE OF OTHER 
WORKERS 


The distribution of gaps along the tracks of ionizing particles in G5 
emulsions has previously been studied by numerous workers. Among 
them are Menon and O’Ceallaigh (1954), O’Ceallaigh (1954), Hull, Happ 
and Morrish (1952), Herz and Davis (1954), Blatt (1954) and Della Corte 
et al. (1953 a, 1953 b, 1954). 

In the papers of Della Corte et al., and Herz and Damas; the models 
assumed constant crystal spacing as described above (§ 9). Such models, 
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as we have seen, fail simply because even when the ionization is very 
large, there is still a residual gap-length distribution, the greatest gap- 
lengths observed being farin excess of the mean spacing between consecutive 
crystals. As mentioned earlier one could however explain such gaps on 
the assumption that 14° of the crystals are completely incapable of 
development. This state of affairs we consider most improbable. 

In the papers of Della Corte et al. the probability of development 
appears to have been mistakenly identified with the reciprocal of the 
mean gap-length (i.e. the coefficient of the exponent of the distribution). 
This led them to assume a dependence of the mean enlargement of a 
crystal after development upon the probability of development. Our 
experimental evidence shows that this is not the case, at least up to 
values of ionization equal to those of relativistic oxygen nuclei. In 
leaving out of consideration the distribution in gelatin gaps between 
crystals, they have also had to appeal to recombination of the ions 
liberated in the crystals, in order to account for the saturation of the 
mean gap length at high values of ionization. 

As mentioned above, the effect of fluctuations in the gelatin gap 
between crystals had previously been considered by O’Ceallaigh. He, 
however, took a value of 1/kp, the reciprocal of the coefficient of the 
exponential gap distribution between successive crystals, of 0-1 uw, in 
place of the value 0-20 we have derived from the observations in heavily- 
ionizing tracks. There also appears to be an error in O’Ceallaigh’s 
calculated gap length distribution. His distribution has a maximum at 
wu ~ 5ky-10k,y instead of increasing continuously down to the minimum 
allowed value of w. 

Blatt (1954) has carried out an extensive theoretical treatment of the 
gap-length distribution according to a number of possible models. 
Comparing his results with a preliminary measurement of a single track, 
he had concluded that models based on variable crystal spacing were 
unacceptable. He has used for w!, the distance between centres of 
grains corresponding to the minimum detectable gap, the developed 
grain diameter as directly measured. The reason for the alleged discre- 
pancy appears to us as likely to have arisen from the difficulty of eliminating 
systematic errors from this direct measurement. Thus, he has taken a 
value of u!=0-65 4, whereas a value of 0-80 ~ would give agreement with 
the variable spacing models. 
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aE ee BiaN LE) SIGX: 
Calculation of Statistical Errors Associated with Blob and Gap Counts 


Let us consider the distribution in length, t, of the interval between 
the centre of the grain terminating one gap, and that of the grain 
beginning the next gap. The density of gaps is H, the minimum length 
measured between centres of the bounding grains is /+-u1. The mean 
length of all such gaps is /-+-u1+1/g=l, say. Let G be the number of 
developed crystals per unit length of track. G—=p/s, and can be found 
from eqn. (12) in terms of g. The interval, ¢, can be composed of any 
number of grains, n. The distribution of the frequency of occurrence of 
different values of n is a geometric progression : 


b(nj—=c(1—e)* 


where € and « are constants given by the conditions 


(ATE eae 
1 


and G=2Lnb(n)=€/e?. 
1 
Hence ea-ij/G and (€= HG, 


Let . be the mean value, and 7 the standard deviation about this mean 
value, of the distance between those adjacent grains, whose separation 
is less than J-+-w1 (i.e. whose separation is not large enough to count as 
a gap). If there are n grains in the interval ¢, the mean value of ¢ is 
(n—1)p, and the standard deviation about this mean value is n1/(n—1). 

We now require the total square deviation about the mean value of 
the quantity t, for all values of x. Let us first sum the squares of the 
individual standard deviations about t=0, over all values of ». The 
contribution from those intervals containing ” grains is 


p(n) I (n—1)n?+ (n—1)"n*]. 
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Hence the total square deviation about t=0 is 
V2=DE(1—e)"™A(n— LP (m— 17 SH) 
1 


2(1—e 2(1—e)(2—e 
nr(1 Jotul ess) 


€ e 


The mean length of the interval ¢ is 


i= Zyn)in— yu) Sp(n)— "C=? 


ih € 
Thus the total square deviation about the mean, f, is 


ee ietahs vectes3! onthe) # 
Riga Vem € + €* ey € (l—e)” 


So far, we have only considered fluctuations in ¢, the distance between 
the end of one gap and the beginning of the next. The total square 
deviation about the mean, of the distance from the beginning of one gap 
to the beginning of the next, is clearly (Y?—#)+1/g?._ This follows because 
the distribution of gap length is exponential, with coefficient g, and the 
standard deviation about the mean gap-length / is 1/g. The mean 
separation between beginnings of successive gaps being 1/H, we have for 
the standard deviation in H, 


dH 1 . 
is hy il Oe Oe an 
H —/ (Nw) v/ ( Ss /9”) 


ee #1, ae) 
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where N;,=number of gaps counted. The quantity 7 depends on the 
model of crystal separation employed, and is zero for the constant- 
spacing model. The term involving 7 for any model is always small in 
comparison with the sum of the other two terms under the square-root. 
e is equal to (//G) from above, and so depends slightly on the model used. 
t is the mean length of an interval and is therefore 1/H—l. 

The above formula, with the term involving 7? neglected, is given in 
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ABSTRACT 


! The following radiation characteristics of the neutron excess tantalum 
isotopes Ta and 1*°Ta and of the isomers 18°™W and 185mW have been 
determined : 


183’, 185 183m \V 185m \W 
Half-life 5-0-+0-ldays 49:54+1:5min 5-5sec 1-62-+0-05 min 
f-energy (Mev) ~0-15, 0-6 ~0-15, 1-72 — = 
y-energy (kev) 110,160,210, 125,175,235 105, 155 130, 165 
240, 320 


185T'a does not decay to any appreciable extent via the isomer 185™W and 
the 175 kev gamma-ray from 1*°Ta is not the same as the 165 key gamma- 
ray from 1™W. The 125 kev gamma-ray from 18>Ta and the 130 kev 
gamma-ray from °™W are probably the same and due to transition 
to the ground state of 18°W. Unless there is an undetected rate-deter- 
mining transition the 165 kev gamma-ray from 1%°™W is then due to the 
isomeric transition and from half-life and K conversion considerations 
it may be classified as E3. Similarly the 155 kev gamma-ray in 13™W is 
probably from the rate-determining transition.; The classification of the 
183mW isomeric transition type on the basis of K conversion measurements 
is not unambiguous. 


§ 1. INTRODUCTION 


THE neutron excess radioactive isotopes Ta and 1°Ta have been 
described previously. Butement (1950 and 1951) obtained half-lives of 
116 hours and 48 minutes for tantalum isotopes produced by (y, p) reactions 
on tungsten, and assigned them to 1°°T'a and 1*°Ta respectively. Duffield 
et al. (1950) irradiated tungsten enriched in the 1**W isotope with x-rays 
of energy up to 23 Mev and detected a 46-minute half-life tantalum activity, 
confirming the assignment to '°Ta. This had a_beta-disintegration 
energy of 1-7 Mev and a conversion electron energy of 75 kev. They also 
found a 1-85-minute half-life in the tungsten which they assigned to an 


* Communicated by the Author. 
+ But see note added in proof (p. 619). 
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isomeric state of 185°W, a 75 kev conversion electron again being associated 
with the decay. Moses and Martin (1951) also found the 48-minute 
tantalum activity after high energy x-ray irradiation of tungsten and 
detected a 150 kev beta-disintegration or conversion electron in addition 
to a 1-6 Mey beta-disintegration. The 150 kev and 1-6 Mey disintegrations 
were found not to be in coincidence. They also found an activity with a 
6-day half-life which had a beta-disintegration energy of 0-6 Mey. 

The work described in this paper was undertaken to examine the 
radiations from 18°Ta and 18°Ta in greater detail. Work on the isomers 
1s3m\W and 185mW is also described. During the course of this investigation 
the assignment of the 5-day activity to 1*°Ta was confirmed by Mihelich 
(1953) and DuMond et al. (1953) who detected it after irradiation of tanta- 
lum in a high flux of thermal neutrons, the *°Ta being produced by double 
neutron capture. The latter authors found the half-life to be 5-2+-0-1 
days and described the gamma-ray spectrum in great detail. 


§ 2. EXPERIMENTAL TECHNIQUES 


The tantalum isotopes were prepared by irradiation of tungstic acid with 
fast neutrons, generated in the Harwell cyclotron, of energy up to 50 Mev 
which produced (n, p) and (n, pn) reactions in the tungsten. The tungstic 
acid was then dissolved in sodium hydroxide solution and the tantalum 
coprecipitated with ferric hydroxide by the addition of a solution con- 
taining a few milligrams of ferric chloride. After washing with sodium 
hydroxide solution the precipitate was dissolved in hydrochloric acid and 
reprecipitated from a solution of inactive sodium tungstate acting as a 
hold-back. 

The radiation characteristics were examined by aluminium absorption, 
beta- and gamma-ray scintillation spectroscopy and coincidence counting 
techniques. The anthracene crystal in the beta-particle scintillation 
spectrometer was a cylinder 0-5 in. long and 1-5in. in diameter. The 
spectrometer had a resolution of 20°, for the 624 kev conversion electron 
from a 7Cs—187™Ba source. The main source of error in the measurement 
of beta-disintegration energies was in the initial calibration with a con- 
version electron of known energy, the final Fermi plot derived from the 
beta-particle spectrum usually being very close to a straight line. 

Conversion electrons were examined by absorption of those electrons 
which were in coincidence with beta-decay. The sample was mounted 
over a small hole in a 1 mm thick brass plate placed between two end- 
window Geiger counters, the brass absorbing most of the particles which 
would otherwise have been scattered from one counter into the other, 
thereby producing a spurious coincidence. 

The gamma-ray scintillation spectrometer had a resolution of 10% for 
the photoelectric absorption peak from the 661 key gamma-ray from 
87Cs-l8mBa, ‘The sodium iodide crystal used was a cylinder 1-0 in. long 
and 1-5 in. in diameter. 
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§ 3. RESULTS 
(i) 188Ta, 

The half-life of this isotope was found to be 5-0-0-1 days from the 
beta- and gamma-ray decay of ten samples followed for several months. 
After subtraction of the long-lived activity, which accounted initially for 
about 5%, of the counting rate, the decay of the 1Ta was accurately 
measurable for about five half-lives. 

Aluminium absorption experiments indicated a rather indeterminate 
low energy component of 0-1—0-2 Mey in addition to the maximum beta- 
particle energy of 0-6 mev. The low energy component accounted for 
30-40%, of the total beta-particle counts. 

Scintillation spectrometer studies of the gamma-ray spectrum of a 
tantalum source after shorter-lived activities had decayed away indicated 
gamma-ray energies of 60, 110, 160, 210, 240 and 320 key with approximate 
relative intensities 25 : 5: 5:3: 10:7, each of which decayed with a 5-day 
half-life. The 60 kev component is very probably the K x-ray. This 
spectrum agrees with the results of DuMond e¢ al. (1953) allowing for the 
relatively much poorer resolution of the scintillation spectrometer. The 
group of gamma-rays of energy close to 145 kev reported by DuMond et al. 
was not resolved but could have been present up to an intensity of 2-5 
relative to the 240 kev gamma-ray intensity of 10. 


(a4) 1°Ta 

The half-life of 1°Ta was measured by following the decay of those 
beta-particles which passed through 500 mg/cm? of aluminium. This was 
sufficient aluminium to absorb all the beta-particles from longer-lived 
isotopes present and the long-lived background of gamma-radiation 
amounted to about 10° of the initial counts. The average half-life 
resulting from six determinations was 49-5 minutes, all the results lying 
within 1-5 minutes of this value. 

The beta-particle spectrum was examined by aluminium absorption and 
scintillation spectroscopy. The former indicated, in agreement with 
previous published work, a maximum beta-particle energy of 1-7 Mev and a 
low energy component of about 0-15 mev which amounted to about 30% | 
of the total count. The absorption curve of those particles which were in 
coincidence with beta-particles of energy >0-2 Mev showed mainly the 
0-15 Mev component. This is therefore a conversion electron from a 
transition in 185W which is in coincidence with the beta-disintegration of 
185Ta, The absorption of those beta-particles which were in coincidence 
with the total gamma-radiation was indistinguishable from the normal 
beta-particle absorption, which is a further indication of a transition in 
185W which is in coincidence with the beta-disintegration of *5°Ta. 

The beta-particle spectrum measured with the scintillation spectrometer 
is shown in fig. 1 together with the corresponding Fermi plot. W is the 
total energy of the electrons in mc? units. N(W) is the number of 
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electrons with energy between W—4x and W+-42, where « is the channel- 
width of the pulse analyser in energy units. is the Fermi beta-disintegra- 
tion function calculated by Bleuler and Ziinti (1946). The plot of 
[N(W)/@W?2]/2 against kinetic energy is very close to a straight line, the 
deviation at the maximum energy being probably due to statistical line 
broadening (Palmer and Laslett 1951). The average of five closely con- 
cordant results for the maximum beta-particle energy was 1-72 Mev. 


Fig. 1 


1-0 1.5 
Energy Mev 


Beta-particle energy of 49-5min 1°Ta (a) N(W)x 10-8, (b) [N(W)/eW2})?, 


The gamma-ray spectrum is shown in fig. 2 and indicates gamma-rays 
of 125, 175 and 235 kev and also a 60 kev component (presumably the 
K x-ray) in the approximate relative intensities, referred to the 175 kev 
gamma-ray, of 4:10: 2-5:14. The sample for this work was prepared 
by irradiation of tungstic acid enriched in 18°W and the spectrum of longer- 
lived activities was weak and easily corrected for. 
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(viz) 185m \V 

The *°™W isomer was prepared by irradiation of tungsten metal with 
neutrons of energy up to 50 Mev in the cyclotron. A half-life of 1-62-+0-05 
minutes was obtained from the gamma-ray decay. 

The mass assignment of this isomer to 18°™W was made by Duffield e¢ 
al. (1950) on the basis of its production from tungsten enriched in 18°W by 
irradiation with x-rays of up to 23 Mev. The possibility of a theoretically 
unlikely assignment to 1*6mW has been tested by measuring the change 


Fig. 2 


I25KeV 
Y 


Intensity in counts/channel 


100 200 
Energy kev 


Gamma-ray energy spectrum of 49-5 min 1°Ta, 


with neutron irradiation energy of the relative yields of 1-6-minute 
activity from two tungstic acid samples, one enriched in 18*W and the 
other in 18°W. The ratio of the amounts of 184W in these two samples was 
58:1 and of 486°W was 1:16. If the 1-6-minute isomer were eee or 
1s6mW the yields from the two samples after irradiating and counting 
under identical conditions would be in one of the ratios 58: 1 or 1: 16, and 
there would be no change in the ratio with change in the energy of the 
irradiating neutrons. If the isomer were *°™W the ratio of the yields 
would be between 58 : 1 and 1 : 16, depending on the relative cross sections 
of the reactions 184W(n, y)!8°™W and 1°°W(n, 2n)*°™W under the condi- 
tions of the irradiation. At low neutron irradiation energies the (n, y) 
reaction will become increasingly more probable relative to the (n, 2n) 
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reaction as the irradiation energy is decreased, and the ratio of the yields 
of the isomer from the two samples will move away from 1: 16 towards 
Sith 

The ratios obtained from neutron energies > 20, +15 and + 10 Mev were 
1:11, 1:47 and 1: <0-3 respectively, the yields having been measured 
by counting only those gamma-rays of energy between 40 and 250 kev. 
(The gamma-ray spectrum of the 1-6-minute isomer is described below.) 
The irradiations were all monitored using the 12-8-hour activity induced 
in copper foils. Irradiation of CaCO, containing a mass of oxygen equal 
to that in the tungstic acid samples showed that no correction for the 
counts from the annihilation radiation from 1°O was necessary. The 
isomer is therefore almost certainly °™W. 

Attempts were made to separate off any of the *°™W isomer formed after 
beta-decay of 48°T'a using a Deacidite F—F anion exchange resin. This was 
used in two ways. In one the cyclotron-produced tantalum activities 
were dissolved in a little 3M HCl: 0-5M HF solution and mounted on a 
column of resin 5 em long and 1 em in diameter. This column was eluted 
at 20-minute intervals with 5 ml. of the HCl : HF solution and the solution 
obtained transferred to the counter in less than one minute. These 
samples were counted in fixed geometry with a scintillation spectrometer 
adjusted to count gamma-rays between 40 kev and 250 kev. In the other 
method the tantalum activities were dissolved in 10ml. of the 
3M HCl: 0-5M HF solution and shaken for 15 minutes with 5 g of the 
resin. The solution was then filtered from the resin and counting begun 
within 75 seconds of the separation. In both methods the total activity 
of the initial sample was measured in the same geometry before mounting 
on the resin and allowance for the decay made so that the 49-5-minute 
185'T'a activity was known at the time of each attempted 1*°™W separation. 

With tracer activities of tungsten and tantalum under the conditions 
described above, not less than 10°% of the tungsten activity and only up to 
0-4°% of the tantalum activity was eluted from the column, and not less 
than 15°, of the tungsten and only up to 0-6% of the tantalum was 
extracted from the resin by shaking for 15 minutes. Similar results 
were obtained with the elution of carrier free 18'W from the column. 

In the experiments with cyclotron produced tantalum activities no short- 
lived decay of the extracted samples was observed. The difference 
between the first count after separation and a count after all the 185™W 
would have decayed was of the same order as the probable errors involved 
in the counting. Knowing the total activity on the resin the activity of 
the samples obtained was that expected if only tantalum activities were 
present. The proportion of gamma-ray counts from 49-5-minute 1Ta 
which are due to decay from the meta-stable 18°™W level must therefore 
be very small. Assuming that 10°% of the available 45™W was eluted 
from the column and that 15% was extracted from the resin by shaking 
with acid, the average of fifteen values for the proportion of 185™W countsin 
Ta is 0-1% with a probable error of +0-5%. The detection efficiency, 
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by gamma-ray counting, of the 1*5™W decay may well be lower than that 
of 18°Ta decay. The proportion of 1°Ta disintegrations which lead to 
isomW may therefore be somewhat higher than the proportion of 1%>Ta 
gamma-ray counts which are due to !8°mW decay. 

Coincidence measurements described above (3 (ii) provide additional 
evidence that the '*°T'a beta-transition does not go directly to the 185™W 
isomer. The gamma-ray spectrum of those gamma-rays which were in 
coincidence with the 1-7 Mev beta-disintegration of 18°Ta was indistin- 
guishable from the normal spectrum. 

The gamma-ray spectrum of a tungsten sample measured two minutes 
after irradiation, using a thirty-channel pulse analyser, is shown in fig. 3. 


Ric. 3 


130 ia 165KeV 


y 


Intensity in counts/channel 


50 100 I50 


Energy kev 


Gamma-ray energy spectrum of 1-6 min 15my, 


There is a 60 kev K x-ray and 130 and 165 kev gamma-rays in the relative 
intensities, after small corrections for background and long-lived activities, 
of 2:1:1. The spectrum of gamma-rays in coincidence with the 165 kev 
gamma-ray was measured with two scintillation counters placed. with the 
two crystals close together. One counter was connected to a single- 
channel pulse analyser which only accepted pulses from 165 kev gamma- 
rays, and the other to a one hundred channel pulse analyser which only 
accepted gamma-ray pulses which were in coincidence with the 165 kev 
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gamma-rays detected by the other counter. The resulting spectrum 
consisted only of the x-ray and the 130 kev gamma-ray. The 165 and 
130 kev gamma-rays are therefore in coincidence. 

Low lying excited states of the tungsten isotopes of mass 182, 183, 184 
and 186 have been shown by Coulomb excitation to be 101, 103, 112 and 
124 kev respectively above the ground state (McClelland, Mark and 
Goodman 1954, Huus and Bjerregaard 1953). It therefore seems probable 
that the 125 kev gamma-ray from 1°Ta and the 130 kev gamma-ray - 
from 185mW are the same, and due to transition to the ground state of 
18s5W, It may then be concluded that the 165 kev gamma-ray from 
IssmW is the rate-determining step in the isomer. The 175 kev gamma- 
ray from 18°Ta must then be from another transition involving other 
excited levels in the 185W nucleus, since it is in coincidence with the 
185Ta, beta-disintegration. It is possible, however, that there may be 
an as yet undetected rate-determining transition in the isomer, in which 
case the 165 kev gamma-ray could account for part of the 175 kev 
gamma-ray observed in the °Ta-spectrum up to an extent equal to the 
intensity of the 125 kev gamma-ray (see figs. 2 and 3). 

Assuming, however, that the 165 kev transition is the rate-determining 
transition in the 185mW isomer, the half-life and energy define it as E3 or 
M3 (Goldhaber and Hill 1952). According to the calculations of Rose et 
al. (1951) the K conversion coefficients of 165 kev E3 and M3 transitions 
are 0-7 and 27 respectively. 

The K conversion coefficient, as indicated by the gamma-ray: X-ray 
intensity ratio, cannot be greater than 2 and the E3 assignment for the 
1-6-minute 18°™W isomeric transition is therefore strongly favoured, 
The 1-6-minute isomer is therefore probably a 165 key E3 transition 
followed by a 130 kev transition to the ground state. 


(2p) Asm 

The half-life, mass assignment and conversion electron energy of this 
isomer have been reported by der Mateosian et al. (1949). The half-life 
and energy identify the transition as E3 or M3 (Goldhaber and Sunyar 
1951), and unpublished work by der Mateosian and Goldhaber indicating 
that mainly L x-rays arise in the conversion supports an E3 transition 
type assignment (Goldhaber and Hill 1952). 

The isomer has been prepared in this laboratory by irradiation of tung- 
sten with fast neutrons in the cyclotron, and the gamma-ray spectrum 
measured using a scintillation spectrometer and a single-channel pulse 
analyser. Pulses from the amplifier were fed into the pulse analyser and 
also directly into a scaler. The samples were removed from the cyclotron 
and counted from ten to twenty, sixty to seventy and eighty to ninety 
seconds after the end of the irradiation. The ratio of the number of counts 
of a given energy from the isomer to the total number of counts from the 
isomer was thereby obtained. Each irradiation gave one point on the 
spectrum which was covered after about twenty such irradiations. 
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The spectrum so obtained indicated gamma-rays of 105 and 155 kev 
together with the K x-ray in the approximate relative intensities 
2-5:1:10. The decay ofeach of these components was followed separately 
and indicated a half-life of about 5-5 seconds. Campbell and Goodrich 
(quoted by Hollander, Perlman and Seaborg 1953) have reported gamma- 
rays of 0-12 and 0-17 mev for this isomer. 

Following similar arguments to those used for 18°™W, the 105 key 
gamma-ray is probably due to transition to the ground state. The 
155 kev gamma-ray is therefore associated with the isomeric transition 
unless there is an undetected rate-determining transition. The K 
conversion coefficients for E3 and M3 transitions of this energy are 
0-8 and 32 respectively (Rose e¢ al. 1951). The K conversion coefficient 
of the 155 kev transition, as indicated by the x-ray: gamma-ray 
intensity ratio, could be as high as 10.* On this basis therefore the 
choice between a transition type assignment of E3 or of M3 cannot be 
made with certainty, although by analogy with 1°™W an E3 assignment 
is favoured in agreement with the conclusions of der Mateosian and 
Goldhaber quoted above. 


ACKNOWLEDGMENTS 


Acknowledgments are gratefully made to Dr. F. D. 8. Butement for 
suggesting this work and for helpful discussions and advice, to the 
members of the Cyclotron Group for carrying out irradiations and to 
Mr. R. B. Thomas for much assistance with the experimental work. 


* Note added in proof. 


A detailed level scheme for 18°W recently given by Murray et al. (Murray, J. J., 
Boehm, F., Marmier, P., and DuMond, J. W. M., 1955, Phys. Rev., 97, 1007) 
shows that decay to the ground state from the 99-07 kev level reached by 
Coulomb excitation can occur by a 52-59 kev-46-48 kev cascade as well as 
directly. The 60 kev component from *°™W must, therefore, be partly due to 
these two low energy gamma-rays. Consideration of the 1*°W level scheme 
_ suggests that the 155 kev gamma-ray from 4°™W may well not be from the rate- 
determining transition whose energy and exact transition type are, therefore, 
still uncertain. The 1%5™W decay scheme is probably also more complex than 
suggested above. 
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SUMMARY 

The theory of the coercivity of non-spherical, single domain particles 
of ferromagnetic material in a non-magnetic matrix has been developed 
independently by Néel (1947) and Stoner and Wohlfarth (1948). Appro- 
priate conditions for the study of the growth of these domains are to be 
found in suitable dilute ferromagnetic alloys where there is an increase 
of solubility of the ferromagnetic constituent with temperature. Such 
characteristics are found in the copper rich alloys of copper—cobalt and 
copper—iron. The magnetic properties of appropriate alloys with increase 
of heat treatment are followed from the normal paramagnetism of the 
solid solution to the ordinary magnetization curve of the precipitated 
ferromagnetic. At a suitable stage a maximum in the coercivity is found 
which the writers believe to be associated with single domain precipitates. 
The characteristics of the hysteresis loops at this stage indicate the 
presence of high coercivities of more than 1000 oe. 


§ 1. INTRODUCTION 

In low fields changes in magnetization in ferromagnetic materials occur 
by two distinct processes; by domain wall movements in which the 
volumes of the domains change and by rotation of the domain magnetiza- 
tion vector. Various theories (e.g. Becker, Kersten, Néel) have been put 
forward to provide the mechanism whereby the free magnetic energy 
controls the position of the domain boundary but in their present form 
these theories cannot easily account for coercivities of more than some 
hundreds of oersteds. However, because of the large energy changes 
accompanying rotational processes, in the absence of wall movements 
high coercivities may be associated with irreversible rotational operations. 

Single domains were first predicted by Frenkel and Dorfman (1930) 
who suggested that in small particles the formation of a domain wall 
may be energetically unfavourable and that such a region would be 
magnetically uniform. Theoretical analyses by Néel (1947) and inde- 
pendently by Stoner and Wohblfarth (1948) have shown that coercivities 
of several thousand oersted can be attributed to shape anisotropy. The 
theoretical maximum coercivities associated with crystal, shape and 
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strain anisotropies with complete orientation of the anisotropy axes are 
shown in table 1 (Kittel 1949); for a body containing randomly 
orientated particles, the values of coercivity are reduced by a factor 
of 0-48. Single domain theories of coercivity are applicable to very 
finely divided powders and to heterogeneous alloys in which magnetic 
‘islands’ are dispersed in a non-magnetic matrix. 


Table 1 


Anisotropy Crystalline Shape 


Fe 500 10700 
Co 6000 8800 
Ni 135 3150 


From the results of experiments on finely divided powders, summarized 
by Weil (1951), it appears that there is good evidence for the crystalline 
and shape anisotropy mechanisms of hysteresis although the values of 
coercivity are less than would be expected from theoretical considerations. 
Little work has been published on suitable heterogeneous alloys, though 
Constant and his collaborators (1939, 1943) investigating the low field 
properties of impure copper and brass reported coercivities of about 
100 oersted. 

A preliminary investigation by Sucksmith (1952) into the variation in 
magnetic properties of stainless steel magnetic recording wire, 12% Cr, 
12%, Ni, rest Fe, with heat treatment and cold work, showed that it 
was desirable to study very dilute two-phase alloys in which interaction 
effects were negligible and the growth of ferromagnetic regions controllable 
by heat treatment. The required characteristics appeared to be present 
in copper rich alloys of copper—cobalt and copper—iron in which the 
solubility of the ferromagnetic component rises from a very low value 
at room temperature to 3-4% at 1000°c. Thus, the rapid quenching of 
alloys of compositions of up to 3°%, cobalt or iron produces a supersaturated 
single phase condition which can be transformed into the two-phase 
equilibrium state by appropriate heat treatment. The present paper 
contains an account of an investigation into the variation with heat 
treatment of the magnetic properties of copper—cobalt and copper-iron 
containing small percentages of ferromagnetic component. 

In the determination of the hysteresis loops of such materials with 
very low intensities of magnetization, induction methods are unduly 
difficult, chiefly due to inadequate sensitivity, and the utilization of 
‘traction ’ methods appears to offer the best prospect of success. Two 
different methods have been investigated: (i) the couple due to a 
transverse field acting upon the moment of a specimen and (ii) the force 
exerted on a magnetized specimen due to a field gradient. Though the 
main interest lies in the low field characteristics, a few high field results 
obtained on a paramagnetic balance are given, 
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§ 2. MEASUREMENT or Hysteresis Loops 
Preparation of the Alloys 


The alloys were prepared by melting the constituents in alumina 
crucibles in a high frequency induction furnace in vacuo or in an 
atmosphere of argon at a pressure of 5 cm of mercury. 


Heat Treatment 


The ingots were forged at 900°C into bars approximately 
lemxlemx7cm. The specimens for magnetic measurements were 
discs 5mm diameter and 0-5 mm thick. <A scratch on the surface of 
_each dise provided a reference mark. Heat treatment was carried out 
with the specimen sealed in a short length of pyrex tubing in a furnace 
maintained at a constant temperature to +1°c. After heat treatment 
the surfaces of the specimens were dull but in no case did appreciable 
oxidation take place and comparison of results from control experiments 
on specimens annealed in vacuo indicated that errors from this source 
were negligible. 


Determination of Hysteresis Loops 


Constant and collaborators (1939, 1943) used a method which consists 
of measuring the couple exerted on a specimen by a small field applied 
perpendicular to the magnetizing field. Preliminary experiments showed 
that in our experiments this method was unsatisfactory, and some 
discussion of this method is presented in the Appendix. On the other 
hand, the technique described below was adopted and found to be 
satisfactory. 

If a ferromagnetic specimen in a field H, is subjected to a field 
gradient dH_,,/dx, the force exerted on the specimen is 


dH 


dius 


With representative values for volume and intensity, in a field 
gradient of 100ecm~!, the force to be measured is of the order of 
10-8 dynes. Such small forces of this magnitude are most easily measured 
by a torsion balance. A schematic arrangement of the apparatus is 
shown in fig. 1. 

The dise specimens were held in a recessed hole at the end of a perspex 
arm 5cm long and counterbalanced by a glass weight. Considerable 
care was taken to keep the moment of inertia of the balance as small 
as possible and by using a simple damping vane immersed in a light oil 
to ensure that the system was near critical damping. 

The balance arm was suspended from the torsion head by a phosphor 
bronze wire 5cm long, the torsion constant of which was determined 
by a separate experiment from the period of oscillation of a metal annulus 
of calculable moment of inertia attached to the suspension wire. To 
cover the range of intensities four gauges of phosphor bronze wire from 
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0-025 mm to 0:06 mm diameter were required. The balance was pro- 
tected from draughts by enclosing the torsion system in a rectangular 
box which also limits the swing of the arm. The rotation of the balance 
and torsion head were determined by the use of small galvanometer 
mirrors attached to the balance arm and torsion head. 


Fig. 1 
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Arrangement of Apparatus. 


The magnetizing field was furnished by two coils rigidly fixed to a 
base plate with the inner faces lin. apart. The field at the mid-point 
of the axis determined by a search coil and ballistic galvanometer was 
101 oe/amp ; a maximum current of 8 amp could be passed for a few 
minutes. The, gradient coils were mounted in the Gaugain—Helmholtz 
position co-axially on the magnetizing coils and symmetrically about 
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the mid-plane of the magnetizing system. The coils connected in 
series opposition were usually excited by a current of 3 amp which by 
calculation produced a field gradient of 7-02 0ecm-! at the position 
occupied by the specimen. Experiment proved that the hysteresis loops 
‘were independent of the magnitude of the field gradient used in the 
determination. 


Hxperimental Procedure 


The specimen was magnetized to saturation in a small electromagnet 
capable of producing about 5000 oe and then transferred to the torsion 
balance. The orientation of the disc was adjusted so that the direction 
of this magnetization and the applied field in the hysteresis loop apparatus 
were collinear. The magnetization curve was obtained by a step by step 
reduction of the field from a large positive value to suitable negative 
values, the field gradient being kept constant. At each value of the 
magnetic field, the torsion head was adjusted to return the specimen 
to the predetermined zero. By the switching on of a suitably small 
field gradient, the specimen was deflected and again brought back to 
the zero by rotation of the torsion head. From a knowledge of the 
rotation of the head the force exerted on the specimen and hence the 
intensity of magnetization could be calculated. No correction was 
necessary for the demagnetizing effects because of the very low intensities. 


Experimental Results 


In each binary system measurements of the variation of magnetic 
properties of a series of alloys 0-7% to 2-0% for copper—cobalt and 
0-3% to 1:5% for copper-iron with time of heat treatment at 350°c, 
400°c and 470°c in the former case and 500°c, 600°c and 800° in the 
latter were carried out. 


300 100 200 
FIELD OF — > 


J-H Curves of 1% Co-Cu alloy. Heat treated at 350°C. 
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In the quenched state, the alloys exhibited paramagnetism almost 
independent of field and with a mass susceptibility of about 10 *sa0n 
heat treatment development of curvature in the J—H graph took place 
whilst the order of magnitude of the susceptibility indicated the presence 
of ferromagnetism. 


Fig. 3 


J-H Curves of 1% Fe-Cu alloy. Heat treated at 650°c. 


Fig. 4 
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Coercivity v. Time at 650°c for 1-:0% Fe-Cu alloy. 


The changes in magnetization curves with time of heat treatment 
for typical alloys 1% Cu-Co heat treated at 350°c and 1% Cu-—Fe 
annealed at 650°C are shown in figs. 2 and 3. In copper-cobalt alloys 
with continued heat treatment the coercivity increased, reached a 
maximum and then decreased. The iron alloys behaved similarly except 
that an initial decrease in coercivity was always observed (fig. 4). 
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Though the largest coercivities observed were 2500e and 400 oe 
for Cu-Co and Cu-Fe respectively, it was found that fields of several 
thousand oersteds were required to develop maximum coercivity and 
remanence. Figure 5 gives the results on a typical specimen of 
- copper—cobalt, a specimen with a coercivity of 74 0e, though the 
hysteresis loops show unusual features. 


Fig. 5 
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J-H Curves for 14°% Co—Cu alloy heat treated for 31} hrs. at 350°c. 
Magnetization fields shown on curves. 


Whilst in normal ferromagnetics the application of a maximum. reverse 
field not much greater than the coercive field is adequate to delineate a 
hysteresis loop and to change the sign of the remanent magnetization, 
such a field applied to these alloys only produced a slight diminution of 
the magnitude of the remanence, and fields up to 1500 oe were 
required to reduce the remanence to a Zero value, whilst to reverse 
the remanence in both magnitude and direction called for fields of 
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about 5000 oe. Figure 6 gives data on a 1% copper-cobalt alloy, the 
coercivity of which under similar conditions to those shown in fig. 5 
were of the order 150 oe. 

The approach to saturation in high fields was relatively slow, and 
we are indebted to Dr. J. Crangle of this department for high field data. 
Above about 3000 oe it was found that the increase of magnetization 
with field was approximately linear up to about 15 000 oe. 


Fig. 6 
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Permanent Magnetization of 1% Co—Cu alloy, heat treated at 350°c. 


§ 3. Discussion 


When compared with the behaviour usually associated with bulk 
ferromagnetic materials some of the characteristics of these alloys appear 
anomalous. One feature of all the magnetization curves is that the 
remanent intensity is a small fraction of the intensity in high fields. 
A second unusual characteristic is that, though in general a field about 
four or five times the coercive force should be adequate to saturate even 
a permanent magnet material, some of these alloys required fields in 
excess of twenty times of coercivity to exhibit maximum remanence and 
coercive force. In addition, in fields of several thousand oe, the magneti- 
. zation of the alloys increases approximately linearly. This field indepen- 
dent incremental susceptibility was also characteristic of similar alloys 
investigated by Bitter and others (1939, 1941). In both cases it was 
found that this susceptibility increased with time of anneal. 

It is possible to draw a generally satisfactory picture of the processes 
involved, though the details are more obscure. As the heat treatment 
commences upon the solid solution, the iron or cobalt atoms are so widely 
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dispersed that the cooperative phenomenon of ferromagnetism is im- 
possible. As the heat treatment progresses, atoms of the ferromagnetic 
element diffuse and aggregate to form regions of higher concentration. 
As these nuclei grow, the threshold of ferromagnetism is reached with 
_ the formation of single domains and the theories of Stoner and Wohlfarth 
and Néel may become applicable. During this stage we get the gradual 
transition of paramagnetism to ferromagnetism with marked increase of 
coercivity owing to the shape anisotropy of non-spherical single domains in 
the non-magnetic matrix of copper. It is to be expected that differences 
in behaviour will be observed according to the crystal structure relation- 
ships of the two phases. Work on other suitable binary alloys of all three 
ferromagnetics are in progress to try and elucidate further this aspect of 
the problem. The case of cobalt binary alloys is particularly complicated 
since this element has two crystalline structures, face-centred-cubic and 
close packed hexagonal, both of which are ferromagnetic and can exist at 
temperatures up to about 400°c. ; 

As heat treatment continues the volume of a ferromagnetic region 
increases until it is energetically more favourable for the region to divide 
into two or more domains. Changes of magnetization of such a multi- 
domain region proceeds by the movement of domain walls such that the 
component of magnetization in the direction of the field increases. This 
process gives rise to low coercivities and in small, approximately spherical 
particles the changes may be almost completely reversible. 

It is obvious that nuclei will develop at different rates owing to local 
variations in concentration so that the processes described above will be 
simultaneous and composite magnetization curves will result from the 
various stages of aggregation, ranging from paramagnetism through a 
wide range of ‘ hard’ and ‘ soft.’ hysteresis curves. From the composite 
magnetization curves it is not possible to separate the hard and soft 
components but the presence of a ‘hard’ magnetic component is 
indicated by the large fields required to produce maximum remanence 
and coercivity and the high reverse fields necessary to reduce to zero 
or change the direction of the remanent magnetization. The low ratio 
of remanent to saturation magnetization can be safely regarded as a 
consequence of the range of properties of the aggregates. Only a 
relatively small percentage of these make substantial contributions to 
the area under the hysteresis loop. 

Thus heat treatment of the supersaturated alloys causes the atoms 
of the ferromagnetic material to aggregate, eventually reaching the 
critical size of ferromagnetism and being sufficiently small, form single 
domains which will exhibit a large coercivity if the regions depart from 
spherical shape. With further heat treatment the particles continue to 
grow and become multidomain. On this basis the variation of coercivity 
with heat treatment is readily explained ; as the proportion of regions 
behaving as single domains increases the coercivity rises, reaches a 
maximum and falls as the fraction of multidomains increases. 
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APPENDIX 


From the definition of intensity of magnetization, the couple M 
exerted on a ferromagnetic specimen of volume v in a field H at an 
angle @ to the direction of magnetization is given by 

M=HvJ sin 0. 

If a specimen is placed with the direction of magnetization, J,, collinear 
with the magnetizing field, H,, the resulting couple is zero. The application 
of a small transverse field, H,, exerts a couple which when the specimen 
is returned to the zero is J,vH,. 


Fig. 7 
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J—H Curves by the two Methods. 1° Fe-Cu—Heat treated for 70 hours at 
650°C. 


A technique for determining hysteresis loops based on this system 
has the merit of simplicity and experiments were carried out to determine 
its validity. On obtaining such curves as shown in fig. 7, in which the 
apparent intensity decreased with increasing field, it was obvious that 
there was a fundamental flaw in the method. Moreover the measured 
intensities were not independent of the magnitude of the transverse field. 

The transverse field method tacitly assumes that the magnetization 
vector remains fixed relative to the specimen. However, Sucksmith, 
Potter and Broadway (1928) in a criticism of a similar technique for the 
determination of the parallel component of magnetization of a single 
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erystal of nickel pointed out that the effective field on the specimen is 
the resultant of the applied fields and the demagnetizing fields. For 
an inhomogeneous alloy in which the ferromagnetic regions are isolated, 
a complete interpretation of the results requires the consideration of 
the demagnetizing field of each individual region. In the dispersed 
systems of aggregates whose shape and distribution is random the above 
method would appear to be inapplicable and was therefore abandoned, 
and all measurements carried out by the method described earlier in 
which the translational force between field and moment are collinear. 
Figure 7 shows curves obtained by both methods on the same specimen. 
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SUMMARY 


A combination of lens spectrometer and cloud chamber was used to 
measure and compare the multiple scattering of positrons and electrons at 
energies in the neighbourhood of 0-4 Mev. The ratio of observed to 
theoretical mean absolute scattering angle was 0-99 for N, and 0-92 for A, 
with a standard error of 3°%%. Comparable results by other workers are 
discussed. In contrast with the conclusions of other workers, no significant 
difference in the behaviour of positrons and electrons could be detected. 


§1. INTRODUCTION 


Tue multiple scattering of charged particles has been the subject of much 
theoretical and practical study. A method of calculating the multiple 
scattering distribution function has been described by Bethe (1953) in a 
paper which contains also a discussion of the relation between the various 
derivations and formulae. Most multiple scattering investigations are 
now based on the cross section for single scattering by a Thomas—Fermi 
field first given by Moliere (1947). Simple formulae for mean value 
quantities of practical utility have been given by Scott (1952) and Scott’s 
formulae reproduce the numerical values given by Moliere’s multiple 
scattering theory (1948) with an error of 1%. 

The experiments fall into three main classes : (i) scattering of collimated 
electron beams by foils with a counter or ionization chamber as detector ; 
(ii) application of Fowler’s or a similar method to individual tracks in 
nuclear emulsions ; (iii) measurement of the scattering produced by the 
gas filling of cloud chambers (Groetzinger et al. 1950). Some of these 
experiments are reviewed by Corson and Hanson (1953) and on the whole 
there is agreement with the theory to within a few per cent. There are 
however several experimental results not fully explained, e.g. (i) the 
accurate experiment by Hanson et al. (1951) on 15-7 Mev electrons shows 
scattering larger than that predicted at the larger angles (Spencer and 
Blanchard 1954); (ii) an unexpectedly large difference between the 
scattering of electrons and positrons of about 1 Mev has been reported by 
Groetzinger et al. (1952) and discussed by Mohr (1954) ; (iii) electrons of 
0-59 Mev in nuclear emulsion show abnormally low scattering according to 
Hisdal (1952) . 
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The mean multiple scattering angle is inversely proportional to the 
product, pv, of the momentum and velocity of the scattered particles. 
The discrepancies just mentioned, having occurred with light particles of 
comparatively low energy involve low values of pv. The present paper 
describes an experiment with both positrons and electrons of the same low 
energy. An improvement in accuracy on previous cloud chamber work 
has been made by (i) avoiding the difficulties of magnetic curvature 
measurement by selecting and measuring the particle energy before entry 
into the chamber ; (ii) investigating in a subsidiary experiment the effect 
of track distortion due to turbulence ; (iii) decreasing the statistical errors, 
and (iv) taking special precautions in respect of comparative measurements 
on electrons and positrons. 


§2. APPARATUS 


The scattering medium was the gas of a conventional horizontal expan- 
sion chamber, 20 cm in diameter, operated at atmospheric pressure. 

The particle beam was provided by an iron-free solenoidal lens spectro- 
meter in which the emitting object was a disc, 3mm in diameter, of Cu®4 
which radiates both electrons and positrons. At the image position was an 
aperture through which the focussed beam passed to the cloud chamber 
along a brass tube, 16cm long, containing collimating slits and anti- 
scattering baffles. The spectrometer vacuum was separated from the 
cloud chamber by a 1-5 mg/cm? distrene foil window supported on a brass 
grid flush with the chamber wall. ‘The beam entered in the plane of the 
illuminated region of the chamber and was photographed with normal 
technique. 

The Cu®4 source was mounted 6mm from the lens axis and, on account of 
the opposite rotation of 90° given to the positron and electron beams by the 
lens field, it produced separate positron and electron images 12 mm apart. 
By rotating the source about the lens axis with the help of a Wilson seal 
either of the images could be brought over the exit hole and either of the 
beams sent into the chamber without altering any other part of the 
apparatus. This enabled alternate batches of ten pictures to be taken with 
electrons and positrons under identical conditions. 

The resolving power of the spectrometer, determined from the shape of 
the Thorium F-line, was such that the beams entering the chamber from 
the Cu®t had a roughly triangular Hp-spectrum with 4Hp/Hp=4%, 
where 4Hp is the width at half-height. The value of Hp was measured by 
reference to the F-line using for the latter the value 1388 gauss-cm given 
by Craig (1952). The momentum of the beams was kept constant through 
the experiment to about 1 in 10 000. 

To make sure that all the tracks measured actually came from the source, 
every tenth photograph was taken with the source covered with a shutter. 
100 such pictures showed no tracks that would have been included in the 
measurements. Further details of the apparatus are given elsewhere 


(Stott, 1955). 
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It was necessary to show to what extent the distortion of the tracks 
would influence the measurement of multiple scattering. On the basis of 
distortion measurements, made by a new method, of which the description 
and results are given elsewhere (Cusack and Stott 1955), it was concluded 
that distortion had a negligible effect on the measured mean scattering 
angle. Additional evidence for the absence of systematic chamber 
produced distortions is the fact that equal numbers of particles were 
scattered to the left and right and the algebraic means of the positron and 
electron scattering angles were very close to zero. 


Fig. 1 


CS nu 


Definition of the measurable quantities t, «, and «. The radius of the circle 
through ABC could also be measured. 


§3. TRACK SELECTION AND MEASUREMENT 


Figure 1 indicates the definition of several measurable quantities. After 
trials it was decided that the most accurate procedure was to measure 
t on the film, using a microscope with a calibrated eyepiece scale consisting 
of a line representing AC and graduations perpendicular to it at the mid- 
point. The magnification was chosen so that AC was equivalent to a 
10-0-+-0:01 cm chord in the chamber. The values of ¢ were converted to 
a’s for comparison with theory. 
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Every track containing a 10 cm chord was measured unless (i) it gave 
rise to a secondary electron of range >3 mm or (ii) it did not come from 
the beam entry window. As only one chord per track was used, all the 
scattering angles were independent and, as no movement of the microscope 
was involved, there was no ‘ stage noise’. Many tracks were, of course, 
longer than 10 cm and the rule was adopted that the chord should begin 
at the nearest clearly identifiable droplet to the entry window. Among 
the positron tracks in argon there occurred. five with very large values of « 
between 50° and 70° resulting from single scattering events with projected 
angles between 50° and 90°. Five such large single scatters is about the 
number to be expected in 4200 cm of track in argon, according to the 
projected Massey formula for positron scattering, and they may therefore 
be regarded as part of the single scattering tail of the distribution. On the 
other hand, one at least suggested by its appearance that it was a rather 
low energy track, perhaps a Compton electron. These five tracks are 
excluded from fig. 2 (b) and the tables. If they were included, | « ], 
for positrons only, would be 10-6, and, for all tracks, 10-4. The figure 
0-92 for argon in table 2 would become 0-945. This would not affect the 
conclusions to be drawn either for the absolute or the comparative 
measurements. 


§4. CALCULATION OF EXPECTED SCATTERING 


The momentum at the mid-point of the tracks was derived from the 
Hp focussed by the lens less a 3°/ correction for loss in the distrene 
window and chamber gas. The error in the corrected value should not 
exceed 0:5%. The momentum thus found is the peak momentum of the 
triangular spectrum referred to above. 2'NZ? was computed from the 
observed temperature and pressure allowing for the water and alcohol 
condensant which contributed about 5°% to the sum for nitrogen and 2%, 
for argon. Pod 

Using the formulae and notation of Scott (1952), the value of | « | was 
calculated from the equation 

(Ta |)? =x,22 (0:299-+ 0-269 logy» 2). 

To allow for electron—electron scattering it is usual to increase Z? by a 
factor (1+1/Z). U. Fano (1954) has recently given a more exact treatment 
of this correction and in the present case his formulae increase the value of 
| | by an additional 2% for N, and 0-6% for A. 

The mean scattering angle also requires correction for the ‘flat-chamber 
effect’ (Moliere 1948). The literature does not contain projected multiple 
scattering distributions for scattering media of finite depth, but it is clear 
that the effect of excluding tracks that reach large distances perpendicular 
to the observation plane will be to sharpen the distribution by reducing the 
large angle tail. “Moliere gives the mean scattering angle for the limiting 
case of an infinitesimally thin chamber. The conditions of the present 
experiment were chosen to approximate to the latter condition as closely 
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as practicable. A comparatively thin illuminated layer of 1 em was used. 
and long tracks, i.e. >10 cm, were selected. An elementary calculation, 
shows that only some 5%, of the incident particles survive to the end of the 
observation length without being scattered out of the light beam. It 
is believed therefore that Moliere’s ‘ flat chamber ’ scattering angle is the 
appropriate one for comparison, though it may slightly underestimate the 
scattering, since the extreme case is obviously not completely realized. 
Equations 10.4 and 10.4’ of Moliere (1948) give 


| | (flat chamber)=0-87 |« | (infinite chamber). 


All the relevant numerical quantities are given in table 1. 


Table 1 

Nitrogen Argon 
Hp (gauss-cm) 1929 2496 
Number of positrons 336 418 
Number of electrons 335 220 
[a] (pos.) 8:3 10-1 
Ja] (elec.) 9-0 10-0 
Total number of tracks 671 638 
ja] (all tracks) 8:6 10-1 
NZ? 2°56 x 10? 7-82 x 1074 
Q 155 198 
Xy 0-85 0-93 
Ja] (theor., flat ch.) 8-9 11-0 


The angles are in degrees; Q is the number of scatters in 10 em of track ; 
Q and x, are computed from eqns. (2) and (3) of Scott (1952). 


§5. ResuLTs AND CoMPARISON WITH THEORY 


The theories based on Moliere’s single scattering law do not differentiate 
between positrons and electrons because they neglect the spin-orbit 
effect, and the calculated value should therefore be compared with the 
results for positrons and electrons taken together. Fig. 2 shows the distri- 
bution of observed scattering angles and a gaussian curve of the same 
area and standard deviation as thesample. The gaussian fits satisfactorily 
for both nitrogen and argon and the deviations from it are such that the 
fit would be improved if higher terms of the Moliere distribution were 
used ; however, experiments of this sort are not capable of demonstrating 
the finer points of the distribution shape. A better test is the value of 
| ~ | (see table 1). The main error in | « | (obs.) is the statistical one of 3°%. 
Microscope setting errors are present in each reading but exert a negligible 
effect on the mean value, and 671 tracks in nitrogen, when measured 
independently by two observers, gave mean values agreeing to better than 
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0:-5°%. The main error in | « | (theor.) is probably due to the small depar- 
ture from the infinitesimal chamber mentioned above. In addition the 
value of the screening angle may have some slight inaccuracy because the 
field of the nitrogen molecule will not be represented exactly by the 
Thomas—Fermi field. The error arising from temperature and pressure 
measurements is unlikely to exceed 0-5%,. 


Fig. 2 
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The measured distribution of |«|. | «| is in degrees. 


For nitrogen, therefore, the agreement with Moliere’s and the equivalent 
theories, supplemented by Fano, is good. Taken by itself the argon 
result could also be called satisfactory agreement but table 2 shows that it 
goes some way towards confirming Hisdal’s result and suggests that the 
deficiency in scattering may be a real effect increasing with Z. As 
Hisdal has pointed out, most though not all experimental values of 
electron scattering are lower than the theoretical value. The most 
likely place to introduce small corrections into the theory is in the value of 
the screening angle. A recént paper by Lenz (1954) gives support to this 
and shows that the single scattering cross section for small scattering angles 
is considerably overestimated by the use of the Thomas-Fermi field because 
the latter does not give a good representation of the falling off of the 
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outer parts of the field. By choosing a more suitable representation, 
Lenz gave another cross section which accorded much better with the 
otherwise unexplained results of Bibermann ef al. (1949). However, 
these results refer to very small angles (~10-? to 10-3 radians) and the 
energies are only 50 to 100 kev. It seems therefore that such deficiencies 
of scattering, of the order of a few per cent, that cannot certainly be 
ascribed to experimental errors, can be qualitatively explained by the 
imperfection of the Thomas—Fermi field. Discrepancies in the scattering 
of electrons by beryllium have been attributed to this cause by Hanson 
et al. (loc. cit.). 


Table 2 


| 
Nitrogen | Argon |G5emulsion* 


Energy (Mev) 0-26 0-40 0-59 
Tal (obs.)/7a1 (theor.) 0:97 0-92 0-62 
Standard error 3% 3% 10% 


* Hisdal (1952). 


To compare positron and electron scattering, certain additional pre- 
cautions were taken to ensure that no conditions were changed other than 
the sign of the particles. This gave a reasonable chance of detecting the 
difference if it existed to the extent of 10° as reported by Groetzinger et 
al. These precautions were: (i) positron and electron pictures were 
taken in alternate batches of 10, taking advantage of the Cu®? source 
which emits both. The lens field was kept constant and any long period 
changes in the cloud chamber operation would affect both particles equally; 
(ii) the energy was such that the energy losses suffered in the window and 
chamber gas were the same for both beams according to Rohrlich and 
Carlson (1954) ; (iii) measurements were made by observers who did not 
know the sign of the particles. Scattering was slightly greater for posi- 
trons than for electrons in argon and, in nitrogen, |[« | was 8°% less for 
positrons than electrons. The significance of the latter difference may 
conveniently be estimated by applying the variance ratio test (F-test) 
which compares the widths of the « distributions. It shows that the 
probability of obtaining the observed positron and electron samples from 
the same distribution is greater than 12°. These results are at variance 
with the conclusions of Groetzinger et al., but the opportunities for statis- 
tical and systematic errors seem less in the present experiment and 
the best conclusion seems to be that the results are quite consistent 
with Mohr’s theory which states that at low Z the spin-orbit effect makes 
the positron r.m.s. scattering angle only 1° less than the electron one. 


To verify this accurately would require a considerable development of 
technique. 
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SUMMARY 

The numbers of neutrons and protons produced in the reactions 
WB(an)83N and 1B(ap)!8C at 1:51 Mev bombarding energy have been 
compared and the ratio of the reduced widths, y,?/y,?, found to be 
5:7--0-5. Between 1:4 and 2-4 Mev bombarding energy, the emission 
of protons is closely similar to that of neutrons, but y,?/y,? has fallen 
to a value nearer unity at a resonance at 2-16 Mev. 

The resonance previously reported when !!B is bombarded by alpha- 
particles of 1-58 Mev energy has been shown to be a doublet and a 
further resonance for neutron emission, but not proton emission, has 
been found at 2-05 Mev. 


$1. INTRODUCTION 

From experiments described in a recent paper (Shire, Wormald, Lindsay- 
Jones, Lunden and Stanley 1953) we found that the reduced width, y,?, 
from the reaction !°B(an)!83N at an alpha-particle energy of 1-51 Mev 
was about 15 times larger than the corresponding reduced width, y,?, 
for the production of protons from the mirror reaction !B(apo)C. In 
the same experiments we showed that the yield of neutrons varied with 
alpha-particle bombarding energy in a way closely similar to that of 
protons giving the ground state of C, except that a resonance for 
neutron production at 2:05 Mev bombarding energy was absent in the 
proton yield curve. 

We describe in the present paper a more careful determination of the 
ratio of the reduced widths y,? and y,? and an examination of the neutron 
and proton yields to higher bombarding energies than were possible at 
the time the previous paper was written. We have also obtained the 
1B(an)44N and !!B(ap)C yield curves between 1:4 and 2-4 Mev 
bombarding energy. 

References to earlier work on the bombardment of boron by alpha- 


particles are given in our previous paper and in Ajzenberg and Lauritsen 
(1955). 


* Communicated by the Authors. 
+ Now at Australian National University, Canberra, A.C.T., Australia. 
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§ 2. COMPARISON OF PRoTron AnD NEuTRON YIELDS 
AT 1-51 Mev AtpHaA-ParticLE ENERGY 
2.1. Hxperimental Method and Apparatus 
We used a beam of a few microamps of singly charged helium ions 
from our electrostatic generator and the separated isotope 1B on thick 
aluminium backings. Since 4B was found to be a copious emitter of 
neutrons near 1-5 and 2-0 Mev we carefully examined the 1B targets 
for 1B contamination and measured the amount of this isotope present. 
In our previous experiments the yield of protons was determined by 
one experiment and the yield of neutrons by another separate one. 
Substantial errors were possible in the measurement of the bombarding 
current, due to electron emission from the target and its surroundings, 
and to differences in the two targets. Even when the same target was 
used for both experiments there was no certainty that the target was 
uniform and unchanged or that exactly the same portions of a non- 
uniform target were employed in the two experiments. We have avoided 


Fig. 1 


K Mn:O,4 sotution 


Apparatus for comparing neutron and proton yields. For details see text. 


these sources of error in the present experiments by measuring the 
protons and neutrons simultaneously. The apparatus we used is shown 
diagrammatically in fig. 1. We supported the target T at the end of 
a long thin-walled brass tube A and, to avoid carbon deposits, heated it 
to a suitable temperature by the heater H before the beginning of an 
experiment. Alpha-particles entered through a defining aperture B ca 
protons from the target were counted by the gridded ionization chamber Cc? 
This we filled with argon at about one atmosphere pressure, at which 
pressure the ground state protons gave output pulses at least ten times 
the amplifier noise. Low energy protons arising from reactions leaving 
13(! in an excited state or from the 11B(«p)!4C reaction were excluded by 
the thin copper foil window W, which we soldered to the ionization 
chamber. Bevelled apertures E were placed at 5 cm intervals between 
the target and the ionization chamber to prevent protons from being 
scattered into the ionization chamber from the tube A. 
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We measured the neutrons by means of the Szilard—Chalmers reaction 
(Szilard and Chalmers 1934) produced in potassium permanganate 
solution contained in the spherical glass vessel, D (an acid carboy), 
surrounding the target. We used 48-7 litres of a solution of 40 grams 
of potassium permanganate per litre. This volume was large enough 
to make the correction for loss of neutrons from the surface reasonably 
small without making the apparatus, when full of solution, too difficult 
to handle. After a bombardment of the target, which took from 10 to 
30 minutes, we recovered the radioactive **MnO, produced by filtering 
the irradiated solution, and suitable washings of the vessel, through 
Whatman 54 filter paper, that had previously been impregnated with 
inactive MnO,. We then carefully dried the filter paper, encased it in 
cellophane, wrapped it round a thin-walled Geiger counter and measured 
the number of counts under standard conditions at a known time after 
the irradiation. We calibrated this detecting system by means of a 
standard RaBeF, source. (See § 2.2 (ii).) 

Corrections were made to the raw data to allow for the escape of 
neutrons from the potassium permanganate solution, for the dead time 
and secular variation of sensitivity of the Geiger counter, and for decay 
of the °*Mn during and after bombardment. We discuss these 
corrections and the accuracy of our measurements in the next section. 


2.2. Errors and Corrections 
. (i) Proton measurements 

The solid angle over which protons entered the ionization chamber 
was defined to better than 4% by one of two alternative apertures 
(0-998 cm and 0-674 cm diameter) placed just in front of the window W 
(fig. 1), and the mean angle between incident alpha-particle beam and 
counted protons was known to better than one degree. The total number 
of protons counted was more than 10000, so that the statistical error 
was less than 1%. We assumed that the angular distribution of the 
Po protons was exactly N= No (1—0-9 cos @) in centre-of-mass space, 
a distribution in close agreement with the experimental results given 
in our earlier paper and with a more recent measurement made with 
somewhat improved apparatus. With this assumption the only signifi- 
cant error arises from the measurement of the number of protons actually 
entering the ionization chamber. 

Figure 2 shows a typical pulse-height distribution obtained with a 
Hutchinson—Scarrott pulse-height analyser (Hutchinson and Scarrott 
1951). The * tail’ on the small-pulse side of the peak raises the question 
as to how much of it, if any, should be included in the proton count. 
These small pulses (about 10°, of the total in number) were not mainly 
due to protons scattered into the ionization chamber from defining 
apertures or other material outside it, for such scattered particles would 
have had less energy than unscattered protons and so would in many 
cases have produced more ionization and larger pulses than unscattered 
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protons. The energy of the unscattered protons was 4:1 Mey with a 
range in air or argon of about 24cm. The copper windows (two 
alternative thicknesses were used) were equivalent to about 4 and 6 cm 
of air, so that unscattered protons after passing through the window 
had a residual range of 18 to 20cm, of which only about 10cm was 
lost in passing through the chamber. Thus all scattered protons with 


Fig. 2 
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energies between about 2 Mev and 4-1 Mev would have given larger 
pulses than the unscattered protons, and the number of such pulses 
was well below 1°% of the number in the main peak. An experiment, 
in which we used a circular stop to prevent protons reaching the ionization 
chamber direct from the target, but which allowed protons scattered 
from the baffles E (fig. 1) to be counted, showed that the number of 
protons scattered into the ionization was less than 1%. 
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We can account for the tail by assuming that it was due to protons 
scattered by the copper window which had hit the walls of the ionization 
chamber without traversing the full 10cm path necessary to give a 
pulse of normal size. A rough calculation showed that this explanation 
gave the right number of small pulses within a factor of two or three. 

- We have therefore taken the numbers indicated by our counting equipment 
as the true number of protons and made no correction for the tail. We 
have not, however, included the portion of the tail consisting of pulses 
smaller than the setting of our scaler discriminator. The correction 
would be only 1 or 2° at most and its omission will tend to compensate 
for the small number of protons which may in fact be scattered into 
the ionization chamber by the apertures E. 


(ii) Neutrons 


The corrections applied for the dead time of the Geiger counter and 
for the decay of the radioactive manganese during and after irradiation 
were straightforward and need not be discussed in detail. We checked 
the sensitivity of the Geiger counter before and after each run by means 
of a standard y-ray source. Errors arising from these corrections appear 
with other random errors in the neutron—proton plot in fig. 4. In an 
experiment carried out with a blank aluminium target we found no 
neutrons, so that no errors arose due to neutrons from the target backing 
or from the apparatus. A comparison of our neutron yield curve with 
that of Jones and Wilkinson (1953) for the !%C(«an)!®O reaction showed 
that our results are not significantly affected by carbon contamination. 

The three possible sources of substantial systematic error arise from 
the estimation of the correction to be applied to neutrons escaping from 
the permanganate solution, from the uncertainty in the value of the 
standard RaBeF, neutron source, and from the presence of “B in the 
0B target. 

In fig. 3 we plot the fraction of neutrons captured against mean neutron 
energy, calculated by using Wallace’s elementary diffusion theory 
(Wallace 1949) and Marshak’s values for the slowing down length 
(Marshak 1947). This theory, which is only approximate, predicts too 
great a loss of neutrons. We have, therefore, compared experimentally 
the loss of neutrons from a RaBeF, source and from a 24Na—D source 
to obtain a value of 0-78-+-0-05 (A, fig. 3) for the fraction of neutrons 
of 4:5 Mev mean energy captured in our experimental vessel. The 
fractions of neutrons captured from the !°B(an)!®N and the 11B(an)14N 
reactions must then be close to the values shown at B and (, fig. 3. 
We adopt values of 0-89--0-05 and 0-96-+0-04 where the errors allow 
for both the estimated error of the determination of point A and the 
errors of interpolation. 

The preparation and calibration of our standard RaBeF, source has 
been described by Bretscher, Cook, Martin and Wilkinson (1949). It has 
recently been recalibrated in terms of the 24Na y-ray photo-disintegration — 
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‘cross section of deuterium and found to give 7:75 x 10*-+4°%, neutrons 
per second.* 

The yield of neutrons from the “B(«n)!4N reaction is substantial at 
1-51 Mev alpha-particle energy (see fig. 5(d)). We measured the amount 
of 11B present in our !°B target by three methods: (a) the “B(px)&Be 
reaction at 0-7 Mev bombarding energy, (b) the 'B(ap)!4C reaction 
at 1:58 Mev, and (c) the B(an)!4N reaction at 2:05 Mev. We used a 
magnet to analyse the alpha-particles emitted in (a) and to separate the 
protons from 14B from those of !°B in (b). Method (c) assumes that the 
peak observed in the 1B target at a bombarding energy of 2-05 mev 
is in fact due entirely to the presence of “B. Our reasons for this 
assumption are given in §3.2(i). The three methods gave values of 


Fig. 3 
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One 2 3 4 5 MeV 
Loss of neutrons from tank. 
Theoretical curve for tank radius 23 cm. 
—---- Curve assumed for calculating corrections. 
A, experimental point for neutrons from RaBeF, source. 
B, C, points assumed for neutrons from 1°B(an)'®N and '1B(an)!4N reactions. 


0-07-++0-02, 0:06+0-02 and 0:08-0:02 for the number of ™B atoms 
in the !°B target relative to the number of ''B atoms in the 'B target. 
We have used the mean value 0-07-+0-01 for this ratio. We measured 
the ratio of the number of neutrons from the same !B and ''B targets 
at 1:51 Mev bombarding energy and also the angular distributions of 
the neutrons (see § 3). From these results we calculated that the number 
of !!B(an)4N neutrons of 1-23 Mev energy entering our permanganate 
solution was 3-6+1-2°% of the number of '°B(an)!®N neutrons. If we 
ignore, as we reasonably may for a correction term, the difference in 
the capture fractions for the two groups of neutrons, the number of 
neutrons measured by means of the permanganate solution must be 
multiplied by 1—0-036, or 0-964-+-0-012, to allow for the presence of 
11B atoms in the 1°B target. 


* D. H. Wilkinson. Private communication. 
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2.3. Hxuperimental Results 


The results of six measurements of the protons and neutrons from 1B 
are shown in fig. 4, before corrections have been applied for the loss of 
neutrons from the permanganate solution or the presence of 1B. 
The mean value of the ratio of neutrons to protons is 2:71-0:05 where 
the results have been weighted in proportion to the number of proton 
counts and the errors shown do not include the errors discussed in 
§ 2.2 (ii). Applying now the corrections for these errors we obtain the 
final corrected value of 2-27--0-21 for the ratio of the number of protons 


Fig. 4 
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Total number of protons and neutrons emitted at 1-51 Mev alpha-particle 
energy. The errors shown are the statistical errors arising from the 
number of particles counted. 
from the reaction !°B(«py)C at H,=1-51 Mev to the number of neutrons 
from the reaction B(an)N at the same bombarding energy. If we 
use for the reaction radii the value 1-45(131/3+-11/3) 10-13 em, ice. 
4-86 x 10- cm, we obtain for y,?/y,? the value 5-7-L0-5. 


§ 3. Resonances, YieLD CURVES AND ANGULAR DISTRIBUTIONS 
3.1. Hxperimental Method and Apparatus 


To count protons we used two ionization chambers of the kind described 
in § 2.1. One was usually left fixed at approximately 90° to the alpha- 
particle beam and the other could be moved to any angle between 9° 
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and 171° to the beam without disturbing the vacuum. By choosing a 
suitable thickness for the copper windows we excluded scattered alpha- 
particles and, in the case of the !°B(ap)!8C reaction, protons other than 
those leaving the residual !°C nucleus unexcited. 


Fig. 5 
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Proton and neutron yields. 
The vertical scales of (a) and (b) have been adjusted to make the ratio of neutrons 
- to protons at 1-51 Mev bombarding energy equal to 2:27. (See § 2.3.) 


Neutrons were measured by a ZnS-plastic scintillator and photo- 
multiplier of the type described by Hornyak (1952). The data eee 
by Hornyak show that the number of counts for a cae ges: 
setting is fairly closely proportional to the energy of the neutrons. ) 
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have assumed exact proportionality in estimating the number of neutrons. 


from our observed counts. For angular distribution measurements we 


used two scintillators, one fixed near 90° to the incident beam, the other 


movable between 0° and 135°. To obtain a reasonable counting rate we 
placed the scintillators close to the boron targets, and we have applied 
corrections for the finite solid angle of the movable scintillator in the 
angular distribution experiments. To compare the angular distribution 
of the neutrons from the !B(an)!8N reaction with theory, we first 
computed the angular distribution to be expected in laboratory space 
for an assumed centre-of-mass distribution and then applied a finite- 
counter correction by the method proposed by Breitenberger (1953), 
using the tables given by Price (1954) to express the point-counter 
laboratory distribution in terms of Legendre polynomials. For other 
distributions we corrected for finite counter size by the method given 
by Deutsch (1951). 


3.2. Experimental Results 
(i) Yeeld Curves 

We give in figs. 5(a), 5(b), 5(c) and 5(d) our yield curves for the reactions 
WB(apy)8C, B(an)18N, 4 B(«p)4C and 14B(an)14N. The 1°B results at 
2:16 Mev bombarding energy lead to values of y,,?/y,? of 1-1 or 1-6 according 
as | for the outgoing particle is taken to be 1 or 2. The y-ray yield which 
we obtained (not shown) for the 1°B(«, p, y)!8C reaction agreed with the 
results of Talbott and Heydenburg (1953). 

A neutron resonance at 2-05 Mev for the !°B(«n)!8C reaction, reported 
in our previous paper, is almost certainly due to the 11B(an)4N reaction 
for which there is a very strong peak at 2:05 Mev. We used three 
different targets of 1°B and found the height of the 2-05 Mev peak 
proportional to the amount of 14B contamination present. The resonance 
which we previously reported at 1:58 Mev for the !B(«p)#C reaction 
consists of two resonances at 1-51 and 1:58 Mev. These can be more 
clearly seen in the results of the 14B(an)!4N reaction (figs. 5(c) and 5(d)). 
The resonances observed for bombarding energies of 1-51, 1:58 and 
2-05 Mev correspond to levels in 1°N at 12-096, 12-147 and 12-491 Mev 
previously found by Roseborough, McCue, Preston and Goodman (1951) 
by means of the “C(pn)!4N reaction. 


(ii) The Angular Distribution of the Neutrons from the Reaction B(an)3N 

We show in fig. 6 the experimental points found at 1-51 Mev bom- 
barding energy together with a theoretical curve obtained from an 
assumed distribution 1—0-9 cos? @ in centre-of-mass space, converted 
to laboratory space and adjusted for the finite detector used. 


(ili) Other Angular Distributions 
At 2-16Mev the protons from !B(«p))"®C and the neutrons from 


1B(an)!83N have similar distributions. Both could be fitted in centre- 
of-mass space to 1—(0-5-L0-1) cos? @. 


. 
— 
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We measured the neutron distribution for the 11B(an)4N reaction at 
1-51 Mev bombarding energy, as part of the determination of the correction 
for contamination of our !°B targets by 11B, and in a preliminary way 
at 1-58 Mev and 2-05 Mey. The strong forward and backward intensities 
observed at 1-51 Mev suggest that the 12-096 mev level in 15N has 
spin 5/2 and even parity, an assignment tentatively made by Hinchey, 


Fig. 6 
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Angular distribution of neutrons from 1B at 1-51 Mev alpha-particle energy. 


The curve gives the theoretical angular distribution corresponding to 
1—0-9 cos? 6 in centre-of-mass space. 


Stelson and Preston (1952) from measurements of the total neutron 
cross section of nitrogen. The distribution at 1:58 Mev bombarding 
energy is nearly isotropic, which suggests that the °N level at 12-147 Mev 
has spin 3/2 and odd parity, but our present results could be fitted by 
other assignments if suitable choices of spin-channel intensities are made. 
The strong forward intensity at 2-05 Mev suggests that the 1°N level at 
12-495 mev has spin 5/2 and even parity, but our measurements indicate 
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that the backward and forward intensities (in centre-of-mass space) are 
not equal, so that interference effects may be occurring. We intend to 
the measure proton angular distributions from the 1+B(«p)“C reaction 
for these levels, since these distributions should determine the spins and 
parities concerned unambiguously. 


§ 4. Discussion AND CONCLUSION 

The main conclusion that we may draw from our results is that the 
emission of protons in the 1B(«py)C reaction and the emission of 
neutrons in the mirror reaction !°B(an)!8N are strikingly similar. The 
yield curves for the two reactions (figs. 5(a) and 5(6)) show similar 
resonances for all bombarding energies between 1-4 and 2-4Mmev; the 
angular distributions of protons and neutrons fit with a high degree of 
accuracy in centre-of-mass space the relation J,=I,(1—0-9 cos? @) for a 
bombarding energy of 1:51 Mey and are similar and perhaps identical 
at higher energies. Nevertheless, the ratio of the reduced widths y,?/y,? 
obtained from our results at 1:51 Mev bombarding energy, though lower 
than the estimate given in our earlier paper (Shire et a/. 1953), is substan- 
tially greater than the value of unity which might be expected, if charge 
independence holds, for mirror reactions. The discrepancy appears to 
be greater than can easily be accounted for by modification of the 
usually accepted formulae for barrier penetration (Bloch ef al. 1951, 
Blatt and Weisskopf 1952), for instance along the lines proposed by Scott 
(1954). Nor can the difference be accounted for by taking into consider- 
ation values of / for the outgoing particles greater than the minimum 
value of two ; higher values make a change of only a few per cent in the 
value of y,?/y,?. 

The inequality of y,? and y,? is, however, not surprising if the compound 
nucleus is not in a pure isotopic spin state. Isotopic spin mixing is to 
be expected when states of the same spin and parity are near to one 
another and this is likely to be so in the present case where many levels 
come close together. 

The behaviour of 'B is quite different from that of 1°B. Though the 
double resonance near 1-58 Mev bombarding energy appears for both 
protons and neutrons, protons are almost completely absent for the very 
strong neutron resonance at 2-05 Mey. Such difference is, of course, not 
surprising, as the reactions 'B(«p)@C and 1B(an)!N are in no sense 
mirror reactions. 
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SUMMARY 


A comparison is made between the nuclear disintegrations produced 
in carbon, nitrogen and oxygen by 140 Mev protons, and the nature of 
the disintegrations expected from the Goldberger model (1948). 


§ 1. INTRODUCTION 


THERE have been two previous communications from this laboratory 
relating to the nuclear interactions produced in photographic emulsions 
by 140 mev protons. Lees et al. (1953), referred to as I, reported on 
measurements of the cross sections for the production of nuclear dis- 
integrations by protons with an average kinetic energy of 130 Mev in 
both the heavy (silver and bromine) and light (carbon, nitrogen and 
oxygen) constituents of nuclear emulsions. The results were shown to 
be in agreement with the ‘ optical model’ of Fernbach, Serber and 
Taylor (1949). 

Morrison, Muirhead and Rosser (1953), referred to as II, reported on 
an experimental study of the nuclear disintegrations produced in silver 
and bromine by protons with an average kinetic energy of 140 Mev. 
The results obtained were shown to be consistent with the assumption 
that the disintegrations proceeded in two main stages. During the 
first stage the incident proton initiated a nucleon cascade leading to the 
direct emission of energetic nucleons (Goldberger 1948). The nucleus 
was then left in an excited state and nuclear particles were subsequently 
evaporated from it (Weisskopf 1937) until the residual nucleus reached 
its ground state. 

In this paper a description will be given of the disintegrations observed 
in nuclear emulsions when protons with an average kinetic energy of 
140 Mev strike the light nuclei carbon, nitrogen and oxygen. The 


experimental results will be compared with a theoretical model similar 
to that used in paper IT. 


* Communicated by Professor P. I. Dee, F.R.S. 
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§ 2. EXPERIMENTAL CRITERION FOR THE SELECTION OF STARS 
PRODUCED IN Ligut Nucuer 


The stars were found by ‘area scanning’ in Ilford G5 emulsions 
which had been exposed to the external beam of the Harwell cyclotron. 
The protons entering the emulsions possessed a kinetic energy of 146 Mev. 

The stars selected for analysis, as representing the disintegration of a 
light nucleus, were those in which an «-particle with a kinetic energy of 
less than 9 Mev had been emitted (Menon et al. 1950, Hodgson 1954). 
It will be shown in § 4, that the stars selected using this criterion do 
not represent all the modes of disintegration of a light nucleus. 


§ 3. THe THEORETICAL MODEL 


3.1. Procedure 


Our theoretical model assumed that the incident protons possessed a 
kinetic energy of 140 Mev and that they initiated nucleon cascades in 
the light nuclei leading to the direct emission of knock-on nucleons. 
The residual nuclei were then left in excited states, and at a later stage 
de-excitation occurred. 

The procedure for obtaining numerical data on the products of a 
nucleon cascade, and on the excitation energy of the residual nucleus 
was the same as that outlined for the heavy nucleus in paper II. This 
procedure will be referred to as the Goldberger method. 

The struck nucleus was assumed to consist of 7 protons and 7 neutrons. 
and to possess a radius of 3-4 10-13 cm; it was also assumed that the 
nucleus possessed an effective coulomb barrier of 2 Mev. 


3.2. Limits of Application 

In carrying out the calculations we assumed that only single 
nucleon-nucleon collisions took place at all energies. However, 
evidence will be offered in § 6, from which it is necessary to conclude 
that the knock-on nucleons may occasionally collide with nuclear 
subgroups, such as «-particles, inside light nuclei. It will be shown 
that the effect is a small one, so that the use of single nucleon-nucleon 
collisions for the Goldberger method is not likely to be grossly incorrect. 

In performing the calculations the possible effects of the wide spacing 
of the energy levels in light nuclei were neglected initially. However, 
the results of Titus and Strauch (1954) have shown that the wide spacings 
of nuclear energy levels in light nuclei play an important role when 
protons of 96 Mev are inelastically scattered by light nuclei. This. 
feature will be discussed more fully in § 7. . : 

Despite these inadequacies it will be shown in subsequent sections 
that the model gives a satisfactory explanation of many of the modes 
of disintegration of a light nucleus when struck by a proton with an 


energy of 140 Mev. 
ZY 2 
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§ 4. Tue Cross Section FoR NUCLEAR INTERACTION 


The total cross section, calculated by us, for the production of nuclear 
interactions in light nuclei (referred to henceforth as stars’) was 
(21-+43)10-26 cm?. This is in very good agreement with the experimental 
value of (16--7)10-26 em2, presented in paper I, and with the value of 
(22+2)10-26 em? obtained by Cassels and Lawson (1954), using counter 
techniques and 134 Mev protons incident on carbon. 

The size-frequency distribution for stars produced in light nuclei, 
presented in paper I, is reproduced in table 1, column 2. This was 
obtained experimentally from a comparison of stars produced in normal 
and diluted Ilford G5 emulsions. Unfortunately this method does not 
allow us to decide whether any individual star was produced in a heavy 
or light nucleus. 


Table 1 


Cross section x 102° em-? 


Type Experimental 
of star Th cnt 
Pies Stars with pas 
P eee Mev 

Op 2-6 1-5 — | 
Ip B45 0-1 
2p EW tas: 0:3 a 
3p 342 3-0 
4p rigse a) ah 6-0 
5p 3:5+1 9+4 4-0 > 11-6 ~ 12 
6p 1540-5 J ey See 

Total 16+7 eal koi > 2143 


In the present work if an «-particle with an energy less than 9 Mev 
was observed to be emitted from a star, then that star was attributed 
to the disintegration of a light nucleus. The size-frequency distribution 
obtained with this criterion is also shown in table 1, column 3. The 
value of 32:5-++-2-5 cm obtained in paper I for the mean free path for 
star production in Ilford G5 emulsions, was used to determine the 
absolute cross sections. It can be seen that with our criterion we select 
most of the 3, 4, 5 and 6 prong stars produced in the light nuclei, but 
that many of the 0, 1 and 2 prong stars are not selected. 

As there is no reliable statistical theory available for the processes 
following the nucleon cascade in light nuclei, it is not possible to 
calculate the size-frequency distribution from the Goldberger model, 
but it is possible to make a qualitative estimate. 
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The distribution of the excitation energy of the residual nucleus 
obtained by the Goldberger method is shown in fig. 1; the average 
value for this quantity was 23 Mev. For high excitation energies the 
nucleus is expected to break up completely, but at low excitation energies 
the residual nucleus need emit only one charged particle or possibly a 
y-ray before losing its excitation energy. The thresholds for the nuclear 
reactions “C(p; p, 3x), “4N(p; n, 2p, 3«) and 1*O(p; p, 4a) are 7-3, 
19-8 and 14-5 Mev respectively. Therefore at excitation energies below 
about 15 Mev very few nuclei will break up completely whereas at higher 
energies the probability for these processes is high. If it were assumed 
arbitrarily that for excitation energies greater than 15 Mev the light 
nucleus broke up completely, giving a 4p, 5p or 6p star and that for 
_ lower excitation energies it did not, then we would expect the cross 
section for stars of (0+1+2-+3) prongs to be about 9x 10~?° em?, and 
about 12x 10-*6 cm? for stars of (4+5+6) prongs. These values are in 
satisfactory agreement with those shown in column 2 of table 1. The 
Goldberger model can, therefore, account qualitatively for the number 
of small and large stars observed experimentally. 


Fig. 1 
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Calculated distribution of excitation energy in residual nucleus 
following the passage of the nucleon cascade. 


A comparison of columns 3 and 4 of table 1 also suggests that the 
application of the criterion discussed in § 2 favours the selection of stars 
with excitation energy greater than ~15 Mev. Further evidence in 
support of this conclusion will be offered in the next section. 
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§ 5. COMPARISON OF THE CALCULATED AND EXPERIMENTAL Data - 
For Protons EmirreD WITH ENERGIES GREATER THAN 30 Mev 


In table 2 the experimental energy distribution of protons emitted 
with energies >30 Mev is compared with the predictions of the Goldberger 
model. The energies of the primary and secondary protons were 
determined by grain counting. The experimental results are based on 
71 stars, all of which had primary protons with energy greater than 
130 Mev, their mean energy being close to 140 Mev; they also had at 
least one «-particle <9 Mev. 

Table 2 


Energy in Mev 30-55 | 55-75 75-95 95-115 >115 


No. of protons :— 
experimental 19 


Goldberger model 
excitation energy 22 14 
>15 Mev 


Goldberger model 
excitation energy 
<15 Mev 


Our calculated results are shown in table 2, and are based on 118 
‘stars’, the data for 59 ‘stars’ being doubled.. Of these 68 had 
excitation energies greater than 15 Mev. An examination of table 2 
shows satisfactory agreement between the experimental data, and the 
energy spectrum calculated for nuclei which are left with excitation 
energy greater than 15 Mev. This agreement lends further support to 
the suggestion made at the end of § 4 that the stars selected experimentally 
by the presence of an «-particle with energy <9 Mev are produced by 
a nuclear excitation greater than ~15 Mev. 


Table 3 


40-60 | 60-80 | 80-100 | 100-120} +120 


Angle in degrees 0-20 | 20-40 


No. of protons :— 
experimental 11 22 7 8 3 1 l 


Goldberger model 
excitation energy 12 24 8 2 — 2 2 
>15 Mev 


Goldberger model 
excitation energy 12 12 + + 
<15 Mev 
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It was suggested in paper II that the angular distribution obtained 
by the Goldberger method should be compared with the distribution of 
the actual angle in space between the primary and secondary protons. 
The calculated angular distribution is compared with the experimental 
results in table 3. It can be seen that the agreement is satisfactory. 


3 6. COMPARISON OF THE CALCULATED AND EXPERIMENTAL Data _ 
FOR THE CHARGED PARTICLES EMITTED WITH ENERGIES LESS THAN 30 Mev 


The majority of the tracks of particles with energy less than 30 Mev 
are emitted in the processes following the nucleon cascade. There is at 
present no reliable statistical theory for the processes occurring in light 
nuclei during the second stage of the disintegration so that no exact 
comparison can be made between the calculated and experimental 
results. The experimental results will be discussed in the light of the 
predictions of the Goldberger type calculations. 

It was found experimentally that there was little variation in either 
the average star size or in the characteristics of the stars produced in 
light nuclei when the energy of the primary protons varied within 
reasonable limits. In this section results are included from all stars 
from which an «-particle <9 Mev had been emitted, and whose primary 
protons had energies greater than 110 Mev. The total number of stars 
examined was 176. 


6.1. Results for «-particles 


The experimental value for the ratio of the number of «-particles to 
the total number of prongs in a star was 0-48+0-03 for the 176 stars 
produced in light nuclei. This high value reflects the ease with which 
light nuclei break up into «-particles at comparatively low excitation 
energies. The value obtained is consistent with the results obtained 
from stars produced by z-mesons and the cosmic radiation (Menon ¢¢ al. 
1950, Harding 1949). 

The measured energy spectrum for the «-particles emitted from stars 
produced in light nuclei is shown in fig. 2. 

The experimental ratio of the number of «-particles emitted in the 
forward direction in the laboratory system to the number emitted in 
the backward direction is shown in table 4 as a function of energy. 


Table 4 

Energy. 0-2 | 24 | 46 | 68 | 8-10 | >10 | Total 
in Mev 

Forward ise) 2°3 2-5 2:3 3°6 2°9 2-4 

to Experimental 

backward. +0-6| 40-5} +0-6| +0-7| +1-4| 41-2] +0°3 
ratio for $$$} $$} $$ |__| 
a-particles| Calculated 2-1 2-4 1-9 1-6 1:5] 1:35 | — 
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The forward excess observed may be attributed mainly to the forward 
momentum the nucleus normally possesses after the residual nucleon 
cascade is over. The average forward velocity of the nucleus calculated 
by the Goldberger method for ‘ stars’ with excitation energies greater 
than 15 Mev was 4x 108 cm/sec. Using this value the forward to the 
backward ratio was calculated assuming that all the «-particles were 


Fig. 2 


NUMBER OF cc PARTICLES /2 MeV INTERVAL 


ENEACY iN Mev 


Measured energy spectrum for «-particles. 


emitted isotropically in the coordinate system of the recoiling nucleus 
during the process following the nucleon cascade. The calculated results 
are also shown in table 4. It can be seen that the forward excess for low 
energy «-particles can be accounted for in this manner, but it is difficult 
to explain the forward excess at higher energies in a similar fashion. 
For example, the forward to backward ratio for all «-particles greater 
than 8 Mev is observed to be 3-2+0-9, whereas the calculated value 
is 1-4. There is some uncertainty in the model used in the Goldberger 
type calculations, but the number and angular distribution of high 
energy protons experimentally observed to leave the nucleus show that, 
after the nucleon cascade is over, the velocity of the recoiling nucleus 
cannot be considerably greater than the calculated value. 

The forward excess could be explained if either the «-particles were 
emitted preferentially in the forward direction relative to the recoiling 
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nucleus during the processes following the nucleon cascade, or if some 
a-particles were emitted in direct collision processes during the nucleon 
cascade. The latter view is supported by the occasional observation of 
a-particles ejected with energies >30 Mev in identified disintegrations 
of light nuclei. These «-particles could arise either from a direct collision 
between a nucleon and an «-particle inside the nucleus, or a ‘ pick-up ’ 
process during the nucleon cascade (Chew and Goldberger 1950). 

One of these two events was of particular interest in that an «-particle 
of energy 43 Mev was emitted at an angle of 26° with respect to a primary 
proton of measured energy 140 Mev, and a proton of energy 65-+20 Mev 
was emitted at an angle of 134° with respect to the primary particle. 
The three tracks were coplanar, the angle between the high energy 
x-particle and the secondary proton being 160°. The only other charged 
particle emitted was an «-particle of energy 2-5 Mev which was probably 
due to the decay of the residual nucleus. The dynamics of the problem 
suggest that the 65+20Mmev proton represents the primary particle 
recoiling after colliding with an «-particle subgroup inside the nucleus. 
If the high energy «-particle were due to a ‘ pick-up’ process, then the 
dynamics of the collision should correspond to a single nucleon-nucleon 
collision, which would be unlikely to give a proton of energy 65 Mev 
in the backward direction even after allowing for possible double 
scattering. One other similar event was observed in which an «-particle 
of energy 26 Mev was emitted at an angle of 38° with respect to a primary 
proton of energy 140 Mev. A proton of energy 55+-10 Mev was emitted 
at an angle of 99° with respect to the primary. As no other tracks were 
emitted, this event was classified as a disintegration of a heavy nucleus. 

We, therefore, conclude tentatively that the primary protons and 
knock-on nucleons sometimes collide with «-particle sub-groups inside 
light nuclei, although part of the forward excess may be due to a 
‘pick-up’ process. An average of ~0-25 «-particles per star, emitted 
in the forward direction during the nucleon cascade, is adequate to. 
account for our experimental results. This number is much lower than 
the average of one knock-on proton and one knock-on neutron per star 
emitted during the nucleon cascade. We need, therefore, assume that 
only a small proportion of the collisions inside light nuclei involve: 
a-particles and the use of only single nucleon-nucleon collisions in the 
Goldberger type calculations is not likely to be grossly incorrect. 


6.2. Results for Protons of Energies less than 30 mev 


The measured energy distribution of protons emerging from 176 stars 
produced in the light nuclei is shown in fig. 3. The shaded area represents 
the calculated number of knock-on protons expected on the Goldberger 
model. These results were obtained from ‘stars’ with calculated 
excitation energies greater than 15 Mev; the results were normalized 
to a total of 176 stars. It can be seen that the broad tail in the energy 
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spectrum can be accounted for as due mainly to knock-on protons emitted 
during the nucleon cascade. The statistical errors on the number of 
protons having energies between 20 and 30 Mey and ending in the 
emulsion are large, but the existence of these protons was confirmed by 
the observation of protons in this energy range by grain counting in 


underdeveloped plates. 
Fig. 3 
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Measured energy spectrum for protons. Shaded area represents 
the calculated contribution to the spectrum from knock-on protons. 


The experimental forward to backward ratio for all identified protons 
of energies less than 30 Mev was 2-3-L0-3. If we assumed that there 
were no knock-on protons with energies less than 30 Mev the velocity 
of the residual nucleus after the passage of the nucleon cascade would 
be increased from 4 108 cm/sec to 5-2 x 108 cm/sec, using the Goldberger 
model. The forward to backward ratio calculated assuming that all 
protons of energy below 30 Mev were emitted isotropically relative to 
a residual nucleus which had a velocity of 5-2 10% em/see was 1-4. 
It can be seen that the forward excess observed experimentally cannot 
be due only to the motion of the recoiling nucleus. If, on the other hand, 
it were assumed that the number of knock-on protons emitted in the 
176 stars produced in the light nucleus is as shown in fig. 3 and that 
they had the angular distribution given by the Goldberger model, and 
that all the other protons were emitted isotropically relative to a residual 
nucleus which had a velocity of 4108 cm/sec, then the forward to 
backward ratio should be 2-4. This value is in good agreement with the 
experimental figure of 2-3+-0-3, showing that the Goldberger model can 
account satisfactorily for the forward excess of protons <30 Mev observed 
experimentally. 
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§ 7. DiscussIon 


The results obtained in this paper lend further support to the suggestion 
made in paper IT that nucleons of low energy (<30 mev) do emerge from 
nuclei during a nucleon cascade, and that they must, therefore, have a 
finite mean free path in nuclear matter before being incorporated in a 
‘compound nucleus. The work of Bernardini, Booth and Lindenbaum 
(1952) and of Hodgson (1954) also supports this assumption. 

Our theoretical model is able to account satisfactorily for the total 
cross section for nuclear interaction in the light nuclei, and qualitatively 
for the size-frequency distribution (§ 4). The model also gives an adequate 
description of the modes of disintegration of the light nuclei when the 
nuclear excitation is greater than ~15 mev (§§ 5 and 6). 

When small quantities of energy are transferred to the nuclei by the 
primary protons, however, the subsequent de-excitation is frequently 
realized by y-emission or by the ejection of a single nucleon. Interactions 
of this type cannot be attributed to the light nuclei when observed in 
a nuclear emulsion. Interactions involving low nuclear excitation may 
be examined, however, in experiments using other techniques, and the 
results obtained show that in thisregion our theoretical model isinadequate. 

The theoretical model is based on the assumption that the nucleon 
cascade proceeds through a series of nucleon-nucleon collisions. We 
have also assumed that, with the exception of the operation of the 
Pauli exclusion principle, these collisions take place in the same manner 
as those between free nucleons. Using these parameters our model 
cannot account for the experimental results of Titus and Strauch (1954), 
who have shown that when 96 Mev protons are incident upon light 
nuclei there are peaks in the energy spectrum of the inelastically 
scattered protons. These peaks occur at energies equal to the primary 
energy minus the excitation energies of known nuclear levels. For 
example, peaks were observed at 4:3, 10-0 and 20-9 mev below the elastic 
peak with carbon. Our theoretical model would predict that the energy 
of the knock-on protons should be continuous. 

It might be expected that the effect of the discrete spacing of energy 
levels in light nuclei would be to forbid all the collisions of the primary 
proton, except when the energy transferred in the collision corresponded 
to the excitation of a nuclear energy level. On our model this would 
reduce the number of allowed collisions for the primary proton, and 
therefore the calculated cross section for star production. The widths 
of the nuclear levels of low excitation energy in light nuclei are generally 
much narrower than the average distance between levels, and the density 
of levels only becomes appreciable at excitation energies of about 25 Mev. 
The Goldberger type calculations were, therefore, repeated assuming that 
all collisions for which the energy transferred from the primary proton 
was less than 25 or greater than 115 Mev were forbidden. It was found 
that the calculated cross section for star production was reduced from 
{21-+3)10-%6 cm? to (10-2)10~* cm?. This value is significantly smaller 
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than the figure of (22+2)10-*cm? determined experimentally by 
Cassels and Lawson (1954) for carbon. Furthermore, with this model 
the number of collisions with the correct energy transfers is not enough 
to account for the number of protons present in the peaks observed by 
Titus and Strauch. We are, therefore, forced to conclude that the 
assumption of single nucleon-nucleon collisions using the parameters 
observed for free collisions does not hold exactly for light nuclei. It is 
possible that the distribution of scattering angles in the centre of 
momentum system of the colliding nucleons may be modified in the 
presence of the potential well in light nuclei such that energy transfers 
corresponding to excited levels of the nucleus are preferred, whilst the 
total cross section for collision remains about the same as that expected 
using our original model. If, on the other hand, the use of the parameters 
for free collisions were correct as a first approximation and collisions 
leading to energy transfers between levels were allowed, then the peaks 
observed by Titus and Strauch would have to be due to a different type 
of collision that occurred occasionally, possibly involving the collision 
of the primary proton with a nuclear subgroup or possibly even the 
nucleus as a whole. 

The present work does not allow us to decide whether the excited 
nucleus, left after the passage of the nucleon cascade, sometimes breaks 
up completely in one act, or always emits single particles or photons 
successively until it has either lost all its excitation energy, or broken 
up completely into its constituents: The work of Livesey and Smith 
(1953) and Goward and Wilkins (1953) on the photodisintegration of 
carbon and oxygen at y-ray energies below 25 Mev suggests that the 
latter view is probably the correct one. 
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ABSTRACT 


The positron emitters 7-4 sec 2>Al and 6:6 sec 26™Al have been produced 
by proton capture in *Mg and Mg. Their subsequent decays have been 
investigated by means of a magnetic lens spectrometer. The maximum 
energies were measured as 3-24-+0-03 mev for %Al and 3-21-+0-03 mev 
for "Al. The spectra show no complexities and yield straight Fermi 
plots down to about 1300 kev. The deviations from linearity at lower 
energies are probably due to scattering in the spectrometer. No gamma- 
rays except annihilation radiation was found. 


§ 1. EXPERIMENTAL 


THE maximum energies of the positron emitters 7-4 sec Al and 6-6 sec 
26mAl have earlier been determined by absorption and cloud chamber 
methods (Hunt, Jones and Churchill 1954; see also Endt and Kluyver 
1954). In the present investigation the positrons have been studied by 
means of a magnetic lens spectrometer, while gamma-rays associated 
with the beta-decays were searched for with a Nal crystal. 

In the study of *°Al, thin targets of natural magnesium were bom- 
barded with 225 kev protons corresponding to the lowest *4Mg (p, y)°Al 
resonance. At this energy the yield of ?®"Al is negligible (Hunt, Jones, 
Churchill and Hancock 1954). ?6™Al was investigated by bombarding 
separated ?°Mg targets and utilizing the 391 kev capture resonance in 
25Me. The protons were accelerated in a Greinacher coupled, high tension 
generator and deflected through 90° by an analysing magnet. Several 
circular slits, 5mm in diameter, between the analysing magnet and the 
spectrometer defined the position of the beam before it hit the magnesium 
target, which was supported by a 2-5 mg/cm? Al foil. The target arrange- 
ment is shown in fig. 1. The Al foil was placed at the source position of 
the lens spectrometer, and a Geiger counter, which was situated behind 
the foil in the spectrometer axis, was used as a beta-monitor. A tungsten 
foil in front, of the monitor absorbed x-rays from the target. 

In the case of natural magnesium, the targets were prepared by 
evaporation in vacuo with subsequent oxidation in air at 450°C, and in 
the case of Mg, by spraying a dilute *Mg(NO3), solution on to a hot Al 
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foil. ®Mg was obtained from Harwell as MgO containing 3-8% *4Mg, 
and very pure %Mg targets were also prepared in the isotope separator 
of the Institute. In all cases rather stable targets were obtained, which 
were able to withstand proton beams of 10-15 amp for many hours. 

The spectrometer used was of the Siegbahn-Slatis type (Slatis and Sieg- 
bahn 1950), which utilizes the focusing properties of a U-shaped field. 
The baffle system consisted of a central baffle, which formed a ring 
aperture 8 mm wide and 20 cm in diameter, and two identical ring apertures. 
near the source end and the detector, respectively. The resolving power 
of the instrument was 3-3°%, for a point source and 3-79, for a source: 
spot 5 mm in diameter. The transmission was about 7°%. 


Ides, tl 


BAFFLE | 


Source end of lens spectrometer. 


The beta-particle detector was an anthracene crystal, which was 3 mm 
thick and 15 mm in diameter. It was mounted on the end of a 25 em 
long lucite rod in optical contact with the photo-cathode of an EMI 5311 
photomultiplier. The multiplier was shielded from the stray magnetic: 
field by means of an iron tube, and it was verified that the counting rate: 
of the detector was independent of the magnetic field when a gamma-ray 
source was placed outside the spectrometer. . 

The iron yoke of the spectrometer made it necessary to check the 
linearity of the current-momentum relationship. Measurements were 
made on several of the conversion lines of Th(B+C’-+-C’) including the: 
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2-53 Mev x-line and on the 624 kev !37Cg line. Furthermore the beta- 
rays from “K were studied and the calibration yielded for this case a 
maximum energy of 3:56-+-0-05 Mev, which is in excellent agreement with 
the accepted value 3-56 Mey (Koerts et al. 1954). 

The counting rate of the Geiger monitor showed some dependence on. 
the magnetic field strength, and correction curves were, therefore, deter- 
mined by measuring the beta-rays from 1-7 Mev =P and from 3-6 Mev 
“K. The curves for the counting rate as function of p/p», were nearly 
identical in the two cases. 

The amount of scattering present in the spectrometer was studied by 
means of the beta-rays from **P. For this purpose a thin source was 
mounted on a 100 ug/em? Zapon film. The resulting Fermi plot deviated. 
from a straight line below 700 kev. We believe that the scattering took 
place mainly at the source end of the spectrometer, where the construction 
was rather compact (see fig. 1). Another *2P source was mounted on a 
2-5 mg/cm? Al foil, but the corresponding spectrum showed no additional 
scattering. Backscattering from our target support was thus negligible 
in this energy region compared with spectrometer scattering. 

The technique used in the measurements of the Al spectra was a 
continuous irradiation and counting over a period long compared with 
the half-life, generally 10-15 minutes. The beam current to the electri- 
cally insulated target was measured by a beam integrator, by means of 
which the target yield versus beam current was controlled. 


§ 2. RESULTS 


Figure 2 shows the obtained spectrum of ?°Al, representing the average 
of 11 runs with a total of 45 000 counts. Owing to the low accelerator 
voltage used, only a machine-independent background of 3 cpm was 
observed when the detector bias corresponded to 200 kev positron 
energy. Figure 3 represents the corresponding Fermi plot from which 
the maximum energy was determined by the method of least squares 
to be 3-24-+0-03 Mev. This figure includes a 4 kev correction for energy 
loss in the target support. 

In the investigation of ?®"Al at the 391 kev resonance, we found a small 
background at lower positron energies. This was caused by 10 min *N 
formed by proton capture in carbon which deposited on the target even 
though a liquid air trap was placed immediately in front of the spectro- 
meter. Figure 4 shows the spectrum after correction for *N. It represents. 
the average of 11 runs with a total of 40000 counts. Both sprayed 
Mg targets and targets prepared in the isotope separator were used for 
the measurements, but they showed no differences, except in the positron 
yield per pcoul. The small field-independent background present at 
higher positron energies was probably produced by capture gamma-rays. 
from 26Al. A similar background was not seen in the measurements on 
25 Al, however, this can beexplained from the fact that the number of capture 
gamma rays per positron is much lower at the 225 kev *4Mg resonance 
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than at the 391 kev Mg resonance. In addition, fluorine contaminations 
might have contributed in the case of 2*Al. For calculation of the Fermi 
plot, this background was subtracted as indicated by the dotted line 
in fig. 4. The resulting Fermi plot is shown in fig. 5 and gives an end- 
point energy of 3-21-L0-03 Mev after correction for the energy loss in the 


target. 
Fig. 2 
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Positron spectrum of ?°Al. 


Fig. 3 
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It can be seen from fig. 2 and fig. 4 that the spectra differ somewhat in 
shape, but no definite conclusions can be drawn in this respect due to the 
serious scattering which took place in the spectrometer. It is, however, 
worth noting that the deviations from the straight line in the Fermi 
plots are of the same order of magnitude as was observed previously 
with *P beta-rays. 
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With a Nal(Tl) crystal and a single channel pulse-height analyser a 
search was undertaken for gamma-rays associated with the positron 
decay. Apart from annihilation radiation, which obscures the region 
below 500 kev, no gamma-rays were found. 
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Positron spectrum of ?6™A]. 
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Fermi plot of ?6™Al. 


§ 3. Discussion 


Our observed end-point energy of the positrons from Al agrees with 
the value 3:17-0:15 Mev which was obtained by absorption measurements 
(Hunt, Jones and Churchill 1954). The maximum energy of a beta- 
transition to the ground state of Mg can be derived from the Q values 
of the *4Mg(d, p)®Mg and the *4Mg(d, n)?°Al reactions. The most 
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accurate determinations (Endt et al. 1952, Goldberg 1953) give 
3:23-0:06 Mev, which is in close agreement with our measurements 
and, therefore, suggest that the observed transition leads to the ground 
state of Mg. Recent determinations of the half-life (Hunt, Jones, 
Churchill and Hancock 1954, Green, Harris and Cooper 1954) have 
given 7:62 sec and 7:20 sec, respectively. If we choose arbitrarily an 
average value of 7:41-0:21 sec, our determination of the maximum 
energy yields the ft value 3660+200 for this superallowed transition. 
It is interesting to compare this value with the estimate /t—3030, 
obtained from the unified theory of nuclear structure (Bohr and Mottel- 
son 1953). 

There is now both theoretical and experimental evidence for the 
suggestion that the two lowest lying levels in 7®Al undergo independent 
beta-decay (Moszkowski and Peaslee 1954, Simanton ef al. 1954). 
A long-lived activity with a half-life of ~10° years originates in the 
lowest 7'=0 state, which probably has spin 5+ (King and Peaslee 1953), 
while the 6:6 sec decay is assigned to the lowest 7’=1, J=0-+ state, 
leading to the ground state of 2*Mg. The latter transition thus belongs | 
to the class of 0-0(no) beta-transitions in 4n-++2 nuclei, for which the 
theory predicts a common /t value directly related to the Fermi coupling 
constant (cf. e.g. Kofoed-Hansen 1953). The most accurate experi- 
mental data in this group of beta-decays have been obtained on the decay 
of 140 (Gerhart 1954), for which the /t value 3275+-75 is reported. Fora 
calculation of the ?°"Al ft value, we take the average 6-56--0-20 see of 
the most accurate determinations of the half-life (Hunt et al., op. cit. 
and Green et al., op. cit.), which are 6-44 sec and 6-68 sec. By using 
our value of the maximum energy, one finds ff=3110+.200. 

The position of the low-lying 7’'=0 levels in ?®Al have recently been 
determined by magnetic analysis of the **Si(d, «)?®Al reaction (Browne 
1954). The two lowest of these are located 4:010+0-008 Mev and 
4-428-+-0-008 Mev over the ground state of ?°Mg and may be interpreted 
as the ground state and a 418 kev excited state in ?°Al. These levels have 
also been reported by Kluyver et al. (1954). Our measurements locate 
the (isomeric) level responsible for the 6-6 sec activity, 4:23--0-03 Mev 
over the ground state of ?®Mg, i.e. 220-+40 kev over the ground state of 
2®Al found by the above mentioned authors. This is in accordance with 
measurements of the ?7Al(y,n) threshold (R. N. H. Haslam, private 
communication quoted by Browne, op. cit.) as this process is unaccom- 
panied by positrons up to 200 kev above the neutron threshold. 


We wish to thank Professor Niels Bohr for excellent working conditions 
at the Institute, and Mr. K. O. Nielsen for having prepared the separated 
**Meg targets. 

Note added in proof—Another beta-spectroscopic determination of the 
maximum energy of 7°™Al has recently been published by R. W. Kavanagh, 


W. RR. Mills and R. Sherr (1955, Phys. Rev., 97, 248). Our work is in excellent 
agreement with their result, 3-20-L0-05 mev. 
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SUMMARY 


The standard excited state approximation for an unsaturated hydro- 
carbon molecule is used to determine the equilibrium configuration of the 
nuclei of the molecule. Each bond length of the molecule is shown to 
depend, to an approximation, on a scalar quantity, called bond order, 
which is characteristic of the bond. This concept of bond order is 
numerically equal to that of 7-bond order used by Coulson and others. 
The quantum mechanical formulation gives the form of the bond order— 
bond length relation, but its actual shape is determined from the known 
bond lengths of naphthalere and anthracene. 


§ 1. THe EQuiciprium CONFIGURATION OF A MOLECULE 


THE quantum mechanical Hamiltonian of a molecule consisting of MW 
nuclei and N electrons, in atomic units, and with spin interactions 
neglected, is 

M M 


NN 
ta Lgt 22 ig Delp t & kre, 
a>B 


ase] 


M N M 
ip eH tibet ip ets 
=1j 


a= t=1 a 


Uy 


(1.1) 


M,, Z, are the mass and charge of the ath nucleus ; ry, 7,; are inter- 
nuclear and interelectronic distances respectively ; and r,, is the distance 
between the ith electron and the «th nucleus. # may be written as 
the sum, #,+ #,, of electronic and nuclear Hamiltonians, 

MM 


M 
Hat FMS ALY Se ee 
a=1 a>fBp ‘ 


Ha 2 HE org) AR ee 
t=1 i>) 
where Mee 
Hy=—fA at 2. ES ee 


The wave function of the whole molecule may be written as Y= VW, . Y,, 
where ¥,, involves only the nuclear co-ordinates, 7,V,=EH . VY, and # isan 
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eigenvalue of the electronic energy. Then the Schroedinger equation for 
the whole molecule, 


Hee Vere Y. 
becomes 
M M 
PHY +EY,Y IS Mop ype & od} 
la=1 a=1 
meme ee ee OCT 6) 


E,, is the total molecular energy. If the terms in the brackets are neglected, 
because Y’, is a slowly varying function of the nuclear co-ordinates (Born— _ 
Oppenheimer approximation), (1.6) becomes 

MM 


M 
{- SMA PE rg Zalig+ BI eee (1.7) 
oa — a | a>B 


Hence, in the equation of the nuclear motion, the potential term is 
oe DD rg Zlig +k. eee ee (1:8) 


& is a function of the nuclear positions through all its terms. If the 
molecule is stable there exists a nuclear configuration for which & is a 
minimum, and about this configuration nuclear vibrations take place. 

In this paper the standard excited state approximation to the molecular 
orbital theory of the electronic structure of an unsaturated hydrocarbon, 
developed by Hall (1952 a), has been used to determine this equilibrium 
configuration. Arising out of this formulation, it appears that a bond in 
an unsaturated hydrocarbon can be characterized by a scalar quantity 
called bond order (so called because it is numerically the same as the 7-bond 
order used previously by Coulson and others), such that the bond lengths 
are, to an approximation, a function of this parameter alone. The form 
of the bond order—bond length relation is derived from the quantum 
mechanical formulation ; its actual shape is determined by fitting this 
curve to the very accurate values of the bond lengths of naphthalene and 
anthracene, obtained by Ahmed and Cruickshank (1952) by a critical 
revision of the very accurate x-ray determinations by Robertson and 
his co-workers (Robertson 1951). 

While it is obvious that this treatment of the problem owes much to 
previous theories, reviewed by Lennard-Jones (1951) and Coulson (1951), 
it is logically independent of them, having a different starting point in the 
standard excited state approximation for unsaturated hydrocarbons. 


§ 2. Tar STanpaRpD Excirep STATE APPROXIMATION 

The molecular orbital theory of molecular structure (Lennard-Jones 

1949, Lennard-Jones and Hall 1951 a) assumes that the total electronic 

wave function of the molecule may be written as a single determinant of 
the form 

Y=(N !)-1? det (64(41), b2(T2), P3(T3)> ---- dy(ty)}> (2.1) 

where t,, T,.... Ty represent the co-ordinates, including spin, of each of 

N electrons, and ¢,, 6... . . dy are an ortho-normal set of functions, called 
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spin orbitals, which are each functions of the co-ordinates of only one 
electron. Each spin orbital, ¢,(7,), is assumed to factorize into a space 
dependent part, the space orbital, y,(x,), and a spin dependent factor, which 
can be either of the two possible electron spin functions «(c,) and B(o,) ; in 
a molecule each space orbital, y,(x,), may give rise to either one or two of 
the spin orbitals, y,(x,) «(o,) and y,(x,) B(o,), and accordingly the space 
orbital is said to be singly or doubly occupied. The expression (2.1) is a 
valid wave function only if all the orbitals are doubly occupied, or if all the 
singly occupied orbitals have the same spin. Further, ¥, satisfies the 
variation equation (e.g. Corson, ch. x), 


8[¥, #,¥,dr=0 7 bie ect ips err 


which gives rise to the integro-differential equations for the orbitals 
(Hall 1952 a): 


d 8 
[H+V+07'+ 7" ]y,(x)— 2 CXG(X)— 2 Cy yXp(X)=9, ith tee laa 
mn Le 


d 8 
[H4+V+4.00'+.0f"lxs(x)— 2 egsya(x)—3 LesxX=0. . . (2.8%) 
be u 


In writing these equations the orbitals have been divided into two sets of 
(i) s singly occupied orbital y,, x3... . and (ii) d doubly occupied orbitals 
Xa Xp ---- Greek subscripts are used for singly occupied orbitals ; 
Greek subscripts overscored for doubly occupied orbitals when such a 
distinction is made; and italic subscripts when no such distinction is 
intended. A sufficient condition for the consistency of the orbital eqns. 
(2.3) is that the doubly and singly occupied orbitals belong to different 
irreducible representations of the symmetry group of the molecule : in the 
case of flat molecules (which are the only molecules considered here) this is 
realized if the doubly occupied orbitals are symmetrical (c-orbitals) and the 
singly occupied ones are anti-symmetrical (z-orbitals) in the plane of the 
molecule. These group theoretic conditions imply that e,;=e7,=0. 
The remaining terms in the equation are defined as usual, except that, in 
this paper, the role of singly and doubly dashed quantities is reversed in 
order that quantities arising from singly occupied orbitals are singly 
dashed, and those from doubly occupied orbitals are doubly dashed. 
Thus : 


oie | Xu(X) LH + V +00’ +00") x, (x) dx, 


7 | xa0)LH+ V+4.0'+.0"] y:(x) dx, 
7 = 1 f Lele of ue 
¥'=229,,(x), W"=2ZLGii(x), (> (Snes Parle 
it 


ll 


a) 


Im n(X,)= | Xm(%,) rij Xn(X;) dx;, 


AY) — © fan) fo) Lu X)=— E Gaul) 00: | 


H is the operator H, (1.4) without the subscript. 
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Neither the form nor the value of Y, is changed if the space orbitals 
are replaced by a new set, 


b= Ltim Xm ee epee oe 1 ha (225) 
m™ 


provided that ¢,3;=t;,=0. The orbital equations for the ;;, orbitals have 
the same form as those for y,, except that Z,,,,, replaces e,,,, where 


Un 2 rurtystre ©. + ~ - + + (2.6) 


The special set of orbitals for rah the #,,-matrix has diagonal frie are 
called molecular orbitals. These orbitals have appreciable values through- 
out the molecule; and —#,, is the vertical ionization potential of the 
electron in the corresponding orbital. A set of orbitals which transform 
into each other under the symmetry covering operations of the molecule 
are said to be equivalent. Equivalent orbitals often have a highly 
localized character, which may lead to very useful approximations. 

The ground state of an unsaturated hydrocarbon (all orbitals doubly 
occupied) can be discussed in the manner above like any other molecule ; 
but, when the electrons are considered as distributed in doubly occupied 
orbitals corresponding to bonds, inner shells and lone pairs, there is usually 
one electron per carbon atom remaining over. These cannot be distri- 
buted in doubly occupied, localized, equivalent orbitals without reducing 
the symmetry of the molecule ; hence, the straight-forward treatment of 
the ground state is more involved for an unsaturated hydrocarbon than 
for a saturated one. 

The standard excited state of an unsaturated hydrocarbon is defined as 
one in which all the doubly occupied orbitals of bonds, inner shells and 
lone pairs are c-orbitals and the one remaining electron per carbon atom is 
associated with a singly occupied z-orbital, highly localized about each 
carbon atom. The z-orbitals usually divide into sets of equivalent 
orbitals.. Molecules in the standard excited state may not have an 
equilibrium configuration, but the stable existence of the molecule in this 
state is not required by the subsequent use made of its orbitals. The 
standard excited state is an example of the general case discussed above, 
The o-orbitals have alla localized character, and therefore may be assumed, 
approximately, to be molecular orbitals. By a similar approximation the 
localized character of z-orbitals gives rise to a simplified e¢,g-matrix : in 
particular, ¢,g=0 approximately 7h a, 8 are neither equal nor are nearest 
neighbours (the labelling of the equivalent orbitals and the nuclei coincides 
conveniently). The z-molecular orbitals are 


8 
Pu= Ubu» Xv» Ga raat ee we pee 2h) 
where the y’s are the localized z-orbitals, and ¢,,, t,2,...¢,, are the 
solutions, for each of the possible values, 1, 2, 3,....s8, of u, of the set of 


equations, 
& 


Be ann) b= On ele eee a eee 2.8) 


y 
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Ey,, Eo, E33,...-, are the solutions for H, arranged in increasing 
order, of the secular equations 

| Cay — Ory EH [=0. . . . . . . . (2.9) 
Each of the eigenvectors (f,,;, t,o, ---- t,;) of (2.8) is normalized to unity, 


so that, by a general algebraic result, [t,g] is a unitary matrix ; and, as 
[e.g] is a real symmetrical matrix, [f,g] is real too. 

The standard excited state approximation to the ground state of an 
unsaturated hydrocarbon molecule makes use of the molecular orbitals 
of the same molecule, with the same nuclear configuration, but in the 
standard excited state : the nuclear configuration need not be an equili- 
brium configuration of either state. The molecular orbitals of the ground 
state are taken to be the o-molecular orbitals of the corresponding standard 
excited state, doubly occupied, together with as many of the z-molecular 
orbitals of the standard excited state as are required, also doubly occupied, 
to accommodate the remaining electrons, the z-orbitals associated with the 
smaller vertical ionization potential being chosen so as to give the most 
stable molecule. Thus, only the first half, say p=4s, of the standard 
excited state z-orbitals are used in the ground state. 


§ 3. THe ENERGY OF AN UNSATURATED HyDROCARBON MOLECULE 


The total electronic energy of an unsaturated hydrocarbon is, on the 
standard excited state approximation, 


d 8/2 
Ba 2 (Hz +E)+S (HtHt). 2. : 2 GD) 


The first summation is the energy due to the o-electrons and the second 
is the energy of the z-electrons. The theory of the approximation (Hall 
1952 a) gives, by a perturbation treatment, 


Pp 
ER=Egt Ma [py PalG' ly, da)—3(H, YalG [da ¥,)] 
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where bse 25M 
(xb, bin |G bn by = [Pax,) Ln (X;)(1/r";;) fb, (X,) hb, (X;) dx, dx; . g . (3.4) 


In (3.2) the two summations, being approximately equal, may be neglected ; 
and similarly, in (3.3) the two summations over the coulomb integrals, 
(1h, |G |v, %.), cancel approximately : so that 

d 


p p 2p 
EX (Hz+Ez)+ 2) (Hat Ban) + a 2 (hy te (Ah ab,) . -. (3.5) 
a a= a=1 y=p+1 
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These contributions to the energy are now considered separately : 
(a) The contribution of the o-electrons is 
6 = 2 (H+ Hs)+ 2 (Az-+ Ez). ae a7) 


(od (inner shells) a (o-bonds) 


=Uot EN y Vip ee eee 28 (38) 
ad 


U4, the contribution of the inner shell orbitals, may be assumed constant. 
U,,=Hj,+ Ex, where A, u are the nuclei at the ends of the o-bond whose 
orbital is #z. N,,—1, when the nuclei A, » are (bonded) nearest neigh- 
bours ; otherwise NV ,,=0. 


(6) The first contribution of the z-electrons is 
Pp 
ee eee Hi) ae ae ot 2 (3.9) 
a=1 
which, expressed in terms of the equivalent orbitals by eqn. (2.7), is 


2p 2p p 
= 3 & (2 an a (ewe ite 


A=1 p=1 \e=1 
2p 2p 
ene Pi (Cxythay)s KEL & oe (ER 
a al peeal 
where ; 
p 
Pr=2 Dd, i Ue . . . . . . . . . (3.11) 
a=1 


is called the bond order of the bond between nuclei A and n». The nuclear 
electrostatic interaction energy, i.e. the summation in eqn. (1.8), may be 
included in &, by replacing h,, by u,, as defined by Hall (1952 b). 
Further, sizeable contributions to (3.10) arose only where either A= 
(which contributions are independent of the nuclear configuration), or 
A, » are nearest neighbours ; so, approximately, 


Gx, OV one ) Nu Pig Vex; . ° . . (3.12) 


A, 
where 2V ;,=(€4,+Uj,) and V, is a constant. 


(c) The second contribution of the 7-electrons, 


ee 22. 
Gea es Hale Natty) Leas ene onl se) 
a=1 B=p+ 
becomes, using (2.7) and ee. (xauxe |G lxyxXs)= (2 |G yd), 
2p 2p 2p 
C= LD (2 ie tu) ( 2 ts, tu) (VA |G |w)oce i (3:14) 
Ayu veo \o=l B=p+1 
But (2.7) is a real unitary transformation, so p,,=p,, and 
p Py 
Steplgate ain les pO we > —«(8.15) 
p=1 b=p+1 
Hence 
2p 2p 
6x O— 4 A 2 Pay, (2 copie —Day)(VA IG |wo). >. = c - (3.16) 


Ayo vw 
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Only the coulomb integrals, i.e. those for which p»=v, A=w, in this 
expression have appreciable values (Pariser and Parr 1953, Hall 1954) ; so 


2p 
C5 e Pi (23 44 —P au) HA | |2A) 
Ae 
2p 
= Wok Dip, Wie ate eee 
Axu 


W, is a constant, being the contribution to (3.16) for A= and 
4W,,=(uA|G|uA). Further coulomb integrals not arising from nearest 
neighbour equivalent orbitals may be neglected, giving 
C$ VEW ot 2 Na Dan Wage sitentitdn sy yp Soke, Gee Oy 
Ayu 


Hence the total energy of the molecule (1.9) is, approximately, 
2p 2p 
€=65+2 DN glU Pn Vigo gala ety atoele) 
ADR 


where &)>=U,+V 5+ W, is a constant. 

It is now assumed that the functions U,,,, V;,,, Wj, are : (i) independent 
of A, »; (ii) independent of the particular conjugated hydrocarbon in 
which they arise ; (iii) functions of only the corresponding bond length : 
then, if they are written as &,, &5, &, respectively, the total energy of any 
unsaturated hydrocarbon molecule is 


je 
E= y+ he aly) dePe E(7 )+pPEsr)]. - - ~ (3.20) 
Here r, and p,, (=1, 2,...L), are the bond lengths and bond orders of the 


L bonds of the molecule. A necessary condition for an equilibrium 
configuration of the molecule is that, for /=1, 2,...L, 


0& 0 é = 
a = = €,(r,)) +915 F2(r 2. €.(r 
ar, or, (7) +2, or, 6 (r,)+P, ar, 6.(r,) 
L Op, 
Ts Pe ee hp, Scie Oe 


The values of 7;, rg,..-.7,, for which these equations are true, are the 
equilibrium bond lengths of the molecule. 


§ 4. Bonp ORDERS AND EQui-Bonb OrRpDERS 


The calculation of the bond orders of all the bonds of a molecule 
requires the solution of the equations (2.8) for some of the roots of the 
secular eqn. (2.9). If the values of e,,; are known for the molecule con- 
cerned, this is a matter of algebra involving no formal difficulty but 
considerable labour, especially if each e,; varies with molecular size. But 
generally the values of e,, are not known. Previous theories, faced with a 
sunilar problem, proceeded on the assumptions that e,,—a, for all «, 
€,j—6 for all «, 8 which are nearest neighbours and €.g—=9 for all other 
cases. This greatly simplifies the solution of the secular equation, since, 
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putting x=(a—F)/b and diving each row by 0, the determinant has x in the 
diagonal and either zero or unity elsewhere. Such secular equations can 
often be solved by special algebraic methods (Rutherford 1947, 1952). 
In any case, the problem is essentially numerical and the values of the 
bond orders are independent of the values of a and b. The bond orders 
of many molecules have been calculated under these assumptions in 
connection with previous theories. 

On the assumptions of the present discussion, bond orders calculated 
in this manner are the bond orders of a molecule all of whose bonds are 
equal, since, in the derivation of (3.21), e,,, etc. have been assumed to be 
the same function of bond length, from bond to bond. For this reason, 
these bond orders are called equi-bond orders. Their values depend on the 
bond in question but are unaffected by the actual size of the molecule, 
provided all the bonds are equal ; they depend on the connectivity of the 
o-bonds of the molecule and are essentially a topological property of the 
molecule. In an actual molecule, bond orders differ from equi-bond 
orders because of the small variation of the bond lengths ; but here it is 
assumed that these variations in bond order are small compared with the 
variation of equi-bond order from bond to bond, and that equi-bond 
orders may be used in (3.21) as good approximations to the actual bond 
orders. This is equivalent to assuming that p, depends only on the bond J, 
and is independent of small variations of the bond length from its equili- 
brium distance. To avoid confusion, it should be pointed out that this 
paper establishes a relation between the equilibrium lengths of different 
bonds and their respective bond orders, while the above assumption is 
concerned with the changes in bond order of one particular bond when all 
the bond lengths are slightly different from the equilibrium configuration. 
The numerical values of the bond orders used here are the same as the 
a-bond orders defined by Coulson (1939), so that not only may the 
previous calculations of them for several molecules be used, but their 
application to other topics, e.g. free valance and carcinogens (Coulson 
1953), remains the same. 


§ 5. Bonp OrpER—Bonp LeneTn RELATION 
Since bond order depends only on the bond, the summation in (3.21) 
vanishes, giving 


dé(r) d&(r)  ,4&s(7) _ ‘e 
Te iy ne gees ames Pe So) ete (Ot) 


as the equation relating the equilibrium length r of a bond to the bond 
order p. This is the new bond order—bond length relation. Its actual 
form depends on the functions 6 (7), & (7), &(r), which, though defined in 
the derivation given here, cannot be expressed easily in analytic form. 
The best that can be done is to adopt the semi-empirical method and find 
the form of these functions by fitting the relation (5.1) to experimentally 
determined values of bond lengths and calculated bond orders. 
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Let 7, be a length approximately equal to the mean length of the bonds 
of unsaturated hydrocarbon molecules. Putting r=r9+p, and expanding 
dé&,/dr, d€,/dr, dé /dr about r=ry as 


a =a,+6,p- 
102 ay} bap t Bi seh ic A eed 
— =A,+bsp+.... 
and neglecting terms of second order in p, (5.1) becomes 
pip eee |. : igi ae tee ae) 


This is the approximate form of the new bond order—bond length relation. 
This relation can be fitted to the experimental data by finding the values 
of a1, Gy, 43, by, by, bs, for which the expression 


S= 2 [(a,+a,p;+agp;7)+p; (1 +0ep;+-b3p7)? . . - (5.4) 
ri 


has a minimum value : p, is the calculated bond order and p, is the experi- 
mentally determined value of p for each of the n bonds used in the curve 
fitting. One of the coefficients in (5.3) can be put equal to unity without 
loss of generality ; here 6,=1. The fitting procedure is essentially the 
‘method of least squares ’ applied to (5.1). 

The bond lengths of naphthalene and anthracene, obtained by Ahmed 
and Cruickshank (1953) from earlier x-ray work by Robertson and others 
(1951), have been used in an attempt to find suitable values of the a’s 
and b’s in this manner. The bond orders of naphthalene have been given 
by Lennard-Jones and Coulson (1939) and of anthracene by Coulson, 
Dandel and Robertson (1951). It was found that this gave an absurd 
result : (5.3) had an asymptote in the middle of the range of p considered. 
The explanation of this failure is the small number of experimental points 
available, the small range of p involved and the experimental error in the 
bond length determinations. The bond lengths of many other unsaturated 
hydrocarbons are known, but with errors about four times those present 
in the naphthalene and anthracene determinations. These errors are not 
large, considered as a fraction of the bond length, but they are considerable 
when regarded as a fraction of the range of variation of the bond lengths 
involved: this is true, to some extent, even of the naphthalene and 
anthracene bond lengths. There is, therefore, little to be gained by 
increasing the number of experimental points and the range of p by 
including the bonds of other molecules, though some bonds with large 
bond order may be useful in spite of the errors. 

If the bond lengths of anthracene and naphthalene are plotted against 
the corresponding bond orders, the shape of the bond order—bond length 
relation is evident (see fig. 1): this is particularly so if the bridge bond of 
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naphthalene could be disregarded. There is very little difficulty in drawing 
a smooth curve through these points. Another point, the ‘double’ 
bond of ethylene was added to extend the range of the curve. This bond 
length, following Coulson, Dandel and Robertson, is taken as 1-344, which 
is a mean of two experimental determinations (1-331 A due to Thompson, 
1-353A due to Gallaway and Baker); the bond order is unity 
(Lennard-Jones and Coulson 1939, Lennard-Jones and Hall 1951). Such a 
smooth curve was drawn and a number of points on it, equally spaced in 
the range p=0-45 to —p=1-0, were used in the fitting procedure outlined 
above, instead of the experimental points. The result is the bond order- 
bond length relation 


r=1-40+ ( 


2:27 p2—17-38p+9-70 an 


176-08p2— 186-08p-+ 100 


which is shown in fig. 1. It is a good fit, with the possible exception of the 
bridge bond of naphthalene, and it follows closely the hand drawn curve. 


Fig. 1 
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The curve is the bond order-bond length relation given by eqn. (5.5). The 
points marked are the experimental determined bond lengths given by 
Ahmed and Cruickshank, with the estimated standard deviation marked 
off on both sides of the value. 


Though this modified procedure is not as satisfactory as a direct fitting, 
it gives a numerical relation in the form required by the theory. The 
resulting bond lengths relation is valid only over the range of experimental 
points used in the fitting. In fact, as p increases (5.5) passes through a 
minimum and then increases towards its horizontal and only asymptote. 
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§ 6. DiscussIon 


An important feature of the present theory of bond lengths is the 
manner in which the notion of bond order occurs in the derivation of the 
bond order-bond length relation. This notion arises essentially as a 
parameter in the expression (3.19) for the total energy of a molecule ; it 
is defined, eqn. (3.11), by the molecular orbitals of the standard excited 
state, and it appears in eqn. (5.1) as the determining parameter for the 
equilibrium positions of the nuclei. Nowhere is bond order related to the 
chemical notions of single and double bond: it is only after calculation 
that it can be asserted that the ethylene ‘ double ’ bond has unit 7-bond 
order, and this result is a consequence of the formally unnecessary 
factor 2 in the definition of bond order. (This factor is included here so 
that the definition of bond order agrees numerically with Coulson’s 
definition (1939) : in his theory the 2 is not without physical significance.) 
Further, it is not evident that diamond can be regarded as having zero 
bond order: it is not an unsaturated hydrocarbon, and the molecular 
orbital treatment of diamond (Hall 1952 c) has no immediate connection 
with the treatment of unsaturated hydrocarbons upon which this paper is 
based. There is, therefore, no evident reason to expect the bond order— 
bond length relation to pass through the diamond bond length for p=0. 
Moreover, the present theory can give no indication of the bond order— 
bond length relation for values of » lying outside the range of p for which 
experimentally determined bond lengths are known: in particular, it 
cannot predict the existence of stable unsaturated hydrocarbons having 
bonds with very small bond orders. The smallest bond order in a stable 
molecule known to the author is 0-45 for the ‘ single ’ bond of butadiene 
(Lennard-Jones and Coulson 1939). If unsaturated hydrocarbons with very 
small bond orders were unstable, the bond order—bond length relation 
would probably have an asymptote, p=const., for some small value of p. 
This does not occur in the present bond order—bond length relation, but 
there is no reason to suppose that it would not occur in a curve based on 
more experimental evidence. Onthe other hand, it is essential to Coulson’s 
theory (1939) that the bond length of diamond lies on the bond order— 
bond length curve for p=0. (In fact, Coulson uses 1+ p as the variable, 
adding | for the o-bond orbital.) 

From this discussion it is clear that bond order is a concept defined with 
respect to a particular theory of molecular structure and it should not be 
assumed that, taken out of this context, it can have any meaning as an 
intrinsic property of the bond (such as, e.g., length). Though numerically 
the same, the two concepts of bond order—the one defined here and the 
other used by Coulson—are essentially different, being imbedded in 
different theories. Both theories are based on the same mathematical 
problem of finding the eigenvalues and eigenvectors of the same matrix, 
but this problem arises for very different reasons in each theory. 
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The Bombardment of Chlorine by Nitrogen Ions 


By G. W. GREENLEES and A. EK. SoucH 
Physics Department, University of Birmingham* 


[Received March 17, 1955} 


§ 1.. INTRODUCTION 


AN investigation has been made of the activities produced by the 
bombardment of chlorine with the internal nitrogen beam of the 60-in. 
Nuffield cyclotron. The beam had a continuous energy distribution of 
maximum energy about 125 Mev with an intensity fall-off by a factor of 
two for each 10 Mev rise of energy above 50 Mey. A target of fused 
Li Cl was used. To avoid confusion by activities due to the lithium and 
oxygen in the target, only activities with mass numbers above chlorine 
were examined. The target was mounted in a screened probe similar 
to that described by Checkett e¢ al. (1954). A copper screening foil 
was used. 


§ 2. EXPERIMENTAL METHOD AND RESULTS 


Chemical separations were made for the elements being investigated 
and the activities were followed with geiger counters. The decay curves 
obtained were analysed and the activities present identified by their 
half-lives. A gross decay curve for the unseparated target was also 
recorded in an attempt to obtain relative yields. 

The identification of the isotopes present by their half-lives was 
unambiguous after chemical separation except in the case of scandium. 
The scandium decay curve could be analysed into three half-lives of 
~ 90 days, 2-55 days and 3-9 hr. The long lived activity was attributed 
to *6Sc (85 days) and the 2-55 day activity to “"Sc (2-4 days). The 
3-9 hr. activity could be either or both of #Se and “Se (3-92 hr. and 
4-0 hr. respectively). 

The beta-spectrum of the scandium fraction was examined using a 
simple semi-circular focusing spectrometer with path radius 10 em. This 
showed a positron activity with H,,,, 1-41-.0-07 Mev in good agreement 
with the value 1-463 Mev for “Se (Bruner and Langer 1950). From the 
spectrum it was not possible to say if any of the less energetic “3Se complex 
spectrum (/,,,, 1:18 Mev) was present. 


* Communicated by the Authors. 
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The gamma spectrum of the scandium fraction was examined using a 
Nal(Tl) crystal and photomultiplier arrangement. This showed a large 
photopeak at 1-15 Mey which decayed with a half-life of 3-9 hr. indicating 
- that most, if not all, of the 3-9 hr. activity was due to “Se (Bruner and 
Langer 1950, Cuffey 1950). A search was made for the photopeak due 
to the 0-375 Mev gamma ray from “Sc (Haskins e¢ al. 1952). A small 
peak was found on the Compton edge of the 0-511 Mev annihilation 
radiation which decayed with a half-life of about 4 hr. Energy measure- 
ments for this peak gave a value 0-385 Mev in good agreement with the 
accepted value 0-375 Mev. The abundance of 43Sc was computed to be 
Ze Of 4Sc: 

“m™Se was identified by its gamma spectrum as well as by half-life 
meaurements. A mean energy of 0-268 Mev was found for the gamma ray 
in good agreement with published data (cf. Hollander et al. 1953). 

Table 1 gives the relative yeilds of the various activities detected. The 
yields are related to an arbitrary value of 1000 for 44Sc. Those values in 
brackets are uncertain because of the small activities associated with 
the long half-lives involved. 


Table 1 

Activity Observed 7’; Yield 
44Sc¢ 3-9 hrs. 1000 
PASC 2-55 days 500 
43Sc¢ 3-9 hrs. 20 
26Sc 90 days (2400) 
38K 10 min. 250 
280 35 min. 330 
39] 52 min. 480 
2Or 39 min. 30-60 
ON 4] min. 100-500 
al 16 days (270) 
4577 3-2 hrs. 50-220 
Ca — 0 


In determining the relative yields no information on the chemical 
extraction efficiencies was available. These were assumed equal and an 
attempt was made to check the yields obtained by an analysis of the gross 
decay curve. Where a range of yield values is given, this indicates 
discrepancies in the results obtained using different extraction methods 
and the analysis of the gross decay curve. In general the accuracy of 
the. yields is not better than a factor of two. 


§ 3. Discussion 
Two possible mechanisms have been suggested. to account for the yields 
of heavy ion reactions: (1) via a compound nucleus (Cohen et al. 1954) 


and (2) by a direct interaction of the ‘ stripping © or ‘ buckshot type 
(Chackett et al. 1954). The present results give a qualitative fit with 


SER. 7, VOL. 46, NO. 377.—JUNE 1955 3A 


684 Correspondence 


either mechanism. More detailed experiments with monoenergetic 
beams are required to determine the relative importance of the two 
processes. 

Using the heavy ion beams available from a 60-in. cyclotron, semi- 
classical wavelength considerations suggest that high values of incident 
angular momentum would be involved in the reactions. This is borne 
out by the scandium results where “4Sc and 44"Sc are produced in approxi- 
mately equal amounts despite the high spin difference between the two 
levels of *4Sc. 

Heavy ion bombardments may prove a useful means of exciting high 
spin isomeric states not readily produced in normal reactions. By such 
means it may be possible to produce as yet undetected isomeric states 
which are expected from shell structure considerations. 
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LXXVIII. REVIEWS OF BOOKS 


X-Ray Diffraction by Polycrystalline Materials. (London: Institute of 

Physics.) [Pp. 725.] Price 63s. 

Tats work, sponsored by the x-ray Analysis Group of the Institute of Physies 
is issued as one of the volumes of the series Physics in Industry published by 
the Institute. Although it deals with only one branch of crystal analysis, 
the study of polycrystalline materials, the volume of 32 chapters and an 
appendix of useful data has been compiled by 28 British and 3 American authors, 
each qualified to write with authority on one or more of the many aspects of the 
subject which are covered. 

Space does not permit of a detailed review of this comprehensive and wholly 
admirable work. Its general nature can be judged from the fact that ten of 
the authors work in industrial laboratories, nine in Universities, five in Jolleges 
of Technology, four in Government Establishments, two in Research Associa- 
tions and one in the Royal Institution. 

The last nine chapters of the book are in the form of essays showing how the 
methods described in the earlier chapters can be applied in many different 
scientific fields. The technical editors Messrs. Peiser, Rooksby and A. J. C. 
Wilson are to be congratulated on the satisfactory completion of a formidable 
task, presented in the customary high quality of publications of the Institute 
of Physics. It is clear that it will remain a standard work of reference for many 
years to come. AS MI" 


Advances in Electronics and Electron Physics, Vol. VI, 1954. [Pp. xi+538.] 
(Academic Press, New York.) $11.80. 


THis volume continues the series Advances in Electronics under an expanded 
title which is thought to describe its contents more accurately. Certainly the 
net is cast wide, to include ** Metallic Conduction at High Frequencies and Low 
Temperatures” (A. B. Pippard); “‘ Relaxation Processes in Ferromagnetism”’ 
(Elihu Abrahams); “‘ Physical Properties of Ferrites” (J. Smit and H. P. J. Wijn); 
‘Space Charge Limited Currents” (H. F. Ivey); ‘A Comparison of 
Analogous Semiconductor and Gaseous Electronics Devices ”’ (W. M. Webster) ; 
“The Electron Microscope—A Review” (M. E. Haine); “Travelling-Wave 
Tubes ” (R. G. E. Hutter); “ Paramagnetism ” (J. Van Den Handel). Indeed 
the variety of topics may be so great that few individuals will care to pay 
this price just for the one or two articles in which they are interested. The 
production is lavish, including name and author indices to the whole volume ; 
one wonders whether it might not be better to publish the articles separately 
and cheaply, perhaps in paper covers, the whole being bound together only in 
a library edition. Jawa Ze 


Superfluids : Macroscopic Theory of Superfluid Helium, Vol. I. By F. Lonnon. 
[Pp. xvi+217.] (Chapman & Hall.) Price 64s. 
Art the time of his death in March 1954 Fritz London was still at the height of 
his intellectual powers. This posthumous second volume of Superfluids is a 
lasting memorial to the importance of his contributions to, and the depth of 
his understanding of, the problem of liquid helium, which he first began to, 
investigate about 1936 and which remained his chief interest until the end. 
Although both volumes purport in their titles to be concerned with macro- 
scopic theory, of superconductivity and superfluid helium respectively, they 
do in fact differ considerably in their approach. The former was indeed almost 
wholly taken up with the macroscopic description of superconductors without 
more than the most general observations on the possible microscopic explanation. 
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of the phenomena. The latter, on the other hand, devotes considerable space 
to genuinely microscopic matters, such as the atomic structure of the liquid 
and the extent to which it may resemble a degenerate Bose—Einstein gas. 
This difference reflects the different approaches which have been made to the 
two problems: repeated failure to achieve a microscopic theory of super- 
conductivity has encouraged the development of a phenomenological picture, 
while with liquid helium the phenomenological treatment has tended to arise 
from more or less successful attempts at formulating a microscopic theory. 
It is only natural that the distinction between microscopic and macroscopic 
should have become somwehat blurred. Nevertheless the emphasis of the second 
volume is still predominantly macroscopic in so far as the most important 
single topic treated is the two-fluid model of helium, originally suggested by 
Tisza. The model is presented with historical truth as springing from London’s 
analogy of the A-point with the temperature at which a Bose-Einstein gas 
becomes degenerate. At the same time the alternative development by Landau 
and later workers, arising out of the quantum theory of a fluid continuum, is 
dealt with in lucid detail. There seems to be little doubt now that the two- 
fluid model gives an adequate description of superfluid helium, and there is 
general agreement on the form of the hydrodynamical equations so far as the 
linear terms are concerned. It is in the non-linear, particularly dissipative, 
effects that discrepancies between theory and experiment become apparent. 
These matters are given a full treatment, which constitutes a very valuable 
section of the book. There is also a detailed discussion of the preperties of 
liquid helium—3 and mixtures of the two isotopes. 

So long as a complete solution of the problem of liquid helium is wanting, 
it would be rash to assert that all the conclusions of the book will stand the 
test of time. Nevertheless so much of the theory developed is so firmly based 
that it is most unlikely to be modified except in matters of detail, and this last 
book of London’s seems destined to occupy for many years a place of honour 
in the literature of low temperature physics. 
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The Language of Taxonomy—an application of symbolic logic to the study of 
classificatory systems. By Joun R. Grea, Associate Professor of Zoology, 
Columbia University. [Pp. ix-+-70.] Columbia University Press (London : 
Geoffrey Cumberlege). Price 20s. 


Mathematics and Plausible Reasoning. Volume I: Induction and Analogy in 
Mathematics ; Volume IT : Patterns of Plausible Inference. By G. Ponya. 
[Vol. I, pp. xvi+280; Vol. Il, pp. x+190.] Princeton University Press 
(London: Geoffrey Cumberlege). Price: Vol. I, 42s., Vol. II, 35s. : 
Vols. [ and II, 70s. . 
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[Pp. x+451.] Princeton University Press (London : Geoffrey Cumberlege). 
Price 63s. 

Contributions to the Theory of Partial Differential Equations. Edited by 
L. Bers, 8S. Bocuner and F. Joun. Papers read at the Conference on 
Partial Differential Equations sponsored by the National Academy of 
Sciences—National Research Council, October 1952. [Pp. 257.] Princeton 
University Press (London : Geoffrey Cumberlege). Price 32s. 
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